      module mo_lu_factor
      private
      public :: lu_fac
      contains
      subroutine lu_fac01( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(1) = 1._r8 / lu(1)
         lu(2) = 1._r8 / lu(2)
         lu(3) = 1._r8 / lu(3)
         lu(4) = 1._r8 / lu(4)
         lu(5) = 1._r8 / lu(5)
         lu(6) = lu(6) * lu(5)
         lu(7) = 1._r8 / lu(7)
         lu(8) = 1._r8 / lu(8)
         lu(9) = lu(9) * lu(8)
         lu(10) = 1._r8 / lu(10)
         lu(11) = 1._r8 / lu(11)
         lu(12) = 1._r8 / lu(12)
         lu(13) = 1._r8 / lu(13)
         lu(14) = 1._r8 / lu(14)
         lu(15) = 1._r8 / lu(15)
         lu(16) = 1._r8 / lu(16)
         lu(17) = 1._r8 / lu(17)
         lu(18) = 1._r8 / lu(18)
         lu(19) = 1._r8 / lu(19)
         lu(20) = lu(20) * lu(19)
         lu(21) = lu(21) * lu(19)
         lu(751) = lu(751) - lu(20) * lu(742)
         lu(760) = lu(760) - lu(21) * lu(742)
         lu(22) = 1._r8 / lu(22)
         lu(23) = lu(23) * lu(22)
         lu(24) = lu(24) * lu(22)
         lu(1008) = lu(1008) - lu(23) * lu(993)
         lu(1078) = lu(1078) - lu(24) * lu(993)
         lu(25) = 1._r8 / lu(25)
         lu(26) = lu(26) * lu(25)
         lu(27) = lu(27) * lu(25)
         lu(1040) = lu(1040) - lu(26) * lu(994)
         lu(1078) = lu(1078) - lu(27) * lu(994)
         lu(28) = 1._r8 / lu(28)
         lu(29) = lu(29) * lu(28)
         lu(30) = lu(30) * lu(28)
         lu(1077) = lu(1077) - lu(29) * lu(995)
         lu(1078) = lu(1078) - lu(30) * lu(995)
         lu(31) = 1._r8 / lu(31)
         lu(32) = lu(32) * lu(31)
         lu(33) = lu(33) * lu(31)
         lu(1073) = lu(1073) - lu(32) * lu(996)
         lu(1078) = lu(1078) - lu(33) * lu(996)
         lu(34) = 1._r8 / lu(34)
         lu(35) = lu(35) * lu(34)
         lu(36) = lu(36) * lu(34)
         lu(37) = lu(37) * lu(34)
         lu(1041) = lu(1041) - lu(35) * lu(997)
         lu(1077) = lu(1077) - lu(36) * lu(997)
         lu(1078) = lu(1078) - lu(37) * lu(997)
         lu(38) = 1._r8 / lu(38)
         lu(39) = lu(39) * lu(38)
         lu(40) = lu(40) * lu(38)
         lu(41) = lu(41) * lu(38)
         lu(1046) = lu(1046) - lu(39) * lu(998)
         lu(1077) = lu(1077) - lu(40) * lu(998)
         lu(1078) = lu(1078) - lu(41) * lu(998)
         lu(43) = 1._r8 / lu(43)
         lu(44) = lu(44) * lu(43)
         lu(87) = lu(87) - lu(44) * lu(82)
         lu(740) = lu(740) - lu(44) * lu(708)
         lu(1078) = lu(1078) - lu(44) * lu(999)
         lu(45) = 1._r8 / lu(45)
         lu(46) = lu(46) * lu(45)
         lu(47) = lu(47) * lu(45)
         lu(48) = lu(48) * lu(45)
         lu(49) = lu(49) * lu(45)
         lu(1001) = lu(1001) - lu(46) * lu(1000)
         lu(1028) = lu(1028) - lu(47) * lu(1000)
         lu(1073) = lu(1073) - lu(48) * lu(1000)
         lu(1078) = lu(1078) - lu(49) * lu(1000)
         lu(50) = 1._r8 / lu(50)
         lu(51) = lu(51) * lu(50)
         lu(52) = lu(52) * lu(50)
         lu(1003) = lu(1003) - lu(51) * lu(1001)
         lu(1078) = lu(1078) - lu(52) * lu(1001)
      end subroutine lu_fac01
      subroutine lu_fac02( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(53) = 1._r8 / lu(53)
         lu(54) = lu(54) * lu(53)
         lu(55) = lu(55) * lu(53)
         lu(751) = lu(751) - lu(54) * lu(743)
         lu(755) = lu(755) - lu(55) * lu(743)
         lu(935) = lu(935) - lu(54) * lu(927)
         lu(939) = lu(939) - lu(55) * lu(927)
         lu(56) = 1._r8 / lu(56)
         lu(57) = lu(57) * lu(56)
         lu(58) = lu(58) * lu(56)
         lu(59) = lu(59) * lu(56)
         lu(60) = lu(60) * lu(56)
         lu(1019) = lu(1019) - lu(57) * lu(1002)
         lu(1027) = lu(1027) - lu(58) * lu(1002)
         lu(1073) = lu(1073) - lu(59) * lu(1002)
         lu(1078) = lu(1078) - lu(60) * lu(1002)
         lu(61) = 1._r8 / lu(61)
         lu(62) = lu(62) * lu(61)
         lu(426) = lu(426) - lu(62) * lu(419)
         lu(444) = lu(444) - lu(62) * lu(433)
         lu(706) = lu(706) - lu(62) * lu(688)
         lu(760) = lu(760) - lu(62) * lu(744)
         lu(1098) = lu(1098) - lu(62) * lu(1081)
         lu(63) = 1._r8 / lu(63)
         lu(64) = lu(64) * lu(63)
         lu(65) = lu(65) * lu(63)
         lu(66) = lu(66) * lu(63)
         lu(148) = - lu(64) * lu(146)
         lu(150) = - lu(65) * lu(146)
         lu(152) = - lu(66) * lu(146)
         lu(810) = lu(810) - lu(64) * lu(798)
         lu(837) = lu(837) - lu(65) * lu(798)
         lu(844) = lu(844) - lu(66) * lu(798)
         lu(67) = 1._r8 / lu(67)
         lu(68) = lu(68) * lu(67)
         lu(69) = lu(69) * lu(67)
         lu(70) = lu(70) * lu(67)
         lu(951) = lu(951) - lu(68) * lu(946)
         lu(969) = lu(969) - lu(69) * lu(946)
         lu(972) = lu(972) - lu(70) * lu(946)
         lu(1026) = - lu(68) * lu(1003)
         lu(1073) = lu(1073) - lu(69) * lu(1003)
         lu(1076) = lu(1076) - lu(70) * lu(1003)
         lu(71) = 1._r8 / lu(71)
         lu(72) = lu(72) * lu(71)
         lu(73) = lu(73) * lu(71)
         lu(74) = lu(74) * lu(71)
         lu(836) = lu(836) - lu(72) * lu(799)
         lu(838) = lu(838) - lu(73) * lu(799)
         lu(840) = lu(840) - lu(74) * lu(799)
         lu(968) = lu(968) - lu(72) * lu(947)
         lu(970) = lu(970) - lu(73) * lu(947)
         lu(972) = lu(972) - lu(74) * lu(947)
         lu(75) = 1._r8 / lu(75)
         lu(76) = lu(76) * lu(75)
         lu(77) = lu(77) * lu(75)
         lu(78) = lu(78) * lu(75)
         lu(348) = lu(348) - lu(76) * lu(347)
         lu(353) = lu(353) - lu(77) * lu(347)
         lu(356) = lu(356) - lu(78) * lu(347)
         lu(1043) = - lu(76) * lu(1004)
         lu(1073) = lu(1073) - lu(77) * lu(1004)
         lu(1078) = lu(1078) - lu(78) * lu(1004)
         lu(79) = 1._r8 / lu(79)
         lu(80) = lu(80) * lu(79)
         lu(81) = lu(81) * lu(79)
         lu(371) = lu(371) - lu(80) * lu(369)
         lu(377) = - lu(81) * lu(369)
         lu(695) = - lu(80) * lu(689)
         lu(706) = lu(706) - lu(81) * lu(689)
         lu(748) = lu(748) - lu(80) * lu(745)
         lu(760) = lu(760) - lu(81) * lu(745)
         lu(932) = lu(932) - lu(80) * lu(928)
         lu(944) = lu(944) - lu(81) * lu(928)
         lu(83) = 1._r8 / lu(83)
         lu(84) = lu(84) * lu(83)
         lu(85) = lu(85) * lu(83)
         lu(86) = lu(86) * lu(83)
         lu(87) = lu(87) * lu(83)
         lu(717) = lu(717) - lu(84) * lu(709)
         lu(733) = lu(733) - lu(85) * lu(709)
         lu(736) = lu(736) - lu(86) * lu(709)
         lu(740) = lu(740) - lu(87) * lu(709)
         lu(1047) = lu(1047) - lu(84) * lu(1005)
         lu(1069) = lu(1069) - lu(85) * lu(1005)
         lu(1073) = lu(1073) - lu(86) * lu(1005)
         lu(1078) = lu(1078) - lu(87) * lu(1005)
         lu(88) = 1._r8 / lu(88)
         lu(89) = lu(89) * lu(88)
         lu(90) = lu(90) * lu(88)
         lu(91) = lu(91) * lu(88)
         lu(334) = lu(334) - lu(89) * lu(330)
         lu(337) = lu(337) - lu(90) * lu(330)
         lu(339) = - lu(91) * lu(330)
         lu(667) = lu(667) - lu(89) * lu(656)
         lu(682) = lu(682) - lu(90) * lu(656)
         lu(686) = - lu(91) * lu(656)
         lu(1053) = lu(1053) - lu(89) * lu(1006)
         lu(1073) = lu(1073) - lu(90) * lu(1006)
         lu(1078) = lu(1078) - lu(91) * lu(1006)
         lu(92) = 1._r8 / lu(92)
         lu(93) = lu(93) * lu(92)
         lu(94) = lu(94) * lu(92)
         lu(95) = lu(95) * lu(92)
         lu(96) = lu(96) * lu(92)
         lu(97) = lu(97) * lu(92)
         lu(717) = lu(717) - lu(93) * lu(710)
         lu(733) = lu(733) - lu(94) * lu(710)
         lu(735) = lu(735) - lu(95) * lu(710)
         lu(737) = lu(737) - lu(96) * lu(710)
         lu(738) = lu(738) - lu(97) * lu(710)
         lu(956) = lu(956) - lu(93) * lu(948)
         lu(965) = lu(965) - lu(94) * lu(948)
         lu(968) = lu(968) - lu(95) * lu(948)
         lu(970) = lu(970) - lu(96) * lu(948)
         lu(972) = lu(972) - lu(97) * lu(948)
      end subroutine lu_fac02
      subroutine lu_fac03( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(98) = 1._r8 / lu(98)
         lu(99) = lu(99) * lu(98)
         lu(100) = lu(100) * lu(98)
         lu(101) = lu(101) * lu(98)
         lu(102) = lu(102) * lu(98)
         lu(333) = lu(333) - lu(99) * lu(331)
         lu(334) = lu(334) - lu(100) * lu(331)
         lu(337) = lu(337) - lu(101) * lu(331)
         lu(339) = lu(339) - lu(102) * lu(331)
         lu(869) = lu(869) - lu(99) * lu(845)
         lu(877) = - lu(100) * lu(845)
         lu(897) = lu(897) - lu(101) * lu(845)
         lu(902) = lu(902) - lu(102) * lu(845)
         lu(1041) = lu(1041) - lu(99) * lu(1007)
         lu(1053) = lu(1053) - lu(100) * lu(1007)
         lu(1073) = lu(1073) - lu(101) * lu(1007)
         lu(1078) = lu(1078) - lu(102) * lu(1007)
         lu(103) = 1._r8 / lu(103)
         lu(104) = lu(104) * lu(103)
         lu(105) = lu(105) * lu(103)
         lu(106) = lu(106) * lu(103)
         lu(107) = lu(107) * lu(103)
         lu(108) = lu(108) * lu(103)
         lu(109) = lu(109) * lu(103)
         lu(814) = lu(814) - lu(104) * lu(800)
         lu(818) = lu(818) - lu(105) * lu(800)
         lu(836) = lu(836) - lu(106) * lu(800)
         lu(837) = lu(837) - lu(107) * lu(800)
         lu(840) = lu(840) - lu(108) * lu(800)
         lu(844) = lu(844) - lu(109) * lu(800)
         lu(1042) = lu(1042) - lu(104) * lu(1008)
         lu(1053) = lu(1053) - lu(105) * lu(1008)
         lu(1072) = - lu(106) * lu(1008)
         lu(1073) = lu(1073) - lu(107) * lu(1008)
         lu(1076) = lu(1076) - lu(108) * lu(1008)
         lu(1080) = lu(1080) - lu(109) * lu(1008)
         lu(110) = 1._r8 / lu(110)
         lu(111) = lu(111) * lu(110)
         lu(112) = lu(112) * lu(110)
         lu(113) = lu(113) * lu(110)
         lu(114) = lu(114) * lu(110)
         lu(308) = lu(308) - lu(111) * lu(307)
         lu(309) = lu(309) - lu(112) * lu(307)
         lu(310) = lu(310) - lu(113) * lu(307)
         lu(314) = - lu(114) * lu(307)
         lu(867) = lu(867) - lu(111) * lu(846)
         lu(877) = lu(877) - lu(112) * lu(846)
         lu(890) = lu(890) - lu(113) * lu(846)
         lu(902) = lu(902) - lu(114) * lu(846)
         lu(1039) = lu(1039) - lu(111) * lu(1009)
         lu(1053) = lu(1053) - lu(112) * lu(1009)
         lu(1066) = lu(1066) - lu(113) * lu(1009)
         lu(1078) = lu(1078) - lu(114) * lu(1009)
         lu(115) = 1._r8 / lu(115)
         lu(116) = lu(116) * lu(115)
         lu(117) = lu(117) * lu(115)
         lu(118) = lu(118) * lu(115)
         lu(119) = lu(119) * lu(115)
         lu(584) = lu(584) - lu(116) * lu(577)
         lu(585) = - lu(117) * lu(577)
         lu(590) = lu(590) - lu(118) * lu(577)
         lu(593) = - lu(119) * lu(577)
         lu(887) = lu(887) - lu(116) * lu(847)
         lu(889) = lu(889) - lu(117) * lu(847)
         lu(897) = lu(897) - lu(118) * lu(847)
         lu(902) = lu(902) - lu(119) * lu(847)
         lu(1063) = lu(1063) - lu(116) * lu(1010)
         lu(1065) = lu(1065) - lu(117) * lu(1010)
         lu(1073) = lu(1073) - lu(118) * lu(1010)
         lu(1078) = lu(1078) - lu(119) * lu(1010)
         lu(120) = 1._r8 / lu(120)
         lu(121) = lu(121) * lu(120)
         lu(122) = lu(122) * lu(120)
         lu(123) = lu(123) * lu(120)
         lu(487) = - lu(121) * lu(483)
         lu(493) = - lu(122) * lu(483)
         lu(494) = - lu(123) * lu(483)
         lu(551) = lu(551) - lu(121) * lu(545)
         lu(558) = - lu(122) * lu(545)
         lu(559) = - lu(123) * lu(545)
         lu(885) = lu(885) - lu(121) * lu(848)
         lu(901) = lu(901) - lu(122) * lu(848)
         lu(902) = lu(902) - lu(123) * lu(848)
         lu(1061) = lu(1061) - lu(121) * lu(1011)
         lu(1077) = lu(1077) - lu(122) * lu(1011)
         lu(1078) = lu(1078) - lu(123) * lu(1011)
         lu(124) = 1._r8 / lu(124)
         lu(125) = lu(125) * lu(124)
         lu(126) = lu(126) * lu(124)
         lu(127) = lu(127) * lu(124)
         lu(128) = lu(128) * lu(124)
         lu(129) = lu(129) * lu(124)
         lu(130) = lu(130) * lu(124)
         lu(763) = - lu(125) * lu(762)
         lu(764) = lu(764) - lu(126) * lu(762)
         lu(788) = lu(788) - lu(127) * lu(762)
         lu(790) = lu(790) - lu(128) * lu(762)
         lu(795) = lu(795) - lu(129) * lu(762)
         lu(797) = lu(797) - lu(130) * lu(762)
         lu(1016) = lu(1016) - lu(125) * lu(1012)
         lu(1019) = lu(1019) - lu(126) * lu(1012)
         lu(1071) = lu(1071) - lu(127) * lu(1012)
         lu(1073) = lu(1073) - lu(128) * lu(1012)
         lu(1078) = lu(1078) - lu(129) * lu(1012)
         lu(1080) = lu(1080) - lu(130) * lu(1012)
         lu(131) = 1._r8 / lu(131)
         lu(132) = lu(132) * lu(131)
         lu(133) = lu(133) * lu(131)
         lu(413) = - lu(132) * lu(411)
         lu(417) = lu(417) - lu(133) * lu(411)
         lu(601) = - lu(132) * lu(595)
         lu(614) = - lu(133) * lu(595)
         lu(641) = lu(641) - lu(132) * lu(632)
         lu(654) = - lu(133) * lu(632)
         lu(672) = lu(672) - lu(132) * lu(657)
         lu(686) = lu(686) - lu(133) * lu(657)
         lu(726) = lu(726) - lu(132) * lu(711)
         lu(740) = lu(740) - lu(133) * lu(711)
         lu(826) = lu(826) - lu(132) * lu(801)
         lu(842) = lu(842) - lu(133) * lu(801)
         lu(1061) = lu(1061) - lu(132) * lu(1013)
         lu(1078) = lu(1078) - lu(133) * lu(1013)
      end subroutine lu_fac03
      subroutine lu_fac04( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(134) = 1._r8 / lu(134)
         lu(135) = lu(135) * lu(134)
         lu(136) = lu(136) * lu(134)
         lu(137) = lu(137) * lu(134)
         lu(138) = lu(138) * lu(134)
         lu(139) = lu(139) * lu(134)
         lu(394) = lu(394) - lu(135) * lu(393)
         lu(395) = lu(395) - lu(136) * lu(393)
         lu(401) = lu(401) - lu(137) * lu(393)
         lu(403) = - lu(138) * lu(393)
         lu(404) = - lu(139) * lu(393)
         lu(870) = - lu(135) * lu(849)
         lu(873) = lu(873) - lu(136) * lu(849)
         lu(897) = lu(897) - lu(137) * lu(849)
         lu(901) = lu(901) - lu(138) * lu(849)
         lu(902) = lu(902) - lu(139) * lu(849)
         lu(1042) = lu(1042) - lu(135) * lu(1014)
         lu(1046) = lu(1046) - lu(136) * lu(1014)
         lu(1073) = lu(1073) - lu(137) * lu(1014)
         lu(1077) = lu(1077) - lu(138) * lu(1014)
         lu(1078) = lu(1078) - lu(139) * lu(1014)
         lu(140) = 1._r8 / lu(140)
         lu(141) = lu(141) * lu(140)
         lu(142) = lu(142) * lu(140)
         lu(143) = lu(143) * lu(140)
         lu(144) = lu(144) * lu(140)
         lu(145) = lu(145) * lu(140)
         lu(472) = lu(472) - lu(141) * lu(470)
         lu(475) = lu(475) - lu(142) * lu(470)
         lu(480) = - lu(143) * lu(470)
         lu(481) = - lu(144) * lu(470)
         lu(482) = lu(482) - lu(145) * lu(470)
         lu(878) = lu(878) - lu(141) * lu(850)
         lu(890) = lu(890) - lu(142) * lu(850)
         lu(901) = lu(901) - lu(143) * lu(850)
         lu(902) = lu(902) - lu(144) * lu(850)
         lu(904) = lu(904) - lu(145) * lu(850)
         lu(1054) = lu(1054) - lu(141) * lu(1015)
         lu(1066) = lu(1066) - lu(142) * lu(1015)
         lu(1077) = lu(1077) - lu(143) * lu(1015)
         lu(1078) = lu(1078) - lu(144) * lu(1015)
         lu(1080) = lu(1080) - lu(145) * lu(1015)
         lu(147) = 1._r8 / lu(147)
         lu(148) = lu(148) * lu(147)
         lu(149) = lu(149) * lu(147)
         lu(150) = lu(150) * lu(147)
         lu(151) = lu(151) * lu(147)
         lu(152) = lu(152) * lu(147)
         lu(765) = - lu(148) * lu(763)
         lu(789) = lu(789) - lu(149) * lu(763)
         lu(790) = lu(790) - lu(150) * lu(763)
         lu(793) = lu(793) - lu(151) * lu(763)
         lu(797) = lu(797) - lu(152) * lu(763)
         lu(810) = lu(810) - lu(148) * lu(802)
         lu(836) = lu(836) - lu(149) * lu(802)
         lu(837) = lu(837) - lu(150) * lu(802)
         lu(840) = lu(840) - lu(151) * lu(802)
         lu(844) = lu(844) - lu(152) * lu(802)
         lu(1034) = lu(1034) - lu(148) * lu(1016)
         lu(1072) = lu(1072) - lu(149) * lu(1016)
         lu(1073) = lu(1073) - lu(150) * lu(1016)
         lu(1076) = lu(1076) - lu(151) * lu(1016)
         lu(1080) = lu(1080) - lu(152) * lu(1016)
         lu(153) = 1._r8 / lu(153)
         lu(154) = lu(154) * lu(153)
         lu(155) = lu(155) * lu(153)
         lu(156) = lu(156) * lu(153)
         lu(157) = lu(157) * lu(153)
         lu(158) = lu(158) * lu(153)
         lu(221) = lu(221) - lu(154) * lu(220)
         lu(222) = lu(222) - lu(155) * lu(220)
         lu(223) = lu(223) - lu(156) * lu(220)
         lu(224) = lu(224) - lu(157) * lu(220)
         lu(229) = - lu(158) * lu(220)
         lu(858) = - lu(154) * lu(851)
         lu(859) = - lu(155) * lu(851)
         lu(860) = lu(860) - lu(156) * lu(851)
         lu(884) = - lu(157) * lu(851)
         lu(902) = lu(902) - lu(158) * lu(851)
         lu(1026) = lu(1026) - lu(154) * lu(1017)
         lu(1027) = lu(1027) - lu(155) * lu(1017)
         lu(1028) = lu(1028) - lu(156) * lu(1017)
         lu(1060) = lu(1060) - lu(157) * lu(1017)
         lu(1078) = lu(1078) - lu(158) * lu(1017)
         lu(159) = 1._r8 / lu(159)
         lu(160) = lu(160) * lu(159)
         lu(161) = lu(161) * lu(159)
         lu(162) = lu(162) * lu(159)
         lu(323) = - lu(160) * lu(315)
         lu(327) = lu(327) - lu(161) * lu(315)
         lu(328) = - lu(162) * lu(315)
         lu(451) = - lu(160) * lu(445)
         lu(457) = - lu(161) * lu(445)
         lu(459) = lu(459) - lu(162) * lu(445)
         lu(725) = lu(725) - lu(160) * lu(712)
         lu(738) = lu(738) - lu(161) * lu(712)
         lu(740) = lu(740) - lu(162) * lu(712)
         lu(825) = lu(825) - lu(160) * lu(803)
         lu(840) = lu(840) - lu(161) * lu(803)
         lu(842) = lu(842) - lu(162) * lu(803)
         lu(1060) = lu(1060) - lu(160) * lu(1018)
         lu(1076) = lu(1076) - lu(161) * lu(1018)
         lu(1078) = lu(1078) - lu(162) * lu(1018)
         lu(163) = 1._r8 / lu(163)
         lu(164) = lu(164) * lu(163)
         lu(165) = lu(165) * lu(163)
         lu(166) = lu(166) * lu(163)
         lu(201) = lu(201) - lu(164) * lu(198)
         lu(203) = - lu(165) * lu(198)
         lu(204) = - lu(166) * lu(198)
         lu(790) = lu(790) - lu(164) * lu(764)
         lu(794) = - lu(165) * lu(764)
         lu(795) = lu(795) - lu(166) * lu(764)
         lu(837) = lu(837) - lu(164) * lu(804)
         lu(841) = - lu(165) * lu(804)
         lu(842) = lu(842) - lu(166) * lu(804)
         lu(897) = lu(897) - lu(164) * lu(852)
         lu(901) = lu(901) - lu(165) * lu(852)
         lu(902) = lu(902) - lu(166) * lu(852)
         lu(1073) = lu(1073) - lu(164) * lu(1019)
         lu(1077) = lu(1077) - lu(165) * lu(1019)
         lu(1078) = lu(1078) - lu(166) * lu(1019)
         lu(167) = 1._r8 / lu(167)
         lu(168) = lu(168) * lu(167)
         lu(169) = lu(169) * lu(167)
         lu(170) = lu(170) * lu(167)
         lu(171) = lu(171) * lu(167)
         lu(172) = lu(172) * lu(167)
         lu(893) = lu(893) - lu(168) * lu(853)
         lu(897) = lu(897) - lu(169) * lu(853)
         lu(900) = lu(900) - lu(170) * lu(853)
         lu(901) = lu(901) - lu(171) * lu(853)
         lu(902) = lu(902) - lu(172) * lu(853)
         lu(965) = lu(965) - lu(168) * lu(949)
         lu(969) = lu(969) - lu(169) * lu(949)
         lu(972) = lu(972) - lu(170) * lu(949)
         lu(973) = - lu(171) * lu(949)
         lu(974) = lu(974) - lu(172) * lu(949)
         lu(1069) = lu(1069) - lu(168) * lu(1020)
         lu(1073) = lu(1073) - lu(169) * lu(1020)
         lu(1076) = lu(1076) - lu(170) * lu(1020)
         lu(1077) = lu(1077) - lu(171) * lu(1020)
         lu(1078) = lu(1078) - lu(172) * lu(1020)
      end subroutine lu_fac04
      subroutine lu_fac05( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(173) = 1._r8 / lu(173)
         lu(174) = lu(174) * lu(173)
         lu(175) = lu(175) * lu(173)
         lu(176) = lu(176) * lu(173)
         lu(177) = lu(177) * lu(173)
         lu(178) = lu(178) * lu(173)
         lu(666) = - lu(174) * lu(658)
         lu(678) = lu(678) - lu(175) * lu(658)
         lu(685) = - lu(176) * lu(658)
         lu(686) = lu(686) - lu(177) * lu(658)
         lu(687) = lu(687) - lu(178) * lu(658)
         lu(876) = lu(876) - lu(174) * lu(854)
         lu(891) = lu(891) - lu(175) * lu(854)
         lu(901) = lu(901) - lu(176) * lu(854)
         lu(902) = lu(902) - lu(177) * lu(854)
         lu(904) = lu(904) - lu(178) * lu(854)
         lu(1050) = lu(1050) - lu(174) * lu(1021)
         lu(1067) = lu(1067) - lu(175) * lu(1021)
         lu(1077) = lu(1077) - lu(176) * lu(1021)
         lu(1078) = lu(1078) - lu(177) * lu(1021)
         lu(1080) = lu(1080) - lu(178) * lu(1021)
         lu(179) = 1._r8 / lu(179)
         lu(180) = lu(180) * lu(179)
         lu(181) = lu(181) * lu(179)
         lu(182) = lu(182) * lu(179)
         lu(183) = lu(183) * lu(179)
         lu(298) = - lu(180) * lu(296)
         lu(300) = - lu(181) * lu(296)
         lu(303) = - lu(182) * lu(296)
         lu(305) = lu(305) - lu(183) * lu(296)
         lu(318) = - lu(180) * lu(316)
         lu(320) = - lu(181) * lu(316)
         lu(324) = - lu(182) * lu(316)
         lu(328) = lu(328) - lu(183) * lu(316)
         lu(811) = lu(811) - lu(180) * lu(805)
         lu(813) = lu(813) - lu(181) * lu(805)
         lu(831) = lu(831) - lu(182) * lu(805)
         lu(842) = lu(842) - lu(183) * lu(805)
         lu(1039) = lu(1039) - lu(180) * lu(1022)
         lu(1041) = lu(1041) - lu(181) * lu(1022)
         lu(1066) = lu(1066) - lu(182) * lu(1022)
         lu(1078) = lu(1078) - lu(183) * lu(1022)
         lu(184) = 1._r8 / lu(184)
         lu(185) = lu(185) * lu(184)
         lu(186) = lu(186) * lu(184)
         lu(187) = lu(187) * lu(184)
         lu(188) = lu(188) * lu(184)
         lu(189) = lu(189) * lu(184)
         lu(190) = lu(190) * lu(184)
         lu(497) = lu(497) - lu(185) * lu(496)
         lu(499) = lu(499) - lu(186) * lu(496)
         lu(501) = lu(501) - lu(187) * lu(496)
         lu(503) = lu(503) - lu(188) * lu(496)
         lu(507) = lu(507) - lu(189) * lu(496)
         lu(510) = - lu(190) * lu(496)
         lu(870) = lu(870) - lu(185) * lu(855)
         lu(880) = lu(880) - lu(186) * lu(855)
         lu(882) = lu(882) - lu(187) * lu(855)
         lu(886) = lu(886) - lu(188) * lu(855)
         lu(897) = lu(897) - lu(189) * lu(855)
         lu(902) = lu(902) - lu(190) * lu(855)
         lu(1042) = lu(1042) - lu(185) * lu(1023)
         lu(1056) = lu(1056) - lu(186) * lu(1023)
         lu(1058) = lu(1058) - lu(187) * lu(1023)
         lu(1062) = lu(1062) - lu(188) * lu(1023)
         lu(1073) = lu(1073) - lu(189) * lu(1023)
         lu(1078) = lu(1078) - lu(190) * lu(1023)
         lu(191) = 1._r8 / lu(191)
         lu(192) = lu(192) * lu(191)
         lu(193) = lu(193) * lu(191)
         lu(194) = lu(194) * lu(191)
         lu(195) = lu(195) * lu(191)
         lu(196) = lu(196) * lu(191)
         lu(197) = lu(197) * lu(191)
         lu(646) = lu(646) - lu(192) * lu(633)
         lu(647) = lu(647) - lu(193) * lu(633)
         lu(648) = - lu(194) * lu(633)
         lu(652) = lu(652) - lu(195) * lu(633)
         lu(654) = lu(654) - lu(196) * lu(633)
         lu(655) = lu(655) - lu(197) * lu(633)
         lu(962) = lu(962) - lu(192) * lu(950)
         lu(963) = - lu(193) * lu(950)
         lu(965) = lu(965) - lu(194) * lu(950)
         lu(972) = lu(972) - lu(195) * lu(950)
         lu(974) = lu(974) - lu(196) * lu(950)
         lu(976) = - lu(197) * lu(950)
         lu(1066) = lu(1066) - lu(192) * lu(1024)
         lu(1067) = lu(1067) - lu(193) * lu(1024)
         lu(1069) = lu(1069) - lu(194) * lu(1024)
         lu(1076) = lu(1076) - lu(195) * lu(1024)
         lu(1078) = lu(1078) - lu(196) * lu(1024)
         lu(1080) = lu(1080) - lu(197) * lu(1024)
         lu(199) = 1._r8 / lu(199)
         lu(200) = lu(200) * lu(199)
         lu(201) = lu(201) * lu(199)
         lu(202) = lu(202) * lu(199)
         lu(203) = lu(203) * lu(199)
         lu(204) = lu(204) * lu(199)
         lu(205) = lu(205) * lu(199)
         lu(836) = lu(836) - lu(200) * lu(806)
         lu(837) = lu(837) - lu(201) * lu(806)
         lu(840) = lu(840) - lu(202) * lu(806)
         lu(841) = lu(841) - lu(203) * lu(806)
         lu(842) = lu(842) - lu(204) * lu(806)
         lu(844) = lu(844) - lu(205) * lu(806)
         lu(896) = lu(896) - lu(200) * lu(856)
         lu(897) = lu(897) - lu(201) * lu(856)
         lu(900) = lu(900) - lu(202) * lu(856)
         lu(901) = lu(901) - lu(203) * lu(856)
         lu(902) = lu(902) - lu(204) * lu(856)
         lu(904) = lu(904) - lu(205) * lu(856)
         lu(1109) = - lu(200) * lu(1100)
         lu(1110) = lu(1110) - lu(201) * lu(1100)
         lu(1113) = - lu(202) * lu(1100)
         lu(1114) = lu(1114) - lu(203) * lu(1100)
         lu(1115) = lu(1115) - lu(204) * lu(1100)
         lu(1117) = lu(1117) - lu(205) * lu(1100)
         lu(206) = 1._r8 / lu(206)
         lu(207) = lu(207) * lu(206)
         lu(208) = lu(208) * lu(206)
         lu(209) = lu(209) * lu(206)
         lu(210) = lu(210) * lu(206)
         lu(211) = lu(211) * lu(206)
         lu(622) = lu(622) - lu(207) * lu(616)
         lu(623) = lu(623) - lu(208) * lu(616)
         lu(629) = - lu(209) * lu(616)
         lu(630) = - lu(210) * lu(616)
         lu(631) = lu(631) - lu(211) * lu(616)
         lu(646) = lu(646) - lu(207) * lu(634)
         lu(647) = lu(647) - lu(208) * lu(634)
         lu(653) = - lu(209) * lu(634)
         lu(654) = lu(654) - lu(210) * lu(634)
         lu(655) = lu(655) - lu(211) * lu(634)
         lu(890) = lu(890) - lu(207) * lu(857)
         lu(891) = lu(891) - lu(208) * lu(857)
         lu(901) = lu(901) - lu(209) * lu(857)
         lu(902) = lu(902) - lu(210) * lu(857)
         lu(904) = lu(904) - lu(211) * lu(857)
         lu(1066) = lu(1066) - lu(207) * lu(1025)
         lu(1067) = lu(1067) - lu(208) * lu(1025)
         lu(1077) = lu(1077) - lu(209) * lu(1025)
         lu(1078) = lu(1078) - lu(210) * lu(1025)
         lu(1080) = lu(1080) - lu(211) * lu(1025)
      end subroutine lu_fac05
      subroutine lu_fac06( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(212) = 1._r8 / lu(212)
         lu(213) = lu(213) * lu(212)
         lu(214) = lu(214) * lu(212)
         lu(215) = lu(215) * lu(212)
         lu(216) = lu(216) * lu(212)
         lu(222) = lu(222) - lu(213) * lu(221)
         lu(224) = lu(224) - lu(214) * lu(221)
         lu(225) = - lu(215) * lu(221)
         lu(227) = lu(227) - lu(216) * lu(221)
         lu(808) = lu(808) - lu(213) * lu(807)
         lu(825) = lu(825) - lu(214) * lu(807)
         lu(831) = lu(831) - lu(215) * lu(807)
         lu(837) = lu(837) - lu(216) * lu(807)
         lu(859) = lu(859) - lu(213) * lu(858)
         lu(884) = lu(884) - lu(214) * lu(858)
         lu(890) = lu(890) - lu(215) * lu(858)
         lu(897) = lu(897) - lu(216) * lu(858)
         lu(952) = - lu(213) * lu(951)
         lu(960) = - lu(214) * lu(951)
         lu(962) = lu(962) - lu(215) * lu(951)
         lu(969) = lu(969) - lu(216) * lu(951)
         lu(1027) = lu(1027) - lu(213) * lu(1026)
         lu(1060) = lu(1060) - lu(214) * lu(1026)
         lu(1066) = lu(1066) - lu(215) * lu(1026)
         lu(1073) = lu(1073) - lu(216) * lu(1026)
         lu(217) = 1._r8 / lu(217)
         lu(218) = lu(218) * lu(217)
         lu(219) = lu(219) * lu(217)
         lu(227) = lu(227) - lu(218) * lu(222)
         lu(229) = lu(229) - lu(219) * lu(222)
         lu(273) = lu(273) - lu(218) * lu(271)
         lu(274) = lu(274) - lu(219) * lu(271)
         lu(556) = lu(556) - lu(218) * lu(546)
         lu(559) = lu(559) - lu(219) * lu(546)
         lu(651) = lu(651) - lu(218) * lu(635)
         lu(654) = lu(654) - lu(219) * lu(635)
         lu(682) = lu(682) - lu(218) * lu(659)
         lu(686) = lu(686) - lu(219) * lu(659)
         lu(736) = lu(736) - lu(218) * lu(713)
         lu(740) = lu(740) - lu(219) * lu(713)
         lu(837) = lu(837) - lu(218) * lu(808)
         lu(842) = lu(842) - lu(219) * lu(808)
         lu(897) = lu(897) - lu(218) * lu(859)
         lu(902) = lu(902) - lu(219) * lu(859)
         lu(969) = lu(969) - lu(218) * lu(952)
         lu(974) = lu(974) - lu(219) * lu(952)
         lu(1073) = lu(1073) - lu(218) * lu(1027)
         lu(1078) = lu(1078) - lu(219) * lu(1027)
         lu(223) = 1._r8 / lu(223)
         lu(224) = lu(224) * lu(223)
         lu(225) = lu(225) * lu(223)
         lu(226) = lu(226) * lu(223)
         lu(227) = lu(227) * lu(223)
         lu(228) = lu(228) * lu(223)
         lu(229) = lu(229) * lu(223)
         lu(825) = lu(825) - lu(224) * lu(809)
         lu(831) = lu(831) - lu(225) * lu(809)
         lu(836) = lu(836) - lu(226) * lu(809)
         lu(837) = lu(837) - lu(227) * lu(809)
         lu(840) = lu(840) - lu(228) * lu(809)
         lu(842) = lu(842) - lu(229) * lu(809)
         lu(884) = lu(884) - lu(224) * lu(860)
         lu(890) = lu(890) - lu(225) * lu(860)
         lu(896) = lu(896) - lu(226) * lu(860)
         lu(897) = lu(897) - lu(227) * lu(860)
         lu(900) = lu(900) - lu(228) * lu(860)
         lu(902) = lu(902) - lu(229) * lu(860)
         lu(1060) = lu(1060) - lu(224) * lu(1028)
         lu(1066) = lu(1066) - lu(225) * lu(1028)
         lu(1072) = lu(1072) - lu(226) * lu(1028)
         lu(1073) = lu(1073) - lu(227) * lu(1028)
         lu(1076) = lu(1076) - lu(228) * lu(1028)
         lu(1078) = lu(1078) - lu(229) * lu(1028)
         lu(230) = 1._r8 / lu(230)
         lu(231) = lu(231) * lu(230)
         lu(232) = lu(232) * lu(230)
         lu(233) = lu(233) * lu(230)
         lu(234) = lu(234) * lu(230)
         lu(529) = lu(529) - lu(231) * lu(528)
         lu(532) = - lu(232) * lu(528)
         lu(534) = - lu(233) * lu(528)
         lu(535) = - lu(234) * lu(528)
         lu(883) = lu(883) - lu(231) * lu(861)
         lu(898) = lu(898) - lu(232) * lu(861)
         lu(901) = lu(901) - lu(233) * lu(861)
         lu(902) = lu(902) - lu(234) * lu(861)
         lu(913) = lu(913) - lu(231) * lu(905)
         lu(920) = lu(920) - lu(232) * lu(905)
         lu(923) = - lu(233) * lu(905)
         lu(924) = lu(924) - lu(234) * lu(905)
         lu(1059) = lu(1059) - lu(231) * lu(1029)
         lu(1074) = lu(1074) - lu(232) * lu(1029)
         lu(1077) = lu(1077) - lu(233) * lu(1029)
         lu(1078) = lu(1078) - lu(234) * lu(1029)
         lu(1104) = lu(1104) - lu(231) * lu(1101)
         lu(1111) = lu(1111) - lu(232) * lu(1101)
         lu(1114) = lu(1114) - lu(233) * lu(1101)
         lu(1115) = lu(1115) - lu(234) * lu(1101)
         lu(235) = 1._r8 / lu(235)
         lu(236) = lu(236) * lu(235)
         lu(237) = lu(237) * lu(235)
         lu(238) = lu(238) * lu(235)
         lu(239) = lu(239) * lu(235)
         lu(240) = lu(240) * lu(235)
         lu(241) = lu(241) * lu(235)
         lu(242) = lu(242) * lu(235)
         lu(360) = lu(360) - lu(236) * lu(359)
         lu(361) = lu(361) - lu(237) * lu(359)
         lu(362) = - lu(238) * lu(359)
         lu(364) = lu(364) - lu(239) * lu(359)
         lu(366) = - lu(240) * lu(359)
         lu(367) = - lu(241) * lu(359)
         lu(368) = lu(368) - lu(242) * lu(359)
         lu(871) = lu(871) - lu(236) * lu(862)
         lu(877) = lu(877) - lu(237) * lu(862)
         lu(881) = - lu(238) * lu(862)
         lu(897) = lu(897) - lu(239) * lu(862)
         lu(901) = lu(901) - lu(240) * lu(862)
         lu(902) = lu(902) - lu(241) * lu(862)
         lu(904) = lu(904) - lu(242) * lu(862)
         lu(1044) = lu(1044) - lu(236) * lu(1030)
         lu(1053) = lu(1053) - lu(237) * lu(1030)
         lu(1057) = lu(1057) - lu(238) * lu(1030)
         lu(1073) = lu(1073) - lu(239) * lu(1030)
         lu(1077) = lu(1077) - lu(240) * lu(1030)
         lu(1078) = lu(1078) - lu(241) * lu(1030)
         lu(1080) = lu(1080) - lu(242) * lu(1030)
         lu(243) = 1._r8 / lu(243)
         lu(244) = lu(244) * lu(243)
         lu(245) = lu(245) * lu(243)
         lu(246) = lu(246) * lu(243)
         lu(247) = lu(247) * lu(243)
         lu(248) = lu(248) * lu(243)
         lu(249) = lu(249) * lu(243)
         lu(250) = lu(250) * lu(243)
         lu(599) = lu(599) - lu(244) * lu(596)
         lu(601) = lu(601) - lu(245) * lu(596)
         lu(602) = lu(602) - lu(246) * lu(596)
         lu(604) = lu(604) - lu(247) * lu(596)
         lu(611) = lu(611) - lu(248) * lu(596)
         lu(614) = lu(614) - lu(249) * lu(596)
         lu(615) = lu(615) - lu(250) * lu(596)
         lu(882) = lu(882) - lu(244) * lu(863)
         lu(885) = lu(885) - lu(245) * lu(863)
         lu(886) = lu(886) - lu(246) * lu(863)
         lu(888) = lu(888) - lu(247) * lu(863)
         lu(897) = lu(897) - lu(248) * lu(863)
         lu(902) = lu(902) - lu(249) * lu(863)
         lu(904) = lu(904) - lu(250) * lu(863)
         lu(1058) = lu(1058) - lu(244) * lu(1031)
         lu(1061) = lu(1061) - lu(245) * lu(1031)
         lu(1062) = lu(1062) - lu(246) * lu(1031)
         lu(1064) = lu(1064) - lu(247) * lu(1031)
         lu(1073) = lu(1073) - lu(248) * lu(1031)
         lu(1078) = lu(1078) - lu(249) * lu(1031)
         lu(1080) = lu(1080) - lu(250) * lu(1031)
         lu(251) = 1._r8 / lu(251)
         lu(252) = lu(252) * lu(251)
         lu(253) = lu(253) * lu(251)
         lu(254) = lu(254) * lu(251)
         lu(255) = lu(255) * lu(251)
         lu(256) = lu(256) * lu(251)
         lu(257) = lu(257) * lu(251)
         lu(258) = lu(258) * lu(251)
         lu(908) = lu(908) - lu(252) * lu(906)
         lu(909) = - lu(253) * lu(906)
         lu(913) = lu(913) - lu(254) * lu(906)
         lu(915) = lu(915) - lu(255) * lu(906)
         lu(920) = lu(920) - lu(256) * lu(906)
         lu(921) = lu(921) - lu(257) * lu(906)
         lu(922) = lu(922) - lu(258) * lu(906)
         lu(930) = lu(930) - lu(252) * lu(929)
         lu(931) = - lu(253) * lu(929)
         lu(932) = lu(932) - lu(254) * lu(929)
         lu(934) = - lu(255) * lu(929)
         lu(939) = lu(939) - lu(256) * lu(929)
         lu(940) = lu(940) - lu(257) * lu(929)
         lu(941) = lu(941) - lu(258) * lu(929)
         lu(955) = - lu(252) * lu(953)
         lu(956) = lu(956) - lu(253) * lu(953)
         lu(959) = - lu(254) * lu(953)
         lu(965) = lu(965) - lu(255) * lu(953)
         lu(970) = lu(970) - lu(256) * lu(953)
         lu(971) = lu(971) - lu(257) * lu(953)
         lu(972) = lu(972) - lu(258) * lu(953)
      end subroutine lu_fac06
      subroutine lu_fac07( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(259) = 1._r8 / lu(259)
         lu(260) = lu(260) * lu(259)
         lu(261) = lu(261) * lu(259)
         lu(262) = lu(262) * lu(259)
         lu(263) = lu(263) * lu(259)
         lu(264) = lu(264) * lu(259)
         lu(265) = lu(265) * lu(259)
         lu(266) = lu(266) * lu(259)
         lu(619) = - lu(260) * lu(617)
         lu(621) = lu(621) - lu(261) * lu(617)
         lu(624) = lu(624) - lu(262) * lu(617)
         lu(627) = lu(627) - lu(263) * lu(617)
         lu(628) = lu(628) - lu(264) * lu(617)
         lu(630) = lu(630) - lu(265) * lu(617)
         lu(631) = lu(631) - lu(266) * lu(617)
         lu(958) = - lu(260) * lu(954)
         lu(961) = lu(961) - lu(261) * lu(954)
         lu(965) = lu(965) - lu(262) * lu(954)
         lu(969) = lu(969) - lu(263) * lu(954)
         lu(972) = lu(972) - lu(264) * lu(954)
         lu(974) = lu(974) - lu(265) * lu(954)
         lu(976) = lu(976) - lu(266) * lu(954)
         lu(1057) = lu(1057) - lu(260) * lu(1032)
         lu(1065) = lu(1065) - lu(261) * lu(1032)
         lu(1069) = lu(1069) - lu(262) * lu(1032)
         lu(1073) = lu(1073) - lu(263) * lu(1032)
         lu(1076) = lu(1076) - lu(264) * lu(1032)
         lu(1078) = lu(1078) - lu(265) * lu(1032)
         lu(1080) = lu(1080) - lu(266) * lu(1032)
         lu(267) = 1._r8 / lu(267)
         lu(268) = lu(268) * lu(267)
         lu(269) = lu(269) * lu(267)
         lu(270) = lu(270) * lu(267)
         lu(337) = lu(337) - lu(268) * lu(332)
         lu(339) = lu(339) - lu(269) * lu(332)
         lu(340) = lu(340) - lu(270) * lu(332)
         lu(478) = lu(478) - lu(268) * lu(471)
         lu(481) = lu(481) - lu(269) * lu(471)
         lu(482) = lu(482) - lu(270) * lu(471)
         lu(556) = lu(556) - lu(268) * lu(547)
         lu(559) = lu(559) - lu(269) * lu(547)
         lu(560) = lu(560) - lu(270) * lu(547)
         lu(590) = lu(590) - lu(268) * lu(578)
         lu(593) = lu(593) - lu(269) * lu(578)
         lu(594) = lu(594) - lu(270) * lu(578)
         lu(611) = lu(611) - lu(268) * lu(597)
         lu(614) = lu(614) - lu(269) * lu(597)
         lu(615) = lu(615) - lu(270) * lu(597)
         lu(682) = lu(682) - lu(268) * lu(660)
         lu(686) = lu(686) - lu(269) * lu(660)
         lu(687) = lu(687) - lu(270) * lu(660)
         lu(1073) = lu(1073) - lu(268) * lu(1033)
         lu(1078) = lu(1078) - lu(269) * lu(1033)
         lu(1080) = lu(1080) - lu(270) * lu(1033)
         lu(272) = 1._r8 / lu(272)
         lu(273) = lu(273) * lu(272)
         lu(274) = lu(274) * lu(272)
         lu(275) = lu(275) * lu(272)
         lu(556) = lu(556) - lu(273) * lu(548)
         lu(559) = lu(559) - lu(274) * lu(548)
         lu(560) = lu(560) - lu(275) * lu(548)
         lu(590) = lu(590) - lu(273) * lu(579)
         lu(593) = lu(593) - lu(274) * lu(579)
         lu(594) = lu(594) - lu(275) * lu(579)
         lu(651) = lu(651) - lu(273) * lu(636)
         lu(654) = lu(654) - lu(274) * lu(636)
         lu(655) = lu(655) - lu(275) * lu(636)
         lu(682) = lu(682) - lu(273) * lu(661)
         lu(686) = lu(686) - lu(274) * lu(661)
         lu(687) = lu(687) - lu(275) * lu(661)
         lu(736) = lu(736) - lu(273) * lu(714)
         lu(740) = lu(740) - lu(274) * lu(714)
         lu(741) = lu(741) - lu(275) * lu(714)
         lu(790) = lu(790) - lu(273) * lu(765)
         lu(795) = lu(795) - lu(274) * lu(765)
         lu(797) = lu(797) - lu(275) * lu(765)
         lu(837) = lu(837) - lu(273) * lu(810)
         lu(842) = lu(842) - lu(274) * lu(810)
         lu(844) = lu(844) - lu(275) * lu(810)
         lu(1073) = lu(1073) - lu(273) * lu(1034)
         lu(1078) = lu(1078) - lu(274) * lu(1034)
         lu(1080) = lu(1080) - lu(275) * lu(1034)
         lu(276) = 1._r8 / lu(276)
         lu(277) = lu(277) * lu(276)
         lu(278) = lu(278) * lu(276)
         lu(279) = lu(279) * lu(276)
         lu(386) = - lu(277) * lu(378)
         lu(390) = - lu(278) * lu(378)
         lu(391) = lu(391) - lu(279) * lu(378)
         lu(452) = lu(452) - lu(277) * lu(446)
         lu(458) = - lu(278) * lu(446)
         lu(459) = lu(459) - lu(279) * lu(446)
         lu(623) = lu(623) - lu(277) * lu(618)
         lu(629) = lu(629) - lu(278) * lu(618)
         lu(630) = lu(630) - lu(279) * lu(618)
         lu(647) = lu(647) - lu(277) * lu(637)
         lu(653) = lu(653) - lu(278) * lu(637)
         lu(654) = lu(654) - lu(279) * lu(637)
         lu(678) = lu(678) - lu(277) * lu(662)
         lu(685) = lu(685) - lu(278) * lu(662)
         lu(686) = lu(686) - lu(279) * lu(662)
         lu(785) = lu(785) - lu(277) * lu(766)
         lu(794) = lu(794) - lu(278) * lu(766)
         lu(795) = lu(795) - lu(279) * lu(766)
         lu(891) = lu(891) - lu(277) * lu(864)
         lu(901) = lu(901) - lu(278) * lu(864)
         lu(902) = lu(902) - lu(279) * lu(864)
         lu(1067) = lu(1067) - lu(277) * lu(1035)
         lu(1077) = lu(1077) - lu(278) * lu(1035)
         lu(1078) = lu(1078) - lu(279) * lu(1035)
         lu(280) = 1._r8 / lu(280)
         lu(281) = lu(281) * lu(280)
         lu(282) = lu(282) * lu(280)
         lu(283) = lu(283) * lu(280)
         lu(284) = lu(284) * lu(280)
         lu(285) = lu(285) * lu(280)
         lu(286) = lu(286) * lu(280)
         lu(287) = lu(287) * lu(280)
         lu(288) = lu(288) * lu(280)
         lu(722) = lu(722) - lu(281) * lu(715)
         lu(724) = lu(724) - lu(282) * lu(715)
         lu(727) = lu(727) - lu(283) * lu(715)
         lu(733) = lu(733) - lu(284) * lu(715)
         lu(734) = - lu(285) * lu(715)
         lu(736) = lu(736) - lu(286) * lu(715)
         lu(738) = lu(738) - lu(287) * lu(715)
         lu(740) = lu(740) - lu(288) * lu(715)
         lu(774) = - lu(281) * lu(767)
         lu(776) = lu(776) - lu(282) * lu(767)
         lu(780) = lu(780) - lu(283) * lu(767)
         lu(786) = lu(786) - lu(284) * lu(767)
         lu(788) = lu(788) - lu(285) * lu(767)
         lu(790) = lu(790) - lu(286) * lu(767)
         lu(793) = lu(793) - lu(287) * lu(767)
         lu(795) = lu(795) - lu(288) * lu(767)
         lu(1056) = lu(1056) - lu(281) * lu(1036)
         lu(1058) = lu(1058) - lu(282) * lu(1036)
         lu(1062) = lu(1062) - lu(283) * lu(1036)
         lu(1069) = lu(1069) - lu(284) * lu(1036)
         lu(1071) = lu(1071) - lu(285) * lu(1036)
         lu(1073) = lu(1073) - lu(286) * lu(1036)
         lu(1076) = lu(1076) - lu(287) * lu(1036)
         lu(1078) = lu(1078) - lu(288) * lu(1036)
         lu(289) = 1._r8 / lu(289)
         lu(290) = lu(290) * lu(289)
         lu(291) = lu(291) * lu(289)
         lu(292) = lu(292) * lu(289)
         lu(293) = lu(293) * lu(289)
         lu(294) = lu(294) * lu(289)
         lu(295) = lu(295) * lu(289)
         lu(892) = lu(892) - lu(290) * lu(865)
         lu(897) = lu(897) - lu(291) * lu(865)
         lu(898) = lu(898) - lu(292) * lu(865)
         lu(901) = lu(901) - lu(293) * lu(865)
         lu(902) = lu(902) - lu(294) * lu(865)
         lu(903) = lu(903) - lu(295) * lu(865)
         lu(914) = lu(914) - lu(290) * lu(907)
         lu(919) = lu(919) - lu(291) * lu(907)
         lu(920) = lu(920) - lu(292) * lu(907)
         lu(923) = lu(923) - lu(293) * lu(907)
         lu(924) = lu(924) - lu(294) * lu(907)
         lu(925) = lu(925) - lu(295) * lu(907)
         lu(1068) = lu(1068) - lu(290) * lu(1037)
         lu(1073) = lu(1073) - lu(291) * lu(1037)
         lu(1074) = lu(1074) - lu(292) * lu(1037)
         lu(1077) = lu(1077) - lu(293) * lu(1037)
         lu(1078) = lu(1078) - lu(294) * lu(1037)
         lu(1079) = lu(1079) - lu(295) * lu(1037)
         lu(1087) = lu(1087) - lu(290) * lu(1082)
         lu(1092) = lu(1092) - lu(291) * lu(1082)
         lu(1093) = - lu(292) * lu(1082)
         lu(1096) = - lu(293) * lu(1082)
         lu(1097) = lu(1097) - lu(294) * lu(1082)
         lu(1098) = lu(1098) - lu(295) * lu(1082)
         lu(297) = 1._r8 / lu(297)
         lu(298) = lu(298) * lu(297)
         lu(299) = lu(299) * lu(297)
         lu(300) = lu(300) * lu(297)
         lu(301) = lu(301) * lu(297)
         lu(302) = lu(302) * lu(297)
         lu(303) = lu(303) * lu(297)
         lu(304) = lu(304) * lu(297)
         lu(305) = lu(305) * lu(297)
         lu(306) = lu(306) * lu(297)
         lu(318) = lu(318) - lu(298) * lu(317)
         lu(319) = lu(319) - lu(299) * lu(317)
         lu(320) = lu(320) - lu(300) * lu(317)
         lu(321) = lu(321) - lu(301) * lu(317)
         lu(322) = lu(322) - lu(302) * lu(317)
         lu(324) = lu(324) - lu(303) * lu(317)
         lu(326) = lu(326) - lu(304) * lu(317)
         lu(328) = lu(328) - lu(305) * lu(317)
         lu(329) = lu(329) - lu(306) * lu(317)
         lu(867) = lu(867) - lu(298) * lu(866)
         lu(868) = lu(868) - lu(299) * lu(866)
         lu(869) = lu(869) - lu(300) * lu(866)
         lu(870) = lu(870) - lu(301) * lu(866)
         lu(877) = lu(877) - lu(302) * lu(866)
         lu(890) = lu(890) - lu(303) * lu(866)
         lu(897) = lu(897) - lu(304) * lu(866)
         lu(902) = lu(902) - lu(305) * lu(866)
         lu(904) = lu(904) - lu(306) * lu(866)
         lu(1039) = lu(1039) - lu(298) * lu(1038)
         lu(1040) = lu(1040) - lu(299) * lu(1038)
         lu(1041) = lu(1041) - lu(300) * lu(1038)
         lu(1042) = lu(1042) - lu(301) * lu(1038)
         lu(1053) = lu(1053) - lu(302) * lu(1038)
         lu(1066) = lu(1066) - lu(303) * lu(1038)
         lu(1073) = lu(1073) - lu(304) * lu(1038)
         lu(1078) = lu(1078) - lu(305) * lu(1038)
         lu(1080) = lu(1080) - lu(306) * lu(1038)
      end subroutine lu_fac07
      subroutine lu_fac08( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(308) = 1._r8 / lu(308)
         lu(309) = lu(309) * lu(308)
         lu(310) = lu(310) * lu(308)
         lu(311) = lu(311) * lu(308)
         lu(312) = lu(312) * lu(308)
         lu(313) = lu(313) * lu(308)
         lu(314) = lu(314) * lu(308)
         lu(322) = lu(322) - lu(309) * lu(318)
         lu(324) = lu(324) - lu(310) * lu(318)
         lu(325) = lu(325) - lu(311) * lu(318)
         lu(326) = lu(326) - lu(312) * lu(318)
         lu(327) = lu(327) - lu(313) * lu(318)
         lu(328) = lu(328) - lu(314) * lu(318)
         lu(818) = lu(818) - lu(309) * lu(811)
         lu(831) = lu(831) - lu(310) * lu(811)
         lu(836) = lu(836) - lu(311) * lu(811)
         lu(837) = lu(837) - lu(312) * lu(811)
         lu(840) = lu(840) - lu(313) * lu(811)
         lu(842) = lu(842) - lu(314) * lu(811)
         lu(877) = lu(877) - lu(309) * lu(867)
         lu(890) = lu(890) - lu(310) * lu(867)
         lu(896) = lu(896) - lu(311) * lu(867)
         lu(897) = lu(897) - lu(312) * lu(867)
         lu(900) = lu(900) - lu(313) * lu(867)
         lu(902) = lu(902) - lu(314) * lu(867)
         lu(1053) = lu(1053) - lu(309) * lu(1039)
         lu(1066) = lu(1066) - lu(310) * lu(1039)
         lu(1072) = lu(1072) - lu(311) * lu(1039)
         lu(1073) = lu(1073) - lu(312) * lu(1039)
         lu(1076) = lu(1076) - lu(313) * lu(1039)
         lu(1078) = lu(1078) - lu(314) * lu(1039)
         lu(319) = 1._r8 / lu(319)
         lu(320) = lu(320) * lu(319)
         lu(321) = lu(321) * lu(319)
         lu(322) = lu(322) * lu(319)
         lu(323) = lu(323) * lu(319)
         lu(324) = lu(324) * lu(319)
         lu(325) = lu(325) * lu(319)
         lu(326) = lu(326) * lu(319)
         lu(327) = lu(327) * lu(319)
         lu(328) = lu(328) * lu(319)
         lu(329) = lu(329) * lu(319)
         lu(813) = lu(813) - lu(320) * lu(812)
         lu(814) = lu(814) - lu(321) * lu(812)
         lu(818) = lu(818) - lu(322) * lu(812)
         lu(825) = lu(825) - lu(323) * lu(812)
         lu(831) = lu(831) - lu(324) * lu(812)
         lu(836) = lu(836) - lu(325) * lu(812)
         lu(837) = lu(837) - lu(326) * lu(812)
         lu(840) = lu(840) - lu(327) * lu(812)
         lu(842) = lu(842) - lu(328) * lu(812)
         lu(844) = lu(844) - lu(329) * lu(812)
         lu(869) = lu(869) - lu(320) * lu(868)
         lu(870) = lu(870) - lu(321) * lu(868)
         lu(877) = lu(877) - lu(322) * lu(868)
         lu(884) = lu(884) - lu(323) * lu(868)
         lu(890) = lu(890) - lu(324) * lu(868)
         lu(896) = lu(896) - lu(325) * lu(868)
         lu(897) = lu(897) - lu(326) * lu(868)
         lu(900) = lu(900) - lu(327) * lu(868)
         lu(902) = lu(902) - lu(328) * lu(868)
         lu(904) = lu(904) - lu(329) * lu(868)
         lu(1041) = lu(1041) - lu(320) * lu(1040)
         lu(1042) = lu(1042) - lu(321) * lu(1040)
         lu(1053) = lu(1053) - lu(322) * lu(1040)
         lu(1060) = lu(1060) - lu(323) * lu(1040)
         lu(1066) = lu(1066) - lu(324) * lu(1040)
         lu(1072) = lu(1072) - lu(325) * lu(1040)
         lu(1073) = lu(1073) - lu(326) * lu(1040)
         lu(1076) = lu(1076) - lu(327) * lu(1040)
         lu(1078) = lu(1078) - lu(328) * lu(1040)
         lu(1080) = lu(1080) - lu(329) * lu(1040)
         lu(333) = 1._r8 / lu(333)
         lu(334) = lu(334) * lu(333)
         lu(335) = lu(335) * lu(333)
         lu(336) = lu(336) * lu(333)
         lu(337) = lu(337) * lu(333)
         lu(338) = lu(338) * lu(333)
         lu(339) = lu(339) * lu(333)
         lu(340) = lu(340) * lu(333)
         lu(667) = lu(667) - lu(334) * lu(663)
         lu(678) = lu(678) - lu(335) * lu(663)
         lu(681) = lu(681) - lu(336) * lu(663)
         lu(682) = lu(682) - lu(337) * lu(663)
         lu(684) = lu(684) - lu(338) * lu(663)
         lu(686) = lu(686) - lu(339) * lu(663)
         lu(687) = lu(687) - lu(340) * lu(663)
         lu(818) = lu(818) - lu(334) * lu(813)
         lu(832) = lu(832) - lu(335) * lu(813)
         lu(836) = lu(836) - lu(336) * lu(813)
         lu(837) = lu(837) - lu(337) * lu(813)
         lu(840) = lu(840) - lu(338) * lu(813)
         lu(842) = lu(842) - lu(339) * lu(813)
         lu(844) = lu(844) - lu(340) * lu(813)
         lu(877) = lu(877) - lu(334) * lu(869)
         lu(891) = lu(891) - lu(335) * lu(869)
         lu(896) = lu(896) - lu(336) * lu(869)
         lu(897) = lu(897) - lu(337) * lu(869)
         lu(900) = lu(900) - lu(338) * lu(869)
         lu(902) = lu(902) - lu(339) * lu(869)
         lu(904) = lu(904) - lu(340) * lu(869)
         lu(1053) = lu(1053) - lu(334) * lu(1041)
         lu(1067) = lu(1067) - lu(335) * lu(1041)
         lu(1072) = lu(1072) - lu(336) * lu(1041)
         lu(1073) = lu(1073) - lu(337) * lu(1041)
         lu(1076) = lu(1076) - lu(338) * lu(1041)
         lu(1078) = lu(1078) - lu(339) * lu(1041)
         lu(1080) = lu(1080) - lu(340) * lu(1041)
         lu(341) = 1._r8 / lu(341)
         lu(342) = lu(342) * lu(341)
         lu(343) = lu(343) * lu(341)
         lu(344) = lu(344) * lu(341)
         lu(345) = lu(345) * lu(341)
         lu(346) = lu(346) * lu(341)
         lu(397) = - lu(342) * lu(394)
         lu(398) = - lu(343) * lu(394)
         lu(399) = lu(399) - lu(344) * lu(394)
         lu(403) = lu(403) - lu(345) * lu(394)
         lu(404) = lu(404) - lu(346) * lu(394)
         lu(498) = - lu(342) * lu(497)
         lu(504) = - lu(343) * lu(497)
         lu(505) = - lu(344) * lu(497)
         lu(509) = - lu(345) * lu(497)
         lu(510) = lu(510) - lu(346) * lu(497)
         lu(668) = lu(668) - lu(342) * lu(664)
         lu(677) = lu(677) - lu(343) * lu(664)
         lu(678) = lu(678) - lu(344) * lu(664)
         lu(685) = lu(685) - lu(345) * lu(664)
         lu(686) = lu(686) - lu(346) * lu(664)
         lu(819) = lu(819) - lu(342) * lu(814)
         lu(831) = lu(831) - lu(343) * lu(814)
         lu(832) = lu(832) - lu(344) * lu(814)
         lu(841) = lu(841) - lu(345) * lu(814)
         lu(842) = lu(842) - lu(346) * lu(814)
         lu(878) = lu(878) - lu(342) * lu(870)
         lu(890) = lu(890) - lu(343) * lu(870)
         lu(891) = lu(891) - lu(344) * lu(870)
         lu(901) = lu(901) - lu(345) * lu(870)
         lu(902) = lu(902) - lu(346) * lu(870)
         lu(1054) = lu(1054) - lu(342) * lu(1042)
         lu(1066) = lu(1066) - lu(343) * lu(1042)
         lu(1067) = lu(1067) - lu(344) * lu(1042)
         lu(1077) = lu(1077) - lu(345) * lu(1042)
         lu(1078) = lu(1078) - lu(346) * lu(1042)
         lu(348) = 1._r8 / lu(348)
         lu(349) = lu(349) * lu(348)
         lu(350) = lu(350) * lu(348)
         lu(351) = lu(351) * lu(348)
         lu(352) = lu(352) * lu(348)
         lu(353) = lu(353) * lu(348)
         lu(354) = lu(354) * lu(348)
         lu(355) = lu(355) * lu(348)
         lu(356) = lu(356) * lu(348)
         lu(357) = lu(357) * lu(348)
         lu(358) = lu(358) * lu(348)
         lu(770) = lu(770) - lu(349) * lu(768)
         lu(777) = lu(777) - lu(350) * lu(768)
         lu(785) = lu(785) - lu(351) * lu(768)
         lu(789) = lu(789) - lu(352) * lu(768)
         lu(790) = lu(790) - lu(353) * lu(768)
         lu(791) = lu(791) - lu(354) * lu(768)
         lu(794) = lu(794) - lu(355) * lu(768)
         lu(795) = lu(795) - lu(356) * lu(768)
         lu(796) = lu(796) - lu(357) * lu(768)
         lu(797) = lu(797) - lu(358) * lu(768)
         lu(978) = lu(978) - lu(349) * lu(977)
         lu(979) = - lu(350) * lu(977)
         lu(980) = - lu(351) * lu(977)
         lu(983) = - lu(352) * lu(977)
         lu(984) = - lu(353) * lu(977)
         lu(985) = lu(985) - lu(354) * lu(977)
         lu(988) = lu(988) - lu(355) * lu(977)
         lu(989) = lu(989) - lu(356) * lu(977)
         lu(990) = - lu(357) * lu(977)
         lu(991) = - lu(358) * lu(977)
         lu(1050) = lu(1050) - lu(349) * lu(1043)
         lu(1059) = lu(1059) - lu(350) * lu(1043)
         lu(1067) = lu(1067) - lu(351) * lu(1043)
         lu(1072) = lu(1072) - lu(352) * lu(1043)
         lu(1073) = lu(1073) - lu(353) * lu(1043)
         lu(1074) = lu(1074) - lu(354) * lu(1043)
         lu(1077) = lu(1077) - lu(355) * lu(1043)
         lu(1078) = lu(1078) - lu(356) * lu(1043)
         lu(1079) = lu(1079) - lu(357) * lu(1043)
         lu(1080) = lu(1080) - lu(358) * lu(1043)
      end subroutine lu_fac08
      subroutine lu_fac09( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(360) = 1._r8 / lu(360)
         lu(361) = lu(361) * lu(360)
         lu(362) = lu(362) * lu(360)
         lu(363) = lu(363) * lu(360)
         lu(364) = lu(364) * lu(360)
         lu(365) = lu(365) * lu(360)
         lu(366) = lu(366) * lu(360)
         lu(367) = lu(367) * lu(360)
         lu(368) = lu(368) * lu(360)
         lu(449) = lu(449) - lu(361) * lu(447)
         lu(450) = - lu(362) * lu(447)
         lu(455) = - lu(363) * lu(447)
         lu(456) = lu(456) - lu(364) * lu(447)
         lu(457) = lu(457) - lu(365) * lu(447)
         lu(458) = lu(458) - lu(366) * lu(447)
         lu(459) = lu(459) - lu(367) * lu(447)
         lu(460) = lu(460) - lu(368) * lu(447)
         lu(818) = lu(818) - lu(361) * lu(815)
         lu(822) = lu(822) - lu(362) * lu(815)
         lu(836) = lu(836) - lu(363) * lu(815)
         lu(837) = lu(837) - lu(364) * lu(815)
         lu(840) = lu(840) - lu(365) * lu(815)
         lu(841) = lu(841) - lu(366) * lu(815)
         lu(842) = lu(842) - lu(367) * lu(815)
         lu(844) = lu(844) - lu(368) * lu(815)
         lu(877) = lu(877) - lu(361) * lu(871)
         lu(881) = lu(881) - lu(362) * lu(871)
         lu(896) = lu(896) - lu(363) * lu(871)
         lu(897) = lu(897) - lu(364) * lu(871)
         lu(900) = lu(900) - lu(365) * lu(871)
         lu(901) = lu(901) - lu(366) * lu(871)
         lu(902) = lu(902) - lu(367) * lu(871)
         lu(904) = lu(904) - lu(368) * lu(871)
         lu(1053) = lu(1053) - lu(361) * lu(1044)
         lu(1057) = lu(1057) - lu(362) * lu(1044)
         lu(1072) = lu(1072) - lu(363) * lu(1044)
         lu(1073) = lu(1073) - lu(364) * lu(1044)
         lu(1076) = lu(1076) - lu(365) * lu(1044)
         lu(1077) = lu(1077) - lu(366) * lu(1044)
         lu(1078) = lu(1078) - lu(367) * lu(1044)
         lu(1080) = lu(1080) - lu(368) * lu(1044)
         lu(370) = 1._r8 / lu(370)
         lu(371) = lu(371) * lu(370)
         lu(372) = lu(372) * lu(370)
         lu(373) = lu(373) * lu(370)
         lu(374) = lu(374) * lu(370)
         lu(375) = lu(375) * lu(370)
         lu(376) = lu(376) * lu(370)
         lu(377) = lu(377) * lu(370)
         lu(695) = lu(695) - lu(371) * lu(690)
         lu(696) = lu(696) - lu(372) * lu(690)
         lu(701) = lu(701) - lu(373) * lu(690)
         lu(702) = - lu(374) * lu(690)
         lu(704) = lu(704) - lu(375) * lu(690)
         lu(705) = lu(705) - lu(376) * lu(690)
         lu(706) = lu(706) - lu(377) * lu(690)
         lu(883) = lu(883) - lu(371) * lu(872)
         lu(892) = lu(892) - lu(372) * lu(872)
         lu(898) = lu(898) - lu(373) * lu(872)
         lu(899) = lu(899) - lu(374) * lu(872)
         lu(901) = lu(901) - lu(375) * lu(872)
         lu(902) = lu(902) - lu(376) * lu(872)
         lu(903) = lu(903) - lu(377) * lu(872)
         lu(913) = lu(913) - lu(371) * lu(908)
         lu(914) = lu(914) - lu(372) * lu(908)
         lu(920) = lu(920) - lu(373) * lu(908)
         lu(921) = lu(921) - lu(374) * lu(908)
         lu(923) = lu(923) - lu(375) * lu(908)
         lu(924) = lu(924) - lu(376) * lu(908)
         lu(925) = lu(925) - lu(377) * lu(908)
         lu(932) = lu(932) - lu(371) * lu(930)
         lu(933) = - lu(372) * lu(930)
         lu(939) = lu(939) - lu(373) * lu(930)
         lu(940) = lu(940) - lu(374) * lu(930)
         lu(942) = - lu(375) * lu(930)
         lu(943) = lu(943) - lu(376) * lu(930)
         lu(944) = lu(944) - lu(377) * lu(930)
         lu(959) = lu(959) - lu(371) * lu(955)
         lu(964) = - lu(372) * lu(955)
         lu(970) = lu(970) - lu(373) * lu(955)
         lu(971) = lu(971) - lu(374) * lu(955)
         lu(973) = lu(973) - lu(375) * lu(955)
         lu(974) = lu(974) - lu(376) * lu(955)
         lu(975) = - lu(377) * lu(955)
         lu(379) = 1._r8 / lu(379)
         lu(380) = lu(380) * lu(379)
         lu(381) = lu(381) * lu(379)
         lu(382) = lu(382) * lu(379)
         lu(383) = lu(383) * lu(379)
         lu(384) = lu(384) * lu(379)
         lu(385) = lu(385) * lu(379)
         lu(386) = lu(386) * lu(379)
         lu(387) = lu(387) * lu(379)
         lu(388) = lu(388) * lu(379)
         lu(389) = lu(389) * lu(379)
         lu(390) = lu(390) * lu(379)
         lu(391) = lu(391) * lu(379)
         lu(392) = lu(392) * lu(379)
         lu(719) = lu(719) - lu(380) * lu(716)
         lu(721) = lu(721) - lu(381) * lu(716)
         lu(724) = lu(724) - lu(382) * lu(716)
         lu(727) = lu(727) - lu(383) * lu(716)
         lu(729) = lu(729) - lu(384) * lu(716)
         lu(730) = lu(730) - lu(385) * lu(716)
         lu(732) = - lu(386) * lu(716)
         lu(733) = lu(733) - lu(387) * lu(716)
         lu(734) = lu(734) - lu(388) * lu(716)
         lu(736) = lu(736) - lu(389) * lu(716)
         lu(739) = - lu(390) * lu(716)
         lu(740) = lu(740) - lu(391) * lu(716)
         lu(741) = lu(741) - lu(392) * lu(716)
         lu(771) = lu(771) - lu(380) * lu(769)
         lu(773) = - lu(381) * lu(769)
         lu(776) = lu(776) - lu(382) * lu(769)
         lu(780) = lu(780) - lu(383) * lu(769)
         lu(782) = - lu(384) * lu(769)
         lu(783) = lu(783) - lu(385) * lu(769)
         lu(785) = lu(785) - lu(386) * lu(769)
         lu(786) = lu(786) - lu(387) * lu(769)
         lu(788) = lu(788) - lu(388) * lu(769)
         lu(790) = lu(790) - lu(389) * lu(769)
         lu(794) = lu(794) - lu(390) * lu(769)
         lu(795) = lu(795) - lu(391) * lu(769)
         lu(797) = lu(797) - lu(392) * lu(769)
         lu(1052) = lu(1052) - lu(380) * lu(1045)
         lu(1055) = - lu(381) * lu(1045)
         lu(1058) = lu(1058) - lu(382) * lu(1045)
         lu(1062) = lu(1062) - lu(383) * lu(1045)
         lu(1064) = lu(1064) - lu(384) * lu(1045)
         lu(1065) = lu(1065) - lu(385) * lu(1045)
         lu(1067) = lu(1067) - lu(386) * lu(1045)
         lu(1069) = lu(1069) - lu(387) * lu(1045)
         lu(1071) = lu(1071) - lu(388) * lu(1045)
         lu(1073) = lu(1073) - lu(389) * lu(1045)
         lu(1077) = lu(1077) - lu(390) * lu(1045)
         lu(1078) = lu(1078) - lu(391) * lu(1045)
         lu(1080) = lu(1080) - lu(392) * lu(1045)
         lu(395) = 1._r8 / lu(395)
         lu(396) = lu(396) * lu(395)
         lu(397) = lu(397) * lu(395)
         lu(398) = lu(398) * lu(395)
         lu(399) = lu(399) * lu(395)
         lu(400) = lu(400) * lu(395)
         lu(401) = lu(401) * lu(395)
         lu(402) = lu(402) * lu(395)
         lu(403) = lu(403) * lu(395)
         lu(404) = lu(404) * lu(395)
         lu(405) = lu(405) * lu(395)
         lu(667) = lu(667) - lu(396) * lu(665)
         lu(668) = lu(668) - lu(397) * lu(665)
         lu(677) = lu(677) - lu(398) * lu(665)
         lu(678) = lu(678) - lu(399) * lu(665)
         lu(681) = lu(681) - lu(400) * lu(665)
         lu(682) = lu(682) - lu(401) * lu(665)
         lu(684) = lu(684) - lu(402) * lu(665)
         lu(685) = lu(685) - lu(403) * lu(665)
         lu(686) = lu(686) - lu(404) * lu(665)
         lu(687) = lu(687) - lu(405) * lu(665)
         lu(818) = lu(818) - lu(396) * lu(816)
         lu(819) = lu(819) - lu(397) * lu(816)
         lu(831) = lu(831) - lu(398) * lu(816)
         lu(832) = lu(832) - lu(399) * lu(816)
         lu(836) = lu(836) - lu(400) * lu(816)
         lu(837) = lu(837) - lu(401) * lu(816)
         lu(840) = lu(840) - lu(402) * lu(816)
         lu(841) = lu(841) - lu(403) * lu(816)
         lu(842) = lu(842) - lu(404) * lu(816)
         lu(844) = lu(844) - lu(405) * lu(816)
         lu(877) = lu(877) - lu(396) * lu(873)
         lu(878) = lu(878) - lu(397) * lu(873)
         lu(890) = lu(890) - lu(398) * lu(873)
         lu(891) = lu(891) - lu(399) * lu(873)
         lu(896) = lu(896) - lu(400) * lu(873)
         lu(897) = lu(897) - lu(401) * lu(873)
         lu(900) = lu(900) - lu(402) * lu(873)
         lu(901) = lu(901) - lu(403) * lu(873)
         lu(902) = lu(902) - lu(404) * lu(873)
         lu(904) = lu(904) - lu(405) * lu(873)
         lu(1053) = lu(1053) - lu(396) * lu(1046)
         lu(1054) = lu(1054) - lu(397) * lu(1046)
         lu(1066) = lu(1066) - lu(398) * lu(1046)
         lu(1067) = lu(1067) - lu(399) * lu(1046)
         lu(1072) = lu(1072) - lu(400) * lu(1046)
         lu(1073) = lu(1073) - lu(401) * lu(1046)
         lu(1076) = lu(1076) - lu(402) * lu(1046)
         lu(1077) = lu(1077) - lu(403) * lu(1046)
         lu(1078) = lu(1078) - lu(404) * lu(1046)
         lu(1080) = lu(1080) - lu(405) * lu(1046)
      end subroutine lu_fac09
      subroutine lu_fac10( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(406) = 1._r8 / lu(406)
         lu(407) = lu(407) * lu(406)
         lu(408) = lu(408) * lu(406)
         lu(409) = lu(409) * lu(406)
         lu(410) = lu(410) * lu(406)
         lu(438) = lu(438) - lu(407) * lu(434)
         lu(441) = lu(441) - lu(408) * lu(434)
         lu(442) = - lu(409) * lu(434)
         lu(443) = lu(443) - lu(410) * lu(434)
         lu(465) = lu(465) - lu(407) * lu(461)
         lu(467) = - lu(408) * lu(461)
         lu(468) = lu(468) - lu(409) * lu(461)
         lu(469) = lu(469) - lu(410) * lu(461)
         lu(540) = lu(540) - lu(407) * lu(537)
         lu(542) = - lu(408) * lu(537)
         lu(543) = lu(543) - lu(409) * lu(537)
         lu(544) = lu(544) - lu(410) * lu(537)
         lu(697) = - lu(407) * lu(691)
         lu(703) = - lu(408) * lu(691)
         lu(704) = lu(704) - lu(409) * lu(691)
         lu(705) = lu(705) - lu(410) * lu(691)
         lu(733) = lu(733) - lu(407) * lu(717)
         lu(738) = lu(738) - lu(408) * lu(717)
         lu(739) = lu(739) - lu(409) * lu(717)
         lu(740) = lu(740) - lu(410) * lu(717)
         lu(915) = lu(915) - lu(407) * lu(909)
         lu(922) = lu(922) - lu(408) * lu(909)
         lu(923) = lu(923) - lu(409) * lu(909)
         lu(924) = lu(924) - lu(410) * lu(909)
         lu(934) = lu(934) - lu(407) * lu(931)
         lu(941) = lu(941) - lu(408) * lu(931)
         lu(942) = lu(942) - lu(409) * lu(931)
         lu(943) = lu(943) - lu(410) * lu(931)
         lu(965) = lu(965) - lu(407) * lu(956)
         lu(972) = lu(972) - lu(408) * lu(956)
         lu(973) = lu(973) - lu(409) * lu(956)
         lu(974) = lu(974) - lu(410) * lu(956)
         lu(1069) = lu(1069) - lu(407) * lu(1047)
         lu(1076) = lu(1076) - lu(408) * lu(1047)
         lu(1077) = lu(1077) - lu(409) * lu(1047)
         lu(1078) = lu(1078) - lu(410) * lu(1047)
         lu(1106) = lu(1106) - lu(407) * lu(1102)
         lu(1113) = lu(1113) - lu(408) * lu(1102)
         lu(1114) = lu(1114) - lu(409) * lu(1102)
         lu(1115) = lu(1115) - lu(410) * lu(1102)
         lu(412) = 1._r8 / lu(412)
         lu(413) = lu(413) * lu(412)
         lu(414) = lu(414) * lu(412)
         lu(415) = lu(415) * lu(412)
         lu(416) = lu(416) * lu(412)
         lu(417) = lu(417) * lu(412)
         lu(418) = lu(418) * lu(412)
         lu(487) = lu(487) - lu(413) * lu(484)
         lu(489) = lu(489) - lu(414) * lu(484)
         lu(491) = lu(491) - lu(415) * lu(484)
         lu(492) = lu(492) - lu(416) * lu(484)
         lu(494) = lu(494) - lu(417) * lu(484)
         lu(495) = lu(495) - lu(418) * lu(484)
         lu(583) = - lu(413) * lu(580)
         lu(588) = lu(588) - lu(414) * lu(580)
         lu(590) = lu(590) - lu(415) * lu(580)
         lu(591) = lu(591) - lu(416) * lu(580)
         lu(593) = lu(593) - lu(417) * lu(580)
         lu(594) = lu(594) - lu(418) * lu(580)
         lu(601) = lu(601) - lu(413) * lu(598)
         lu(608) = lu(608) - lu(414) * lu(598)
         lu(611) = lu(611) - lu(415) * lu(598)
         lu(612) = lu(612) - lu(416) * lu(598)
         lu(614) = lu(614) - lu(417) * lu(598)
         lu(615) = lu(615) - lu(418) * lu(598)
         lu(726) = lu(726) - lu(413) * lu(718)
         lu(733) = lu(733) - lu(414) * lu(718)
         lu(736) = lu(736) - lu(415) * lu(718)
         lu(738) = lu(738) - lu(416) * lu(718)
         lu(740) = lu(740) - lu(417) * lu(718)
         lu(741) = lu(741) - lu(418) * lu(718)
         lu(826) = lu(826) - lu(413) * lu(817)
         lu(833) = lu(833) - lu(414) * lu(817)
         lu(837) = lu(837) - lu(415) * lu(817)
         lu(840) = lu(840) - lu(416) * lu(817)
         lu(842) = lu(842) - lu(417) * lu(817)
         lu(844) = lu(844) - lu(418) * lu(817)
         lu(885) = lu(885) - lu(413) * lu(874)
         lu(893) = lu(893) - lu(414) * lu(874)
         lu(897) = lu(897) - lu(415) * lu(874)
         lu(900) = lu(900) - lu(416) * lu(874)
         lu(902) = lu(902) - lu(417) * lu(874)
         lu(904) = lu(904) - lu(418) * lu(874)
         lu(1061) = lu(1061) - lu(413) * lu(1048)
         lu(1069) = lu(1069) - lu(414) * lu(1048)
         lu(1073) = lu(1073) - lu(415) * lu(1048)
         lu(1076) = lu(1076) - lu(416) * lu(1048)
         lu(1078) = lu(1078) - lu(417) * lu(1048)
         lu(1080) = lu(1080) - lu(418) * lu(1048)
         lu(420) = 1._r8 / lu(420)
         lu(421) = lu(421) * lu(420)
         lu(422) = lu(422) * lu(420)
         lu(423) = lu(423) * lu(420)
         lu(424) = lu(424) * lu(420)
         lu(425) = lu(425) * lu(420)
         lu(426) = lu(426) * lu(420)
         lu(437) = lu(437) - lu(421) * lu(435)
         lu(439) = lu(439) - lu(422) * lu(435)
         lu(440) = lu(440) - lu(423) * lu(435)
         lu(442) = lu(442) - lu(424) * lu(435)
         lu(443) = lu(443) - lu(425) * lu(435)
         lu(444) = lu(444) - lu(426) * lu(435)
         lu(696) = lu(696) - lu(421) * lu(692)
         lu(698) = - lu(422) * lu(692)
         lu(701) = lu(701) - lu(423) * lu(692)
         lu(704) = lu(704) - lu(424) * lu(692)
         lu(705) = lu(705) - lu(425) * lu(692)
         lu(706) = lu(706) - lu(426) * lu(692)
         lu(749) = lu(749) - lu(421) * lu(746)
         lu(751) = lu(751) - lu(422) * lu(746)
         lu(755) = lu(755) - lu(423) * lu(746)
         lu(758) = - lu(424) * lu(746)
         lu(759) = lu(759) - lu(425) * lu(746)
         lu(760) = lu(760) - lu(426) * lu(746)
         lu(892) = lu(892) - lu(421) * lu(875)
         lu(894) = lu(894) - lu(422) * lu(875)
         lu(898) = lu(898) - lu(423) * lu(875)
         lu(901) = lu(901) - lu(424) * lu(875)
         lu(902) = lu(902) - lu(425) * lu(875)
         lu(903) = lu(903) - lu(426) * lu(875)
         lu(914) = lu(914) - lu(421) * lu(910)
         lu(916) = lu(916) - lu(422) * lu(910)
         lu(920) = lu(920) - lu(423) * lu(910)
         lu(923) = lu(923) - lu(424) * lu(910)
         lu(924) = lu(924) - lu(425) * lu(910)
         lu(925) = lu(925) - lu(426) * lu(910)
         lu(1068) = lu(1068) - lu(421) * lu(1049)
         lu(1070) = lu(1070) - lu(422) * lu(1049)
         lu(1074) = lu(1074) - lu(423) * lu(1049)
         lu(1077) = lu(1077) - lu(424) * lu(1049)
         lu(1078) = lu(1078) - lu(425) * lu(1049)
         lu(1079) = lu(1079) - lu(426) * lu(1049)
         lu(1087) = lu(1087) - lu(421) * lu(1083)
         lu(1089) = lu(1089) - lu(422) * lu(1083)
         lu(1093) = lu(1093) - lu(423) * lu(1083)
         lu(1096) = lu(1096) - lu(424) * lu(1083)
         lu(1097) = lu(1097) - lu(425) * lu(1083)
         lu(1098) = lu(1098) - lu(426) * lu(1083)
         lu(427) = 1._r8 / lu(427)
         lu(428) = lu(428) * lu(427)
         lu(429) = lu(429) * lu(427)
         lu(430) = lu(430) * lu(427)
         lu(431) = lu(431) * lu(427)
         lu(432) = lu(432) * lu(427)
         lu(680) = - lu(428) * lu(666)
         lu(682) = lu(682) - lu(429) * lu(666)
         lu(683) = - lu(430) * lu(666)
         lu(685) = lu(685) - lu(431) * lu(666)
         lu(686) = lu(686) - lu(432) * lu(666)
         lu(699) = - lu(428) * lu(693)
         lu(700) = - lu(429) * lu(693)
         lu(701) = lu(701) - lu(430) * lu(693)
         lu(704) = lu(704) - lu(431) * lu(693)
         lu(705) = lu(705) - lu(432) * lu(693)
         lu(788) = lu(788) - lu(428) * lu(770)
         lu(790) = lu(790) - lu(429) * lu(770)
         lu(791) = lu(791) - lu(430) * lu(770)
         lu(794) = lu(794) - lu(431) * lu(770)
         lu(795) = lu(795) - lu(432) * lu(770)
         lu(895) = lu(895) - lu(428) * lu(876)
         lu(897) = lu(897) - lu(429) * lu(876)
         lu(898) = lu(898) - lu(430) * lu(876)
         lu(901) = lu(901) - lu(431) * lu(876)
         lu(902) = lu(902) - lu(432) * lu(876)
         lu(917) = lu(917) - lu(428) * lu(911)
         lu(919) = lu(919) - lu(429) * lu(911)
         lu(920) = lu(920) - lu(430) * lu(911)
         lu(923) = lu(923) - lu(431) * lu(911)
         lu(924) = lu(924) - lu(432) * lu(911)
         lu(982) = - lu(428) * lu(978)
         lu(984) = lu(984) - lu(429) * lu(978)
         lu(985) = lu(985) - lu(430) * lu(978)
         lu(988) = lu(988) - lu(431) * lu(978)
         lu(989) = lu(989) - lu(432) * lu(978)
         lu(1071) = lu(1071) - lu(428) * lu(1050)
         lu(1073) = lu(1073) - lu(429) * lu(1050)
         lu(1074) = lu(1074) - lu(430) * lu(1050)
         lu(1077) = lu(1077) - lu(431) * lu(1050)
         lu(1078) = lu(1078) - lu(432) * lu(1050)
         lu(1090) = lu(1090) - lu(428) * lu(1084)
         lu(1092) = lu(1092) - lu(429) * lu(1084)
         lu(1093) = lu(1093) - lu(430) * lu(1084)
         lu(1096) = lu(1096) - lu(431) * lu(1084)
         lu(1097) = lu(1097) - lu(432) * lu(1084)
         lu(1108) = - lu(428) * lu(1103)
         lu(1110) = lu(1110) - lu(429) * lu(1103)
         lu(1111) = lu(1111) - lu(430) * lu(1103)
         lu(1114) = lu(1114) - lu(431) * lu(1103)
         lu(1115) = lu(1115) - lu(432) * lu(1103)
         lu(436) = 1._r8 / lu(436)
         lu(437) = lu(437) * lu(436)
         lu(438) = lu(438) * lu(436)
         lu(439) = lu(439) * lu(436)
         lu(440) = lu(440) * lu(436)
         lu(441) = lu(441) * lu(436)
         lu(442) = lu(442) * lu(436)
         lu(443) = lu(443) * lu(436)
         lu(444) = lu(444) * lu(436)
         lu(696) = lu(696) - lu(437) * lu(694)
         lu(697) = lu(697) - lu(438) * lu(694)
         lu(698) = lu(698) - lu(439) * lu(694)
         lu(701) = lu(701) - lu(440) * lu(694)
         lu(703) = lu(703) - lu(441) * lu(694)
         lu(704) = lu(704) - lu(442) * lu(694)
         lu(705) = lu(705) - lu(443) * lu(694)
         lu(706) = lu(706) - lu(444) * lu(694)
         lu(749) = lu(749) - lu(437) * lu(747)
         lu(750) = - lu(438) * lu(747)
         lu(751) = lu(751) - lu(439) * lu(747)
         lu(755) = lu(755) - lu(440) * lu(747)
         lu(757) = lu(757) - lu(441) * lu(747)
         lu(758) = lu(758) - lu(442) * lu(747)
         lu(759) = lu(759) - lu(443) * lu(747)
         lu(760) = lu(760) - lu(444) * lu(747)
         lu(914) = lu(914) - lu(437) * lu(912)
         lu(915) = lu(915) - lu(438) * lu(912)
         lu(916) = lu(916) - lu(439) * lu(912)
         lu(920) = lu(920) - lu(440) * lu(912)
         lu(922) = lu(922) - lu(441) * lu(912)
         lu(923) = lu(923) - lu(442) * lu(912)
         lu(924) = lu(924) - lu(443) * lu(912)
         lu(925) = lu(925) - lu(444) * lu(912)
         lu(964) = lu(964) - lu(437) * lu(957)
         lu(965) = lu(965) - lu(438) * lu(957)
         lu(966) = lu(966) - lu(439) * lu(957)
         lu(970) = lu(970) - lu(440) * lu(957)
         lu(972) = lu(972) - lu(441) * lu(957)
         lu(973) = lu(973) - lu(442) * lu(957)
         lu(974) = lu(974) - lu(443) * lu(957)
         lu(975) = lu(975) - lu(444) * lu(957)
         lu(1068) = lu(1068) - lu(437) * lu(1051)
         lu(1069) = lu(1069) - lu(438) * lu(1051)
         lu(1070) = lu(1070) - lu(439) * lu(1051)
         lu(1074) = lu(1074) - lu(440) * lu(1051)
         lu(1076) = lu(1076) - lu(441) * lu(1051)
         lu(1077) = lu(1077) - lu(442) * lu(1051)
         lu(1078) = lu(1078) - lu(443) * lu(1051)
         lu(1079) = lu(1079) - lu(444) * lu(1051)
         lu(1087) = lu(1087) - lu(437) * lu(1085)
         lu(1088) = lu(1088) - lu(438) * lu(1085)
         lu(1089) = lu(1089) - lu(439) * lu(1085)
         lu(1093) = lu(1093) - lu(440) * lu(1085)
         lu(1095) = - lu(441) * lu(1085)
         lu(1096) = lu(1096) - lu(442) * lu(1085)
         lu(1097) = lu(1097) - lu(443) * lu(1085)
         lu(1098) = lu(1098) - lu(444) * lu(1085)
         lu(448) = 1._r8 / lu(448)
         lu(449) = lu(449) * lu(448)
         lu(450) = lu(450) * lu(448)
         lu(451) = lu(451) * lu(448)
         lu(452) = lu(452) * lu(448)
         lu(453) = lu(453) * lu(448)
         lu(454) = lu(454) * lu(448)
         lu(455) = lu(455) * lu(448)
         lu(456) = lu(456) * lu(448)
         lu(457) = lu(457) * lu(448)
         lu(458) = lu(458) * lu(448)
         lu(459) = lu(459) * lu(448)
         lu(460) = lu(460) * lu(448)
         lu(562) = lu(562) - lu(449) * lu(561)
         lu(563) = - lu(450) * lu(561)
         lu(564) = lu(564) - lu(451) * lu(561)
         lu(568) = lu(568) - lu(452) * lu(561)
         lu(569) = - lu(453) * lu(561)
         lu(570) = lu(570) - lu(454) * lu(561)
         lu(571) = - lu(455) * lu(561)
         lu(572) = lu(572) - lu(456) * lu(561)
         lu(573) = - lu(457) * lu(561)
         lu(574) = - lu(458) * lu(561)
         lu(575) = lu(575) - lu(459) * lu(561)
         lu(576) = lu(576) - lu(460) * lu(561)
         lu(720) = lu(720) - lu(449) * lu(719)
         lu(723) = lu(723) - lu(450) * lu(719)
         lu(725) = lu(725) - lu(451) * lu(719)
         lu(732) = lu(732) - lu(452) * lu(719)
         lu(733) = lu(733) - lu(453) * lu(719)
         lu(734) = lu(734) - lu(454) * lu(719)
         lu(735) = lu(735) - lu(455) * lu(719)
         lu(736) = lu(736) - lu(456) * lu(719)
         lu(738) = lu(738) - lu(457) * lu(719)
         lu(739) = lu(739) - lu(458) * lu(719)
         lu(740) = lu(740) - lu(459) * lu(719)
         lu(741) = lu(741) - lu(460) * lu(719)
         lu(772) = lu(772) - lu(449) * lu(771)
         lu(775) = - lu(450) * lu(771)
         lu(778) = lu(778) - lu(451) * lu(771)
         lu(785) = lu(785) - lu(452) * lu(771)
         lu(786) = lu(786) - lu(453) * lu(771)
         lu(788) = lu(788) - lu(454) * lu(771)
         lu(789) = lu(789) - lu(455) * lu(771)
         lu(790) = lu(790) - lu(456) * lu(771)
         lu(793) = lu(793) - lu(457) * lu(771)
         lu(794) = lu(794) - lu(458) * lu(771)
         lu(795) = lu(795) - lu(459) * lu(771)
         lu(797) = lu(797) - lu(460) * lu(771)
         lu(1053) = lu(1053) - lu(449) * lu(1052)
         lu(1057) = lu(1057) - lu(450) * lu(1052)
         lu(1060) = lu(1060) - lu(451) * lu(1052)
         lu(1067) = lu(1067) - lu(452) * lu(1052)
         lu(1069) = lu(1069) - lu(453) * lu(1052)
         lu(1071) = lu(1071) - lu(454) * lu(1052)
         lu(1072) = lu(1072) - lu(455) * lu(1052)
         lu(1073) = lu(1073) - lu(456) * lu(1052)
         lu(1076) = lu(1076) - lu(457) * lu(1052)
         lu(1077) = lu(1077) - lu(458) * lu(1052)
         lu(1078) = lu(1078) - lu(459) * lu(1052)
         lu(1080) = lu(1080) - lu(460) * lu(1052)
      end subroutine lu_fac10
      subroutine lu_fac11( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(462) = 1._r8 / lu(462)
         lu(463) = lu(463) * lu(462)
         lu(464) = lu(464) * lu(462)
         lu(465) = lu(465) * lu(462)
         lu(466) = lu(466) * lu(462)
         lu(467) = lu(467) * lu(462)
         lu(468) = lu(468) * lu(462)
         lu(469) = lu(469) * lu(462)
         lu(567) = lu(567) - lu(463) * lu(562)
         lu(568) = lu(568) - lu(464) * lu(562)
         lu(569) = lu(569) - lu(465) * lu(562)
         lu(572) = lu(572) - lu(466) * lu(562)
         lu(573) = lu(573) - lu(467) * lu(562)
         lu(574) = lu(574) - lu(468) * lu(562)
         lu(575) = lu(575) - lu(469) * lu(562)
         lu(677) = lu(677) - lu(463) * lu(667)
         lu(678) = lu(678) - lu(464) * lu(667)
         lu(679) = - lu(465) * lu(667)
         lu(682) = lu(682) - lu(466) * lu(667)
         lu(684) = lu(684) - lu(467) * lu(667)
         lu(685) = lu(685) - lu(468) * lu(667)
         lu(686) = lu(686) - lu(469) * lu(667)
         lu(731) = lu(731) - lu(463) * lu(720)
         lu(732) = lu(732) - lu(464) * lu(720)
         lu(733) = lu(733) - lu(465) * lu(720)
         lu(736) = lu(736) - lu(466) * lu(720)
         lu(738) = lu(738) - lu(467) * lu(720)
         lu(739) = lu(739) - lu(468) * lu(720)
         lu(740) = lu(740) - lu(469) * lu(720)
         lu(784) = - lu(463) * lu(772)
         lu(785) = lu(785) - lu(464) * lu(772)
         lu(786) = lu(786) - lu(465) * lu(772)
         lu(790) = lu(790) - lu(466) * lu(772)
         lu(793) = lu(793) - lu(467) * lu(772)
         lu(794) = lu(794) - lu(468) * lu(772)
         lu(795) = lu(795) - lu(469) * lu(772)
         lu(831) = lu(831) - lu(463) * lu(818)
         lu(832) = lu(832) - lu(464) * lu(818)
         lu(833) = lu(833) - lu(465) * lu(818)
         lu(837) = lu(837) - lu(466) * lu(818)
         lu(840) = lu(840) - lu(467) * lu(818)
         lu(841) = lu(841) - lu(468) * lu(818)
         lu(842) = lu(842) - lu(469) * lu(818)
         lu(890) = lu(890) - lu(463) * lu(877)
         lu(891) = lu(891) - lu(464) * lu(877)
         lu(893) = lu(893) - lu(465) * lu(877)
         lu(897) = lu(897) - lu(466) * lu(877)
         lu(900) = lu(900) - lu(467) * lu(877)
         lu(901) = lu(901) - lu(468) * lu(877)
         lu(902) = lu(902) - lu(469) * lu(877)
         lu(1066) = lu(1066) - lu(463) * lu(1053)
         lu(1067) = lu(1067) - lu(464) * lu(1053)
         lu(1069) = lu(1069) - lu(465) * lu(1053)
         lu(1073) = lu(1073) - lu(466) * lu(1053)
         lu(1076) = lu(1076) - lu(467) * lu(1053)
         lu(1077) = lu(1077) - lu(468) * lu(1053)
         lu(1078) = lu(1078) - lu(469) * lu(1053)
         lu(472) = 1._r8 / lu(472)
         lu(473) = lu(473) * lu(472)
         lu(474) = lu(474) * lu(472)
         lu(475) = lu(475) * lu(472)
         lu(476) = lu(476) * lu(472)
         lu(477) = lu(477) * lu(472)
         lu(478) = lu(478) * lu(472)
         lu(479) = lu(479) * lu(472)
         lu(480) = lu(480) * lu(472)
         lu(481) = lu(481) * lu(472)
         lu(482) = lu(482) * lu(472)
         lu(500) = - lu(473) * lu(498)
         lu(502) = - lu(474) * lu(498)
         lu(504) = lu(504) - lu(475) * lu(498)
         lu(505) = lu(505) - lu(476) * lu(498)
         lu(506) = lu(506) - lu(477) * lu(498)
         lu(507) = lu(507) - lu(478) * lu(498)
         lu(508) = lu(508) - lu(479) * lu(498)
         lu(509) = lu(509) - lu(480) * lu(498)
         lu(510) = lu(510) - lu(481) * lu(498)
         lu(511) = - lu(482) * lu(498)
         lu(669) = lu(669) - lu(473) * lu(668)
         lu(671) = lu(671) - lu(474) * lu(668)
         lu(677) = lu(677) - lu(475) * lu(668)
         lu(678) = lu(678) - lu(476) * lu(668)
         lu(681) = lu(681) - lu(477) * lu(668)
         lu(682) = lu(682) - lu(478) * lu(668)
         lu(684) = lu(684) - lu(479) * lu(668)
         lu(685) = lu(685) - lu(480) * lu(668)
         lu(686) = lu(686) - lu(481) * lu(668)
         lu(687) = lu(687) - lu(482) * lu(668)
         lu(822) = lu(822) - lu(473) * lu(819)
         lu(825) = lu(825) - lu(474) * lu(819)
         lu(831) = lu(831) - lu(475) * lu(819)
         lu(832) = lu(832) - lu(476) * lu(819)
         lu(836) = lu(836) - lu(477) * lu(819)
         lu(837) = lu(837) - lu(478) * lu(819)
         lu(840) = lu(840) - lu(479) * lu(819)
         lu(841) = lu(841) - lu(480) * lu(819)
         lu(842) = lu(842) - lu(481) * lu(819)
         lu(844) = lu(844) - lu(482) * lu(819)
         lu(881) = lu(881) - lu(473) * lu(878)
         lu(884) = lu(884) - lu(474) * lu(878)
         lu(890) = lu(890) - lu(475) * lu(878)
         lu(891) = lu(891) - lu(476) * lu(878)
         lu(896) = lu(896) - lu(477) * lu(878)
         lu(897) = lu(897) - lu(478) * lu(878)
         lu(900) = lu(900) - lu(479) * lu(878)
         lu(901) = lu(901) - lu(480) * lu(878)
         lu(902) = lu(902) - lu(481) * lu(878)
         lu(904) = lu(904) - lu(482) * lu(878)
         lu(1057) = lu(1057) - lu(473) * lu(1054)
         lu(1060) = lu(1060) - lu(474) * lu(1054)
         lu(1066) = lu(1066) - lu(475) * lu(1054)
         lu(1067) = lu(1067) - lu(476) * lu(1054)
         lu(1072) = lu(1072) - lu(477) * lu(1054)
         lu(1073) = lu(1073) - lu(478) * lu(1054)
         lu(1076) = lu(1076) - lu(479) * lu(1054)
         lu(1077) = lu(1077) - lu(480) * lu(1054)
         lu(1078) = lu(1078) - lu(481) * lu(1054)
         lu(1080) = lu(1080) - lu(482) * lu(1054)
         lu(485) = 1._r8 / lu(485)
         lu(486) = lu(486) * lu(485)
         lu(487) = lu(487) * lu(485)
         lu(488) = lu(488) * lu(485)
         lu(489) = lu(489) * lu(485)
         lu(490) = lu(490) * lu(485)
         lu(491) = lu(491) * lu(485)
         lu(492) = lu(492) * lu(485)
         lu(493) = lu(493) * lu(485)
         lu(494) = lu(494) * lu(485)
         lu(495) = lu(495) * lu(485)
         lu(724) = lu(724) - lu(486) * lu(721)
         lu(726) = lu(726) - lu(487) * lu(721)
         lu(727) = lu(727) - lu(488) * lu(721)
         lu(733) = lu(733) - lu(489) * lu(721)
         lu(735) = lu(735) - lu(490) * lu(721)
         lu(736) = lu(736) - lu(491) * lu(721)
         lu(738) = lu(738) - lu(492) * lu(721)
         lu(739) = lu(739) - lu(493) * lu(721)
         lu(740) = lu(740) - lu(494) * lu(721)
         lu(741) = lu(741) - lu(495) * lu(721)
         lu(776) = lu(776) - lu(486) * lu(773)
         lu(779) = - lu(487) * lu(773)
         lu(780) = lu(780) - lu(488) * lu(773)
         lu(786) = lu(786) - lu(489) * lu(773)
         lu(789) = lu(789) - lu(490) * lu(773)
         lu(790) = lu(790) - lu(491) * lu(773)
         lu(793) = lu(793) - lu(492) * lu(773)
         lu(794) = lu(794) - lu(493) * lu(773)
         lu(795) = lu(795) - lu(494) * lu(773)
         lu(797) = lu(797) - lu(495) * lu(773)
         lu(823) = lu(823) - lu(486) * lu(820)
         lu(826) = lu(826) - lu(487) * lu(820)
         lu(827) = lu(827) - lu(488) * lu(820)
         lu(833) = lu(833) - lu(489) * lu(820)
         lu(836) = lu(836) - lu(490) * lu(820)
         lu(837) = lu(837) - lu(491) * lu(820)
         lu(840) = lu(840) - lu(492) * lu(820)
         lu(841) = lu(841) - lu(493) * lu(820)
         lu(842) = lu(842) - lu(494) * lu(820)
         lu(844) = lu(844) - lu(495) * lu(820)
         lu(882) = lu(882) - lu(486) * lu(879)
         lu(885) = lu(885) - lu(487) * lu(879)
         lu(886) = lu(886) - lu(488) * lu(879)
         lu(893) = lu(893) - lu(489) * lu(879)
         lu(896) = lu(896) - lu(490) * lu(879)
         lu(897) = lu(897) - lu(491) * lu(879)
         lu(900) = lu(900) - lu(492) * lu(879)
         lu(901) = lu(901) - lu(493) * lu(879)
         lu(902) = lu(902) - lu(494) * lu(879)
         lu(904) = lu(904) - lu(495) * lu(879)
         lu(1058) = lu(1058) - lu(486) * lu(1055)
         lu(1061) = lu(1061) - lu(487) * lu(1055)
         lu(1062) = lu(1062) - lu(488) * lu(1055)
         lu(1069) = lu(1069) - lu(489) * lu(1055)
         lu(1072) = lu(1072) - lu(490) * lu(1055)
         lu(1073) = lu(1073) - lu(491) * lu(1055)
         lu(1076) = lu(1076) - lu(492) * lu(1055)
         lu(1077) = lu(1077) - lu(493) * lu(1055)
         lu(1078) = lu(1078) - lu(494) * lu(1055)
         lu(1080) = lu(1080) - lu(495) * lu(1055)
         lu(499) = 1._r8 / lu(499)
         lu(500) = lu(500) * lu(499)
         lu(501) = lu(501) * lu(499)
         lu(502) = lu(502) * lu(499)
         lu(503) = lu(503) * lu(499)
         lu(504) = lu(504) * lu(499)
         lu(505) = lu(505) * lu(499)
         lu(506) = lu(506) * lu(499)
         lu(507) = lu(507) * lu(499)
         lu(508) = lu(508) * lu(499)
         lu(509) = lu(509) * lu(499)
         lu(510) = lu(510) * lu(499)
         lu(511) = lu(511) * lu(499)
         lu(723) = lu(723) - lu(500) * lu(722)
         lu(724) = lu(724) - lu(501) * lu(722)
         lu(725) = lu(725) - lu(502) * lu(722)
         lu(727) = lu(727) - lu(503) * lu(722)
         lu(731) = lu(731) - lu(504) * lu(722)
         lu(732) = lu(732) - lu(505) * lu(722)
         lu(735) = lu(735) - lu(506) * lu(722)
         lu(736) = lu(736) - lu(507) * lu(722)
         lu(738) = lu(738) - lu(508) * lu(722)
         lu(739) = lu(739) - lu(509) * lu(722)
         lu(740) = lu(740) - lu(510) * lu(722)
         lu(741) = lu(741) - lu(511) * lu(722)
         lu(775) = lu(775) - lu(500) * lu(774)
         lu(776) = lu(776) - lu(501) * lu(774)
         lu(778) = lu(778) - lu(502) * lu(774)
         lu(780) = lu(780) - lu(503) * lu(774)
         lu(784) = lu(784) - lu(504) * lu(774)
         lu(785) = lu(785) - lu(505) * lu(774)
         lu(789) = lu(789) - lu(506) * lu(774)
         lu(790) = lu(790) - lu(507) * lu(774)
         lu(793) = lu(793) - lu(508) * lu(774)
         lu(794) = lu(794) - lu(509) * lu(774)
         lu(795) = lu(795) - lu(510) * lu(774)
         lu(797) = lu(797) - lu(511) * lu(774)
         lu(822) = lu(822) - lu(500) * lu(821)
         lu(823) = lu(823) - lu(501) * lu(821)
         lu(825) = lu(825) - lu(502) * lu(821)
         lu(827) = lu(827) - lu(503) * lu(821)
         lu(831) = lu(831) - lu(504) * lu(821)
         lu(832) = lu(832) - lu(505) * lu(821)
         lu(836) = lu(836) - lu(506) * lu(821)
         lu(837) = lu(837) - lu(507) * lu(821)
         lu(840) = lu(840) - lu(508) * lu(821)
         lu(841) = lu(841) - lu(509) * lu(821)
         lu(842) = lu(842) - lu(510) * lu(821)
         lu(844) = lu(844) - lu(511) * lu(821)
         lu(881) = lu(881) - lu(500) * lu(880)
         lu(882) = lu(882) - lu(501) * lu(880)
         lu(884) = lu(884) - lu(502) * lu(880)
         lu(886) = lu(886) - lu(503) * lu(880)
         lu(890) = lu(890) - lu(504) * lu(880)
         lu(891) = lu(891) - lu(505) * lu(880)
         lu(896) = lu(896) - lu(506) * lu(880)
         lu(897) = lu(897) - lu(507) * lu(880)
         lu(900) = lu(900) - lu(508) * lu(880)
         lu(901) = lu(901) - lu(509) * lu(880)
         lu(902) = lu(902) - lu(510) * lu(880)
         lu(904) = lu(904) - lu(511) * lu(880)
         lu(1057) = lu(1057) - lu(500) * lu(1056)
         lu(1058) = lu(1058) - lu(501) * lu(1056)
         lu(1060) = lu(1060) - lu(502) * lu(1056)
         lu(1062) = lu(1062) - lu(503) * lu(1056)
         lu(1066) = lu(1066) - lu(504) * lu(1056)
         lu(1067) = lu(1067) - lu(505) * lu(1056)
         lu(1072) = lu(1072) - lu(506) * lu(1056)
         lu(1073) = lu(1073) - lu(507) * lu(1056)
         lu(1076) = lu(1076) - lu(508) * lu(1056)
         lu(1077) = lu(1077) - lu(509) * lu(1056)
         lu(1078) = lu(1078) - lu(510) * lu(1056)
         lu(1080) = lu(1080) - lu(511) * lu(1056)
         lu(512) = 1._r8 / lu(512)
         lu(513) = lu(513) * lu(512)
         lu(514) = lu(514) * lu(512)
         lu(515) = lu(515) * lu(512)
         lu(516) = lu(516) * lu(512)
         lu(517) = lu(517) * lu(512)
         lu(550) = lu(550) - lu(513) * lu(549)
         lu(552) = lu(552) - lu(514) * lu(549)
         lu(556) = lu(556) - lu(515) * lu(549)
         lu(559) = lu(559) - lu(516) * lu(549)
         lu(560) = lu(560) - lu(517) * lu(549)
         lu(564) = lu(564) - lu(513) * lu(563)
         lu(567) = lu(567) - lu(514) * lu(563)
         lu(572) = lu(572) - lu(515) * lu(563)
         lu(575) = lu(575) - lu(516) * lu(563)
         lu(576) = lu(576) - lu(517) * lu(563)
         lu(582) = lu(582) - lu(513) * lu(581)
         lu(586) = lu(586) - lu(514) * lu(581)
         lu(590) = lu(590) - lu(515) * lu(581)
         lu(593) = lu(593) - lu(516) * lu(581)
         lu(594) = lu(594) - lu(517) * lu(581)
         lu(620) = - lu(513) * lu(619)
         lu(622) = lu(622) - lu(514) * lu(619)
         lu(627) = lu(627) - lu(515) * lu(619)
         lu(630) = lu(630) - lu(516) * lu(619)
         lu(631) = lu(631) - lu(517) * lu(619)
         lu(640) = lu(640) - lu(513) * lu(638)
         lu(646) = lu(646) - lu(514) * lu(638)
         lu(651) = lu(651) - lu(515) * lu(638)
         lu(654) = lu(654) - lu(516) * lu(638)
         lu(655) = lu(655) - lu(517) * lu(638)
         lu(671) = lu(671) - lu(513) * lu(669)
         lu(677) = lu(677) - lu(514) * lu(669)
         lu(682) = lu(682) - lu(515) * lu(669)
         lu(686) = lu(686) - lu(516) * lu(669)
         lu(687) = lu(687) - lu(517) * lu(669)
         lu(725) = lu(725) - lu(513) * lu(723)
         lu(731) = lu(731) - lu(514) * lu(723)
         lu(736) = lu(736) - lu(515) * lu(723)
         lu(740) = lu(740) - lu(516) * lu(723)
         lu(741) = lu(741) - lu(517) * lu(723)
         lu(778) = lu(778) - lu(513) * lu(775)
         lu(784) = lu(784) - lu(514) * lu(775)
         lu(790) = lu(790) - lu(515) * lu(775)
         lu(795) = lu(795) - lu(516) * lu(775)
         lu(797) = lu(797) - lu(517) * lu(775)
         lu(825) = lu(825) - lu(513) * lu(822)
         lu(831) = lu(831) - lu(514) * lu(822)
         lu(837) = lu(837) - lu(515) * lu(822)
         lu(842) = lu(842) - lu(516) * lu(822)
         lu(844) = lu(844) - lu(517) * lu(822)
         lu(884) = lu(884) - lu(513) * lu(881)
         lu(890) = lu(890) - lu(514) * lu(881)
         lu(897) = lu(897) - lu(515) * lu(881)
         lu(902) = lu(902) - lu(516) * lu(881)
         lu(904) = lu(904) - lu(517) * lu(881)
         lu(960) = lu(960) - lu(513) * lu(958)
         lu(962) = lu(962) - lu(514) * lu(958)
         lu(969) = lu(969) - lu(515) * lu(958)
         lu(974) = lu(974) - lu(516) * lu(958)
         lu(976) = lu(976) - lu(517) * lu(958)
         lu(1060) = lu(1060) - lu(513) * lu(1057)
         lu(1066) = lu(1066) - lu(514) * lu(1057)
         lu(1073) = lu(1073) - lu(515) * lu(1057)
         lu(1078) = lu(1078) - lu(516) * lu(1057)
         lu(1080) = lu(1080) - lu(517) * lu(1057)
      end subroutine lu_fac11
      subroutine lu_fac12( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(518) = 1._r8 / lu(518)
         lu(519) = lu(519) * lu(518)
         lu(520) = lu(520) * lu(518)
         lu(521) = lu(521) * lu(518)
         lu(522) = lu(522) * lu(518)
         lu(523) = lu(523) * lu(518)
         lu(524) = lu(524) * lu(518)
         lu(525) = lu(525) * lu(518)
         lu(526) = lu(526) * lu(518)
         lu(527) = lu(527) * lu(518)
         lu(600) = - lu(519) * lu(599)
         lu(603) = - lu(520) * lu(599)
         lu(605) = - lu(521) * lu(599)
         lu(606) = lu(606) - lu(522) * lu(599)
         lu(609) = - lu(523) * lu(599)
         lu(611) = lu(611) - lu(524) * lu(599)
         lu(613) = - lu(525) * lu(599)
         lu(614) = lu(614) - lu(526) * lu(599)
         lu(615) = lu(615) - lu(527) * lu(599)
         lu(640) = lu(640) - lu(519) * lu(639)
         lu(643) = lu(643) - lu(520) * lu(639)
         lu(645) = lu(645) - lu(521) * lu(639)
         lu(646) = lu(646) - lu(522) * lu(639)
         lu(649) = lu(649) - lu(523) * lu(639)
         lu(651) = lu(651) - lu(524) * lu(639)
         lu(653) = lu(653) - lu(525) * lu(639)
         lu(654) = lu(654) - lu(526) * lu(639)
         lu(655) = lu(655) - lu(527) * lu(639)
         lu(671) = lu(671) - lu(519) * lu(670)
         lu(674) = lu(674) - lu(520) * lu(670)
         lu(676) = lu(676) - lu(521) * lu(670)
         lu(677) = lu(677) - lu(522) * lu(670)
         lu(680) = lu(680) - lu(523) * lu(670)
         lu(682) = lu(682) - lu(524) * lu(670)
         lu(685) = lu(685) - lu(525) * lu(670)
         lu(686) = lu(686) - lu(526) * lu(670)
         lu(687) = lu(687) - lu(527) * lu(670)
         lu(725) = lu(725) - lu(519) * lu(724)
         lu(728) = lu(728) - lu(520) * lu(724)
         lu(730) = lu(730) - lu(521) * lu(724)
         lu(731) = lu(731) - lu(522) * lu(724)
         lu(734) = lu(734) - lu(523) * lu(724)
         lu(736) = lu(736) - lu(524) * lu(724)
         lu(739) = lu(739) - lu(525) * lu(724)
         lu(740) = lu(740) - lu(526) * lu(724)
         lu(741) = lu(741) - lu(527) * lu(724)
         lu(778) = lu(778) - lu(519) * lu(776)
         lu(781) = - lu(520) * lu(776)
         lu(783) = lu(783) - lu(521) * lu(776)
         lu(784) = lu(784) - lu(522) * lu(776)
         lu(788) = lu(788) - lu(523) * lu(776)
         lu(790) = lu(790) - lu(524) * lu(776)
         lu(794) = lu(794) - lu(525) * lu(776)
         lu(795) = lu(795) - lu(526) * lu(776)
         lu(797) = lu(797) - lu(527) * lu(776)
         lu(825) = lu(825) - lu(519) * lu(823)
         lu(828) = lu(828) - lu(520) * lu(823)
         lu(830) = lu(830) - lu(521) * lu(823)
         lu(831) = lu(831) - lu(522) * lu(823)
         lu(835) = lu(835) - lu(523) * lu(823)
         lu(837) = lu(837) - lu(524) * lu(823)
         lu(841) = lu(841) - lu(525) * lu(823)
         lu(842) = lu(842) - lu(526) * lu(823)
         lu(844) = lu(844) - lu(527) * lu(823)
         lu(884) = lu(884) - lu(519) * lu(882)
         lu(887) = lu(887) - lu(520) * lu(882)
         lu(889) = lu(889) - lu(521) * lu(882)
         lu(890) = lu(890) - lu(522) * lu(882)
         lu(895) = lu(895) - lu(523) * lu(882)
         lu(897) = lu(897) - lu(524) * lu(882)
         lu(901) = lu(901) - lu(525) * lu(882)
         lu(902) = lu(902) - lu(526) * lu(882)
         lu(904) = lu(904) - lu(527) * lu(882)
         lu(1060) = lu(1060) - lu(519) * lu(1058)
         lu(1063) = lu(1063) - lu(520) * lu(1058)
         lu(1065) = lu(1065) - lu(521) * lu(1058)
         lu(1066) = lu(1066) - lu(522) * lu(1058)
         lu(1071) = lu(1071) - lu(523) * lu(1058)
         lu(1073) = lu(1073) - lu(524) * lu(1058)
         lu(1077) = lu(1077) - lu(525) * lu(1058)
         lu(1078) = lu(1078) - lu(526) * lu(1058)
         lu(1080) = lu(1080) - lu(527) * lu(1058)
         lu(529) = 1._r8 / lu(529)
         lu(530) = lu(530) * lu(529)
         lu(531) = lu(531) * lu(529)
         lu(532) = lu(532) * lu(529)
         lu(533) = lu(533) * lu(529)
         lu(534) = lu(534) * lu(529)
         lu(535) = lu(535) * lu(529)
         lu(536) = lu(536) * lu(529)
         lu(699) = lu(699) - lu(530) * lu(695)
         lu(700) = lu(700) - lu(531) * lu(695)
         lu(701) = lu(701) - lu(532) * lu(695)
         lu(702) = lu(702) - lu(533) * lu(695)
         lu(704) = lu(704) - lu(534) * lu(695)
         lu(705) = lu(705) - lu(535) * lu(695)
         lu(707) = - lu(536) * lu(695)
         lu(752) = - lu(530) * lu(748)
         lu(754) = lu(754) - lu(531) * lu(748)
         lu(755) = lu(755) - lu(532) * lu(748)
         lu(756) = lu(756) - lu(533) * lu(748)
         lu(758) = lu(758) - lu(534) * lu(748)
         lu(759) = lu(759) - lu(535) * lu(748)
         lu(761) = - lu(536) * lu(748)
         lu(788) = lu(788) - lu(530) * lu(777)
         lu(790) = lu(790) - lu(531) * lu(777)
         lu(791) = lu(791) - lu(532) * lu(777)
         lu(792) = lu(792) - lu(533) * lu(777)
         lu(794) = lu(794) - lu(534) * lu(777)
         lu(795) = lu(795) - lu(535) * lu(777)
         lu(797) = lu(797) - lu(536) * lu(777)
         lu(835) = lu(835) - lu(530) * lu(824)
         lu(837) = lu(837) - lu(531) * lu(824)
         lu(838) = lu(838) - lu(532) * lu(824)
         lu(839) = lu(839) - lu(533) * lu(824)
         lu(841) = lu(841) - lu(534) * lu(824)
         lu(842) = lu(842) - lu(535) * lu(824)
         lu(844) = lu(844) - lu(536) * lu(824)
         lu(895) = lu(895) - lu(530) * lu(883)
         lu(897) = lu(897) - lu(531) * lu(883)
         lu(898) = lu(898) - lu(532) * lu(883)
         lu(899) = lu(899) - lu(533) * lu(883)
         lu(901) = lu(901) - lu(534) * lu(883)
         lu(902) = lu(902) - lu(535) * lu(883)
         lu(904) = lu(904) - lu(536) * lu(883)
         lu(917) = lu(917) - lu(530) * lu(913)
         lu(919) = lu(919) - lu(531) * lu(913)
         lu(920) = lu(920) - lu(532) * lu(913)
         lu(921) = lu(921) - lu(533) * lu(913)
         lu(923) = lu(923) - lu(534) * lu(913)
         lu(924) = lu(924) - lu(535) * lu(913)
         lu(926) = lu(926) - lu(536) * lu(913)
         lu(936) = - lu(530) * lu(932)
         lu(938) = lu(938) - lu(531) * lu(932)
         lu(939) = lu(939) - lu(532) * lu(932)
         lu(940) = lu(940) - lu(533) * lu(932)
         lu(942) = lu(942) - lu(534) * lu(932)
         lu(943) = lu(943) - lu(535) * lu(932)
         lu(945) = - lu(536) * lu(932)
         lu(967) = lu(967) - lu(530) * lu(959)
         lu(969) = lu(969) - lu(531) * lu(959)
         lu(970) = lu(970) - lu(532) * lu(959)
         lu(971) = lu(971) - lu(533) * lu(959)
         lu(973) = lu(973) - lu(534) * lu(959)
         lu(974) = lu(974) - lu(535) * lu(959)
         lu(976) = lu(976) - lu(536) * lu(959)
         lu(982) = lu(982) - lu(530) * lu(979)
         lu(984) = lu(984) - lu(531) * lu(979)
         lu(985) = lu(985) - lu(532) * lu(979)
         lu(986) = - lu(533) * lu(979)
         lu(988) = lu(988) - lu(534) * lu(979)
         lu(989) = lu(989) - lu(535) * lu(979)
         lu(991) = lu(991) - lu(536) * lu(979)
         lu(1071) = lu(1071) - lu(530) * lu(1059)
         lu(1073) = lu(1073) - lu(531) * lu(1059)
         lu(1074) = lu(1074) - lu(532) * lu(1059)
         lu(1075) = lu(1075) - lu(533) * lu(1059)
         lu(1077) = lu(1077) - lu(534) * lu(1059)
         lu(1078) = lu(1078) - lu(535) * lu(1059)
         lu(1080) = lu(1080) - lu(536) * lu(1059)
         lu(1108) = lu(1108) - lu(530) * lu(1104)
         lu(1110) = lu(1110) - lu(531) * lu(1104)
         lu(1111) = lu(1111) - lu(532) * lu(1104)
         lu(1112) = - lu(533) * lu(1104)
         lu(1114) = lu(1114) - lu(534) * lu(1104)
         lu(1115) = lu(1115) - lu(535) * lu(1104)
         lu(1117) = lu(1117) - lu(536) * lu(1104)
         lu(538) = 1._r8 / lu(538)
         lu(539) = lu(539) * lu(538)
         lu(540) = lu(540) * lu(538)
         lu(541) = lu(541) * lu(538)
         lu(542) = lu(542) * lu(538)
         lu(543) = lu(543) * lu(538)
         lu(544) = lu(544) * lu(538)
         lu(552) = lu(552) - lu(539) * lu(550)
         lu(554) = lu(554) - lu(540) * lu(550)
         lu(556) = lu(556) - lu(541) * lu(550)
         lu(557) = lu(557) - lu(542) * lu(550)
         lu(558) = lu(558) - lu(543) * lu(550)
         lu(559) = lu(559) - lu(544) * lu(550)
         lu(567) = lu(567) - lu(539) * lu(564)
         lu(569) = lu(569) - lu(540) * lu(564)
         lu(572) = lu(572) - lu(541) * lu(564)
         lu(573) = lu(573) - lu(542) * lu(564)
         lu(574) = lu(574) - lu(543) * lu(564)
         lu(575) = lu(575) - lu(544) * lu(564)
         lu(586) = lu(586) - lu(539) * lu(582)
         lu(588) = lu(588) - lu(540) * lu(582)
         lu(590) = lu(590) - lu(541) * lu(582)
         lu(591) = lu(591) - lu(542) * lu(582)
         lu(592) = - lu(543) * lu(582)
         lu(593) = lu(593) - lu(544) * lu(582)
         lu(606) = lu(606) - lu(539) * lu(600)
         lu(608) = lu(608) - lu(540) * lu(600)
         lu(611) = lu(611) - lu(541) * lu(600)
         lu(612) = lu(612) - lu(542) * lu(600)
         lu(613) = lu(613) - lu(543) * lu(600)
         lu(614) = lu(614) - lu(544) * lu(600)
         lu(622) = lu(622) - lu(539) * lu(620)
         lu(624) = lu(624) - lu(540) * lu(620)
         lu(627) = lu(627) - lu(541) * lu(620)
         lu(628) = lu(628) - lu(542) * lu(620)
         lu(629) = lu(629) - lu(543) * lu(620)
         lu(630) = lu(630) - lu(544) * lu(620)
         lu(646) = lu(646) - lu(539) * lu(640)
         lu(648) = lu(648) - lu(540) * lu(640)
         lu(651) = lu(651) - lu(541) * lu(640)
         lu(652) = lu(652) - lu(542) * lu(640)
         lu(653) = lu(653) - lu(543) * lu(640)
         lu(654) = lu(654) - lu(544) * lu(640)
         lu(677) = lu(677) - lu(539) * lu(671)
         lu(679) = lu(679) - lu(540) * lu(671)
         lu(682) = lu(682) - lu(541) * lu(671)
         lu(684) = lu(684) - lu(542) * lu(671)
         lu(685) = lu(685) - lu(543) * lu(671)
         lu(686) = lu(686) - lu(544) * lu(671)
         lu(731) = lu(731) - lu(539) * lu(725)
         lu(733) = lu(733) - lu(540) * lu(725)
         lu(736) = lu(736) - lu(541) * lu(725)
         lu(738) = lu(738) - lu(542) * lu(725)
         lu(739) = lu(739) - lu(543) * lu(725)
         lu(740) = lu(740) - lu(544) * lu(725)
         lu(784) = lu(784) - lu(539) * lu(778)
         lu(786) = lu(786) - lu(540) * lu(778)
         lu(790) = lu(790) - lu(541) * lu(778)
         lu(793) = lu(793) - lu(542) * lu(778)
         lu(794) = lu(794) - lu(543) * lu(778)
         lu(795) = lu(795) - lu(544) * lu(778)
         lu(831) = lu(831) - lu(539) * lu(825)
         lu(833) = lu(833) - lu(540) * lu(825)
         lu(837) = lu(837) - lu(541) * lu(825)
         lu(840) = lu(840) - lu(542) * lu(825)
         lu(841) = lu(841) - lu(543) * lu(825)
         lu(842) = lu(842) - lu(544) * lu(825)
         lu(890) = lu(890) - lu(539) * lu(884)
         lu(893) = lu(893) - lu(540) * lu(884)
         lu(897) = lu(897) - lu(541) * lu(884)
         lu(900) = lu(900) - lu(542) * lu(884)
         lu(901) = lu(901) - lu(543) * lu(884)
         lu(902) = lu(902) - lu(544) * lu(884)
         lu(962) = lu(962) - lu(539) * lu(960)
         lu(965) = lu(965) - lu(540) * lu(960)
         lu(969) = lu(969) - lu(541) * lu(960)
         lu(972) = lu(972) - lu(542) * lu(960)
         lu(973) = lu(973) - lu(543) * lu(960)
         lu(974) = lu(974) - lu(544) * lu(960)
         lu(1066) = lu(1066) - lu(539) * lu(1060)
         lu(1069) = lu(1069) - lu(540) * lu(1060)
         lu(1073) = lu(1073) - lu(541) * lu(1060)
         lu(1076) = lu(1076) - lu(542) * lu(1060)
         lu(1077) = lu(1077) - lu(543) * lu(1060)
         lu(1078) = lu(1078) - lu(544) * lu(1060)
         lu(551) = 1._r8 / lu(551)
         lu(552) = lu(552) * lu(551)
         lu(553) = lu(553) * lu(551)
         lu(554) = lu(554) * lu(551)
         lu(555) = lu(555) * lu(551)
         lu(556) = lu(556) * lu(551)
         lu(557) = lu(557) * lu(551)
         lu(558) = lu(558) * lu(551)
         lu(559) = lu(559) * lu(551)
         lu(560) = lu(560) * lu(551)
         lu(586) = lu(586) - lu(552) * lu(583)
         lu(587) = lu(587) - lu(553) * lu(583)
         lu(588) = lu(588) - lu(554) * lu(583)
         lu(589) = lu(589) - lu(555) * lu(583)
         lu(590) = lu(590) - lu(556) * lu(583)
         lu(591) = lu(591) - lu(557) * lu(583)
         lu(592) = lu(592) - lu(558) * lu(583)
         lu(593) = lu(593) - lu(559) * lu(583)
         lu(594) = lu(594) - lu(560) * lu(583)
         lu(606) = lu(606) - lu(552) * lu(601)
         lu(607) = lu(607) - lu(553) * lu(601)
         lu(608) = lu(608) - lu(554) * lu(601)
         lu(610) = lu(610) - lu(555) * lu(601)
         lu(611) = lu(611) - lu(556) * lu(601)
         lu(612) = lu(612) - lu(557) * lu(601)
         lu(613) = lu(613) - lu(558) * lu(601)
         lu(614) = lu(614) - lu(559) * lu(601)
         lu(615) = lu(615) - lu(560) * lu(601)
         lu(646) = lu(646) - lu(552) * lu(641)
         lu(647) = lu(647) - lu(553) * lu(641)
         lu(648) = lu(648) - lu(554) * lu(641)
         lu(650) = lu(650) - lu(555) * lu(641)
         lu(651) = lu(651) - lu(556) * lu(641)
         lu(652) = lu(652) - lu(557) * lu(641)
         lu(653) = lu(653) - lu(558) * lu(641)
         lu(654) = lu(654) - lu(559) * lu(641)
         lu(655) = lu(655) - lu(560) * lu(641)
         lu(677) = lu(677) - lu(552) * lu(672)
         lu(678) = lu(678) - lu(553) * lu(672)
         lu(679) = lu(679) - lu(554) * lu(672)
         lu(681) = lu(681) - lu(555) * lu(672)
         lu(682) = lu(682) - lu(556) * lu(672)
         lu(684) = lu(684) - lu(557) * lu(672)
         lu(685) = lu(685) - lu(558) * lu(672)
         lu(686) = lu(686) - lu(559) * lu(672)
         lu(687) = lu(687) - lu(560) * lu(672)
         lu(731) = lu(731) - lu(552) * lu(726)
         lu(732) = lu(732) - lu(553) * lu(726)
         lu(733) = lu(733) - lu(554) * lu(726)
         lu(735) = lu(735) - lu(555) * lu(726)
         lu(736) = lu(736) - lu(556) * lu(726)
         lu(738) = lu(738) - lu(557) * lu(726)
         lu(739) = lu(739) - lu(558) * lu(726)
         lu(740) = lu(740) - lu(559) * lu(726)
         lu(741) = lu(741) - lu(560) * lu(726)
         lu(784) = lu(784) - lu(552) * lu(779)
         lu(785) = lu(785) - lu(553) * lu(779)
         lu(786) = lu(786) - lu(554) * lu(779)
         lu(789) = lu(789) - lu(555) * lu(779)
         lu(790) = lu(790) - lu(556) * lu(779)
         lu(793) = lu(793) - lu(557) * lu(779)
         lu(794) = lu(794) - lu(558) * lu(779)
         lu(795) = lu(795) - lu(559) * lu(779)
         lu(797) = lu(797) - lu(560) * lu(779)
         lu(831) = lu(831) - lu(552) * lu(826)
         lu(832) = lu(832) - lu(553) * lu(826)
         lu(833) = lu(833) - lu(554) * lu(826)
         lu(836) = lu(836) - lu(555) * lu(826)
         lu(837) = lu(837) - lu(556) * lu(826)
         lu(840) = lu(840) - lu(557) * lu(826)
         lu(841) = lu(841) - lu(558) * lu(826)
         lu(842) = lu(842) - lu(559) * lu(826)
         lu(844) = lu(844) - lu(560) * lu(826)
         lu(890) = lu(890) - lu(552) * lu(885)
         lu(891) = lu(891) - lu(553) * lu(885)
         lu(893) = lu(893) - lu(554) * lu(885)
         lu(896) = lu(896) - lu(555) * lu(885)
         lu(897) = lu(897) - lu(556) * lu(885)
         lu(900) = lu(900) - lu(557) * lu(885)
         lu(901) = lu(901) - lu(558) * lu(885)
         lu(902) = lu(902) - lu(559) * lu(885)
         lu(904) = lu(904) - lu(560) * lu(885)
         lu(1066) = lu(1066) - lu(552) * lu(1061)
         lu(1067) = lu(1067) - lu(553) * lu(1061)
         lu(1069) = lu(1069) - lu(554) * lu(1061)
         lu(1072) = lu(1072) - lu(555) * lu(1061)
         lu(1073) = lu(1073) - lu(556) * lu(1061)
         lu(1076) = lu(1076) - lu(557) * lu(1061)
         lu(1077) = lu(1077) - lu(558) * lu(1061)
         lu(1078) = lu(1078) - lu(559) * lu(1061)
         lu(1080) = lu(1080) - lu(560) * lu(1061)
         lu(565) = 1._r8 / lu(565)
         lu(566) = lu(566) * lu(565)
         lu(567) = lu(567) * lu(565)
         lu(568) = lu(568) * lu(565)
         lu(569) = lu(569) * lu(565)
         lu(570) = lu(570) * lu(565)
         lu(571) = lu(571) * lu(565)
         lu(572) = lu(572) * lu(565)
         lu(573) = lu(573) * lu(565)
         lu(574) = lu(574) * lu(565)
         lu(575) = lu(575) * lu(565)
         lu(576) = lu(576) * lu(565)
         lu(603) = lu(603) - lu(566) * lu(602)
         lu(606) = lu(606) - lu(567) * lu(602)
         lu(607) = lu(607) - lu(568) * lu(602)
         lu(608) = lu(608) - lu(569) * lu(602)
         lu(609) = lu(609) - lu(570) * lu(602)
         lu(610) = lu(610) - lu(571) * lu(602)
         lu(611) = lu(611) - lu(572) * lu(602)
         lu(612) = lu(612) - lu(573) * lu(602)
         lu(613) = lu(613) - lu(574) * lu(602)
         lu(614) = lu(614) - lu(575) * lu(602)
         lu(615) = lu(615) - lu(576) * lu(602)
         lu(643) = lu(643) - lu(566) * lu(642)
         lu(646) = lu(646) - lu(567) * lu(642)
         lu(647) = lu(647) - lu(568) * lu(642)
         lu(648) = lu(648) - lu(569) * lu(642)
         lu(649) = lu(649) - lu(570) * lu(642)
         lu(650) = lu(650) - lu(571) * lu(642)
         lu(651) = lu(651) - lu(572) * lu(642)
         lu(652) = lu(652) - lu(573) * lu(642)
         lu(653) = lu(653) - lu(574) * lu(642)
         lu(654) = lu(654) - lu(575) * lu(642)
         lu(655) = lu(655) - lu(576) * lu(642)
         lu(674) = lu(674) - lu(566) * lu(673)
         lu(677) = lu(677) - lu(567) * lu(673)
         lu(678) = lu(678) - lu(568) * lu(673)
         lu(679) = lu(679) - lu(569) * lu(673)
         lu(680) = lu(680) - lu(570) * lu(673)
         lu(681) = lu(681) - lu(571) * lu(673)
         lu(682) = lu(682) - lu(572) * lu(673)
         lu(684) = lu(684) - lu(573) * lu(673)
         lu(685) = lu(685) - lu(574) * lu(673)
         lu(686) = lu(686) - lu(575) * lu(673)
         lu(687) = lu(687) - lu(576) * lu(673)
         lu(728) = lu(728) - lu(566) * lu(727)
         lu(731) = lu(731) - lu(567) * lu(727)
         lu(732) = lu(732) - lu(568) * lu(727)
         lu(733) = lu(733) - lu(569) * lu(727)
         lu(734) = lu(734) - lu(570) * lu(727)
         lu(735) = lu(735) - lu(571) * lu(727)
         lu(736) = lu(736) - lu(572) * lu(727)
         lu(738) = lu(738) - lu(573) * lu(727)
         lu(739) = lu(739) - lu(574) * lu(727)
         lu(740) = lu(740) - lu(575) * lu(727)
         lu(741) = lu(741) - lu(576) * lu(727)
         lu(781) = lu(781) - lu(566) * lu(780)
         lu(784) = lu(784) - lu(567) * lu(780)
         lu(785) = lu(785) - lu(568) * lu(780)
         lu(786) = lu(786) - lu(569) * lu(780)
         lu(788) = lu(788) - lu(570) * lu(780)
         lu(789) = lu(789) - lu(571) * lu(780)
         lu(790) = lu(790) - lu(572) * lu(780)
         lu(793) = lu(793) - lu(573) * lu(780)
         lu(794) = lu(794) - lu(574) * lu(780)
         lu(795) = lu(795) - lu(575) * lu(780)
         lu(797) = lu(797) - lu(576) * lu(780)
         lu(828) = lu(828) - lu(566) * lu(827)
         lu(831) = lu(831) - lu(567) * lu(827)
         lu(832) = lu(832) - lu(568) * lu(827)
         lu(833) = lu(833) - lu(569) * lu(827)
         lu(835) = lu(835) - lu(570) * lu(827)
         lu(836) = lu(836) - lu(571) * lu(827)
         lu(837) = lu(837) - lu(572) * lu(827)
         lu(840) = lu(840) - lu(573) * lu(827)
         lu(841) = lu(841) - lu(574) * lu(827)
         lu(842) = lu(842) - lu(575) * lu(827)
         lu(844) = lu(844) - lu(576) * lu(827)
         lu(887) = lu(887) - lu(566) * lu(886)
         lu(890) = lu(890) - lu(567) * lu(886)
         lu(891) = lu(891) - lu(568) * lu(886)
         lu(893) = lu(893) - lu(569) * lu(886)
         lu(895) = lu(895) - lu(570) * lu(886)
         lu(896) = lu(896) - lu(571) * lu(886)
         lu(897) = lu(897) - lu(572) * lu(886)
         lu(900) = lu(900) - lu(573) * lu(886)
         lu(901) = lu(901) - lu(574) * lu(886)
         lu(902) = lu(902) - lu(575) * lu(886)
         lu(904) = lu(904) - lu(576) * lu(886)
         lu(1063) = lu(1063) - lu(566) * lu(1062)
         lu(1066) = lu(1066) - lu(567) * lu(1062)
         lu(1067) = lu(1067) - lu(568) * lu(1062)
         lu(1069) = lu(1069) - lu(569) * lu(1062)
         lu(1071) = lu(1071) - lu(570) * lu(1062)
         lu(1072) = lu(1072) - lu(571) * lu(1062)
         lu(1073) = lu(1073) - lu(572) * lu(1062)
         lu(1076) = lu(1076) - lu(573) * lu(1062)
         lu(1077) = lu(1077) - lu(574) * lu(1062)
         lu(1078) = lu(1078) - lu(575) * lu(1062)
         lu(1080) = lu(1080) - lu(576) * lu(1062)
      end subroutine lu_fac12
      subroutine lu_fac13( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(584) = 1._r8 / lu(584)
         lu(585) = lu(585) * lu(584)
         lu(586) = lu(586) * lu(584)
         lu(587) = lu(587) * lu(584)
         lu(588) = lu(588) * lu(584)
         lu(589) = lu(589) * lu(584)
         lu(590) = lu(590) * lu(584)
         lu(591) = lu(591) * lu(584)
         lu(592) = lu(592) * lu(584)
         lu(593) = lu(593) * lu(584)
         lu(594) = lu(594) * lu(584)
         lu(605) = lu(605) - lu(585) * lu(603)
         lu(606) = lu(606) - lu(586) * lu(603)
         lu(607) = lu(607) - lu(587) * lu(603)
         lu(608) = lu(608) - lu(588) * lu(603)
         lu(610) = lu(610) - lu(589) * lu(603)
         lu(611) = lu(611) - lu(590) * lu(603)
         lu(612) = lu(612) - lu(591) * lu(603)
         lu(613) = lu(613) - lu(592) * lu(603)
         lu(614) = lu(614) - lu(593) * lu(603)
         lu(615) = lu(615) - lu(594) * lu(603)
         lu(645) = lu(645) - lu(585) * lu(643)
         lu(646) = lu(646) - lu(586) * lu(643)
         lu(647) = lu(647) - lu(587) * lu(643)
         lu(648) = lu(648) - lu(588) * lu(643)
         lu(650) = lu(650) - lu(589) * lu(643)
         lu(651) = lu(651) - lu(590) * lu(643)
         lu(652) = lu(652) - lu(591) * lu(643)
         lu(653) = lu(653) - lu(592) * lu(643)
         lu(654) = lu(654) - lu(593) * lu(643)
         lu(655) = lu(655) - lu(594) * lu(643)
         lu(676) = lu(676) - lu(585) * lu(674)
         lu(677) = lu(677) - lu(586) * lu(674)
         lu(678) = lu(678) - lu(587) * lu(674)
         lu(679) = lu(679) - lu(588) * lu(674)
         lu(681) = lu(681) - lu(589) * lu(674)
         lu(682) = lu(682) - lu(590) * lu(674)
         lu(684) = lu(684) - lu(591) * lu(674)
         lu(685) = lu(685) - lu(592) * lu(674)
         lu(686) = lu(686) - lu(593) * lu(674)
         lu(687) = lu(687) - lu(594) * lu(674)
         lu(730) = lu(730) - lu(585) * lu(728)
         lu(731) = lu(731) - lu(586) * lu(728)
         lu(732) = lu(732) - lu(587) * lu(728)
         lu(733) = lu(733) - lu(588) * lu(728)
         lu(735) = lu(735) - lu(589) * lu(728)
         lu(736) = lu(736) - lu(590) * lu(728)
         lu(738) = lu(738) - lu(591) * lu(728)
         lu(739) = lu(739) - lu(592) * lu(728)
         lu(740) = lu(740) - lu(593) * lu(728)
         lu(741) = lu(741) - lu(594) * lu(728)
         lu(783) = lu(783) - lu(585) * lu(781)
         lu(784) = lu(784) - lu(586) * lu(781)
         lu(785) = lu(785) - lu(587) * lu(781)
         lu(786) = lu(786) - lu(588) * lu(781)
         lu(789) = lu(789) - lu(589) * lu(781)
         lu(790) = lu(790) - lu(590) * lu(781)
         lu(793) = lu(793) - lu(591) * lu(781)
         lu(794) = lu(794) - lu(592) * lu(781)
         lu(795) = lu(795) - lu(593) * lu(781)
         lu(797) = lu(797) - lu(594) * lu(781)
         lu(830) = lu(830) - lu(585) * lu(828)
         lu(831) = lu(831) - lu(586) * lu(828)
         lu(832) = lu(832) - lu(587) * lu(828)
         lu(833) = lu(833) - lu(588) * lu(828)
         lu(836) = lu(836) - lu(589) * lu(828)
         lu(837) = lu(837) - lu(590) * lu(828)
         lu(840) = lu(840) - lu(591) * lu(828)
         lu(841) = lu(841) - lu(592) * lu(828)
         lu(842) = lu(842) - lu(593) * lu(828)
         lu(844) = lu(844) - lu(594) * lu(828)
         lu(889) = lu(889) - lu(585) * lu(887)
         lu(890) = lu(890) - lu(586) * lu(887)
         lu(891) = lu(891) - lu(587) * lu(887)
         lu(893) = lu(893) - lu(588) * lu(887)
         lu(896) = lu(896) - lu(589) * lu(887)
         lu(897) = lu(897) - lu(590) * lu(887)
         lu(900) = lu(900) - lu(591) * lu(887)
         lu(901) = lu(901) - lu(592) * lu(887)
         lu(902) = lu(902) - lu(593) * lu(887)
         lu(904) = lu(904) - lu(594) * lu(887)
         lu(1065) = lu(1065) - lu(585) * lu(1063)
         lu(1066) = lu(1066) - lu(586) * lu(1063)
         lu(1067) = lu(1067) - lu(587) * lu(1063)
         lu(1069) = lu(1069) - lu(588) * lu(1063)
         lu(1072) = lu(1072) - lu(589) * lu(1063)
         lu(1073) = lu(1073) - lu(590) * lu(1063)
         lu(1076) = lu(1076) - lu(591) * lu(1063)
         lu(1077) = lu(1077) - lu(592) * lu(1063)
         lu(1078) = lu(1078) - lu(593) * lu(1063)
         lu(1080) = lu(1080) - lu(594) * lu(1063)
         lu(604) = 1._r8 / lu(604)
         lu(605) = lu(605) * lu(604)
         lu(606) = lu(606) * lu(604)
         lu(607) = lu(607) * lu(604)
         lu(608) = lu(608) * lu(604)
         lu(609) = lu(609) * lu(604)
         lu(610) = lu(610) * lu(604)
         lu(611) = lu(611) * lu(604)
         lu(612) = lu(612) * lu(604)
         lu(613) = lu(613) * lu(604)
         lu(614) = lu(614) * lu(604)
         lu(615) = lu(615) * lu(604)
         lu(645) = lu(645) - lu(605) * lu(644)
         lu(646) = lu(646) - lu(606) * lu(644)
         lu(647) = lu(647) - lu(607) * lu(644)
         lu(648) = lu(648) - lu(608) * lu(644)
         lu(649) = lu(649) - lu(609) * lu(644)
         lu(650) = lu(650) - lu(610) * lu(644)
         lu(651) = lu(651) - lu(611) * lu(644)
         lu(652) = lu(652) - lu(612) * lu(644)
         lu(653) = lu(653) - lu(613) * lu(644)
         lu(654) = lu(654) - lu(614) * lu(644)
         lu(655) = lu(655) - lu(615) * lu(644)
         lu(676) = lu(676) - lu(605) * lu(675)
         lu(677) = lu(677) - lu(606) * lu(675)
         lu(678) = lu(678) - lu(607) * lu(675)
         lu(679) = lu(679) - lu(608) * lu(675)
         lu(680) = lu(680) - lu(609) * lu(675)
         lu(681) = lu(681) - lu(610) * lu(675)
         lu(682) = lu(682) - lu(611) * lu(675)
         lu(684) = lu(684) - lu(612) * lu(675)
         lu(685) = lu(685) - lu(613) * lu(675)
         lu(686) = lu(686) - lu(614) * lu(675)
         lu(687) = lu(687) - lu(615) * lu(675)
         lu(730) = lu(730) - lu(605) * lu(729)
         lu(731) = lu(731) - lu(606) * lu(729)
         lu(732) = lu(732) - lu(607) * lu(729)
         lu(733) = lu(733) - lu(608) * lu(729)
         lu(734) = lu(734) - lu(609) * lu(729)
         lu(735) = lu(735) - lu(610) * lu(729)
         lu(736) = lu(736) - lu(611) * lu(729)
         lu(738) = lu(738) - lu(612) * lu(729)
         lu(739) = lu(739) - lu(613) * lu(729)
         lu(740) = lu(740) - lu(614) * lu(729)
         lu(741) = lu(741) - lu(615) * lu(729)
         lu(783) = lu(783) - lu(605) * lu(782)
         lu(784) = lu(784) - lu(606) * lu(782)
         lu(785) = lu(785) - lu(607) * lu(782)
         lu(786) = lu(786) - lu(608) * lu(782)
         lu(788) = lu(788) - lu(609) * lu(782)
         lu(789) = lu(789) - lu(610) * lu(782)
         lu(790) = lu(790) - lu(611) * lu(782)
         lu(793) = lu(793) - lu(612) * lu(782)
         lu(794) = lu(794) - lu(613) * lu(782)
         lu(795) = lu(795) - lu(614) * lu(782)
         lu(797) = lu(797) - lu(615) * lu(782)
         lu(830) = lu(830) - lu(605) * lu(829)
         lu(831) = lu(831) - lu(606) * lu(829)
         lu(832) = lu(832) - lu(607) * lu(829)
         lu(833) = lu(833) - lu(608) * lu(829)
         lu(835) = lu(835) - lu(609) * lu(829)
         lu(836) = lu(836) - lu(610) * lu(829)
         lu(837) = lu(837) - lu(611) * lu(829)
         lu(840) = lu(840) - lu(612) * lu(829)
         lu(841) = lu(841) - lu(613) * lu(829)
         lu(842) = lu(842) - lu(614) * lu(829)
         lu(844) = lu(844) - lu(615) * lu(829)
         lu(889) = lu(889) - lu(605) * lu(888)
         lu(890) = lu(890) - lu(606) * lu(888)
         lu(891) = lu(891) - lu(607) * lu(888)
         lu(893) = lu(893) - lu(608) * lu(888)
         lu(895) = lu(895) - lu(609) * lu(888)
         lu(896) = lu(896) - lu(610) * lu(888)
         lu(897) = lu(897) - lu(611) * lu(888)
         lu(900) = lu(900) - lu(612) * lu(888)
         lu(901) = lu(901) - lu(613) * lu(888)
         lu(902) = lu(902) - lu(614) * lu(888)
         lu(904) = lu(904) - lu(615) * lu(888)
         lu(1065) = lu(1065) - lu(605) * lu(1064)
         lu(1066) = lu(1066) - lu(606) * lu(1064)
         lu(1067) = lu(1067) - lu(607) * lu(1064)
         lu(1069) = lu(1069) - lu(608) * lu(1064)
         lu(1071) = lu(1071) - lu(609) * lu(1064)
         lu(1072) = lu(1072) - lu(610) * lu(1064)
         lu(1073) = lu(1073) - lu(611) * lu(1064)
         lu(1076) = lu(1076) - lu(612) * lu(1064)
         lu(1077) = lu(1077) - lu(613) * lu(1064)
         lu(1078) = lu(1078) - lu(614) * lu(1064)
         lu(1080) = lu(1080) - lu(615) * lu(1064)
         lu(621) = 1._r8 / lu(621)
         lu(622) = lu(622) * lu(621)
         lu(623) = lu(623) * lu(621)
         lu(624) = lu(624) * lu(621)
         lu(625) = lu(625) * lu(621)
         lu(626) = lu(626) * lu(621)
         lu(627) = lu(627) * lu(621)
         lu(628) = lu(628) * lu(621)
         lu(629) = lu(629) * lu(621)
         lu(630) = lu(630) * lu(621)
         lu(631) = lu(631) * lu(621)
         lu(646) = lu(646) - lu(622) * lu(645)
         lu(647) = lu(647) - lu(623) * lu(645)
         lu(648) = lu(648) - lu(624) * lu(645)
         lu(649) = lu(649) - lu(625) * lu(645)
         lu(650) = lu(650) - lu(626) * lu(645)
         lu(651) = lu(651) - lu(627) * lu(645)
         lu(652) = lu(652) - lu(628) * lu(645)
         lu(653) = lu(653) - lu(629) * lu(645)
         lu(654) = lu(654) - lu(630) * lu(645)
         lu(655) = lu(655) - lu(631) * lu(645)
         lu(677) = lu(677) - lu(622) * lu(676)
         lu(678) = lu(678) - lu(623) * lu(676)
         lu(679) = lu(679) - lu(624) * lu(676)
         lu(680) = lu(680) - lu(625) * lu(676)
         lu(681) = lu(681) - lu(626) * lu(676)
         lu(682) = lu(682) - lu(627) * lu(676)
         lu(684) = lu(684) - lu(628) * lu(676)
         lu(685) = lu(685) - lu(629) * lu(676)
         lu(686) = lu(686) - lu(630) * lu(676)
         lu(687) = lu(687) - lu(631) * lu(676)
         lu(731) = lu(731) - lu(622) * lu(730)
         lu(732) = lu(732) - lu(623) * lu(730)
         lu(733) = lu(733) - lu(624) * lu(730)
         lu(734) = lu(734) - lu(625) * lu(730)
         lu(735) = lu(735) - lu(626) * lu(730)
         lu(736) = lu(736) - lu(627) * lu(730)
         lu(738) = lu(738) - lu(628) * lu(730)
         lu(739) = lu(739) - lu(629) * lu(730)
         lu(740) = lu(740) - lu(630) * lu(730)
         lu(741) = lu(741) - lu(631) * lu(730)
         lu(784) = lu(784) - lu(622) * lu(783)
         lu(785) = lu(785) - lu(623) * lu(783)
         lu(786) = lu(786) - lu(624) * lu(783)
         lu(788) = lu(788) - lu(625) * lu(783)
         lu(789) = lu(789) - lu(626) * lu(783)
         lu(790) = lu(790) - lu(627) * lu(783)
         lu(793) = lu(793) - lu(628) * lu(783)
         lu(794) = lu(794) - lu(629) * lu(783)
         lu(795) = lu(795) - lu(630) * lu(783)
         lu(797) = lu(797) - lu(631) * lu(783)
         lu(831) = lu(831) - lu(622) * lu(830)
         lu(832) = lu(832) - lu(623) * lu(830)
         lu(833) = lu(833) - lu(624) * lu(830)
         lu(835) = lu(835) - lu(625) * lu(830)
         lu(836) = lu(836) - lu(626) * lu(830)
         lu(837) = lu(837) - lu(627) * lu(830)
         lu(840) = lu(840) - lu(628) * lu(830)
         lu(841) = lu(841) - lu(629) * lu(830)
         lu(842) = lu(842) - lu(630) * lu(830)
         lu(844) = lu(844) - lu(631) * lu(830)
         lu(890) = lu(890) - lu(622) * lu(889)
         lu(891) = lu(891) - lu(623) * lu(889)
         lu(893) = lu(893) - lu(624) * lu(889)
         lu(895) = lu(895) - lu(625) * lu(889)
         lu(896) = lu(896) - lu(626) * lu(889)
         lu(897) = lu(897) - lu(627) * lu(889)
         lu(900) = lu(900) - lu(628) * lu(889)
         lu(901) = lu(901) - lu(629) * lu(889)
         lu(902) = lu(902) - lu(630) * lu(889)
         lu(904) = lu(904) - lu(631) * lu(889)
         lu(962) = lu(962) - lu(622) * lu(961)
         lu(963) = lu(963) - lu(623) * lu(961)
         lu(965) = lu(965) - lu(624) * lu(961)
         lu(967) = lu(967) - lu(625) * lu(961)
         lu(968) = lu(968) - lu(626) * lu(961)
         lu(969) = lu(969) - lu(627) * lu(961)
         lu(972) = lu(972) - lu(628) * lu(961)
         lu(973) = lu(973) - lu(629) * lu(961)
         lu(974) = lu(974) - lu(630) * lu(961)
         lu(976) = lu(976) - lu(631) * lu(961)
         lu(1066) = lu(1066) - lu(622) * lu(1065)
         lu(1067) = lu(1067) - lu(623) * lu(1065)
         lu(1069) = lu(1069) - lu(624) * lu(1065)
         lu(1071) = lu(1071) - lu(625) * lu(1065)
         lu(1072) = lu(1072) - lu(626) * lu(1065)
         lu(1073) = lu(1073) - lu(627) * lu(1065)
         lu(1076) = lu(1076) - lu(628) * lu(1065)
         lu(1077) = lu(1077) - lu(629) * lu(1065)
         lu(1078) = lu(1078) - lu(630) * lu(1065)
         lu(1080) = lu(1080) - lu(631) * lu(1065)
         lu(646) = 1._r8 / lu(646)
         lu(647) = lu(647) * lu(646)
         lu(648) = lu(648) * lu(646)
         lu(649) = lu(649) * lu(646)
         lu(650) = lu(650) * lu(646)
         lu(651) = lu(651) * lu(646)
         lu(652) = lu(652) * lu(646)
         lu(653) = lu(653) * lu(646)
         lu(654) = lu(654) * lu(646)
         lu(655) = lu(655) * lu(646)
         lu(678) = lu(678) - lu(647) * lu(677)
         lu(679) = lu(679) - lu(648) * lu(677)
         lu(680) = lu(680) - lu(649) * lu(677)
         lu(681) = lu(681) - lu(650) * lu(677)
         lu(682) = lu(682) - lu(651) * lu(677)
         lu(684) = lu(684) - lu(652) * lu(677)
         lu(685) = lu(685) - lu(653) * lu(677)
         lu(686) = lu(686) - lu(654) * lu(677)
         lu(687) = lu(687) - lu(655) * lu(677)
         lu(732) = lu(732) - lu(647) * lu(731)
         lu(733) = lu(733) - lu(648) * lu(731)
         lu(734) = lu(734) - lu(649) * lu(731)
         lu(735) = lu(735) - lu(650) * lu(731)
         lu(736) = lu(736) - lu(651) * lu(731)
         lu(738) = lu(738) - lu(652) * lu(731)
         lu(739) = lu(739) - lu(653) * lu(731)
         lu(740) = lu(740) - lu(654) * lu(731)
         lu(741) = lu(741) - lu(655) * lu(731)
         lu(785) = lu(785) - lu(647) * lu(784)
         lu(786) = lu(786) - lu(648) * lu(784)
         lu(788) = lu(788) - lu(649) * lu(784)
         lu(789) = lu(789) - lu(650) * lu(784)
         lu(790) = lu(790) - lu(651) * lu(784)
         lu(793) = lu(793) - lu(652) * lu(784)
         lu(794) = lu(794) - lu(653) * lu(784)
         lu(795) = lu(795) - lu(654) * lu(784)
         lu(797) = lu(797) - lu(655) * lu(784)
         lu(832) = lu(832) - lu(647) * lu(831)
         lu(833) = lu(833) - lu(648) * lu(831)
         lu(835) = lu(835) - lu(649) * lu(831)
         lu(836) = lu(836) - lu(650) * lu(831)
         lu(837) = lu(837) - lu(651) * lu(831)
         lu(840) = lu(840) - lu(652) * lu(831)
         lu(841) = lu(841) - lu(653) * lu(831)
         lu(842) = lu(842) - lu(654) * lu(831)
         lu(844) = lu(844) - lu(655) * lu(831)
         lu(891) = lu(891) - lu(647) * lu(890)
         lu(893) = lu(893) - lu(648) * lu(890)
         lu(895) = lu(895) - lu(649) * lu(890)
         lu(896) = lu(896) - lu(650) * lu(890)
         lu(897) = lu(897) - lu(651) * lu(890)
         lu(900) = lu(900) - lu(652) * lu(890)
         lu(901) = lu(901) - lu(653) * lu(890)
         lu(902) = lu(902) - lu(654) * lu(890)
         lu(904) = lu(904) - lu(655) * lu(890)
         lu(963) = lu(963) - lu(647) * lu(962)
         lu(965) = lu(965) - lu(648) * lu(962)
         lu(967) = lu(967) - lu(649) * lu(962)
         lu(968) = lu(968) - lu(650) * lu(962)
         lu(969) = lu(969) - lu(651) * lu(962)
         lu(972) = lu(972) - lu(652) * lu(962)
         lu(973) = lu(973) - lu(653) * lu(962)
         lu(974) = lu(974) - lu(654) * lu(962)
         lu(976) = lu(976) - lu(655) * lu(962)
         lu(1067) = lu(1067) - lu(647) * lu(1066)
         lu(1069) = lu(1069) - lu(648) * lu(1066)
         lu(1071) = lu(1071) - lu(649) * lu(1066)
         lu(1072) = lu(1072) - lu(650) * lu(1066)
         lu(1073) = lu(1073) - lu(651) * lu(1066)
         lu(1076) = lu(1076) - lu(652) * lu(1066)
         lu(1077) = lu(1077) - lu(653) * lu(1066)
         lu(1078) = lu(1078) - lu(654) * lu(1066)
         lu(1080) = lu(1080) - lu(655) * lu(1066)
         lu(678) = 1._r8 / lu(678)
         lu(679) = lu(679) * lu(678)
         lu(680) = lu(680) * lu(678)
         lu(681) = lu(681) * lu(678)
         lu(682) = lu(682) * lu(678)
         lu(683) = lu(683) * lu(678)
         lu(684) = lu(684) * lu(678)
         lu(685) = lu(685) * lu(678)
         lu(686) = lu(686) * lu(678)
         lu(687) = lu(687) * lu(678)
         lu(733) = lu(733) - lu(679) * lu(732)
         lu(734) = lu(734) - lu(680) * lu(732)
         lu(735) = lu(735) - lu(681) * lu(732)
         lu(736) = lu(736) - lu(682) * lu(732)
         lu(737) = lu(737) - lu(683) * lu(732)
         lu(738) = lu(738) - lu(684) * lu(732)
         lu(739) = lu(739) - lu(685) * lu(732)
         lu(740) = lu(740) - lu(686) * lu(732)
         lu(741) = lu(741) - lu(687) * lu(732)
         lu(786) = lu(786) - lu(679) * lu(785)
         lu(788) = lu(788) - lu(680) * lu(785)
         lu(789) = lu(789) - lu(681) * lu(785)
         lu(790) = lu(790) - lu(682) * lu(785)
         lu(791) = lu(791) - lu(683) * lu(785)
         lu(793) = lu(793) - lu(684) * lu(785)
         lu(794) = lu(794) - lu(685) * lu(785)
         lu(795) = lu(795) - lu(686) * lu(785)
         lu(797) = lu(797) - lu(687) * lu(785)
         lu(833) = lu(833) - lu(679) * lu(832)
         lu(835) = lu(835) - lu(680) * lu(832)
         lu(836) = lu(836) - lu(681) * lu(832)
         lu(837) = lu(837) - lu(682) * lu(832)
         lu(838) = lu(838) - lu(683) * lu(832)
         lu(840) = lu(840) - lu(684) * lu(832)
         lu(841) = lu(841) - lu(685) * lu(832)
         lu(842) = lu(842) - lu(686) * lu(832)
         lu(844) = lu(844) - lu(687) * lu(832)
         lu(893) = lu(893) - lu(679) * lu(891)
         lu(895) = lu(895) - lu(680) * lu(891)
         lu(896) = lu(896) - lu(681) * lu(891)
         lu(897) = lu(897) - lu(682) * lu(891)
         lu(898) = lu(898) - lu(683) * lu(891)
         lu(900) = lu(900) - lu(684) * lu(891)
         lu(901) = lu(901) - lu(685) * lu(891)
         lu(902) = lu(902) - lu(686) * lu(891)
         lu(904) = lu(904) - lu(687) * lu(891)
         lu(965) = lu(965) - lu(679) * lu(963)
         lu(967) = lu(967) - lu(680) * lu(963)
         lu(968) = lu(968) - lu(681) * lu(963)
         lu(969) = lu(969) - lu(682) * lu(963)
         lu(970) = lu(970) - lu(683) * lu(963)
         lu(972) = lu(972) - lu(684) * lu(963)
         lu(973) = lu(973) - lu(685) * lu(963)
         lu(974) = lu(974) - lu(686) * lu(963)
         lu(976) = lu(976) - lu(687) * lu(963)
         lu(981) = - lu(679) * lu(980)
         lu(982) = lu(982) - lu(680) * lu(980)
         lu(983) = lu(983) - lu(681) * lu(980)
         lu(984) = lu(984) - lu(682) * lu(980)
         lu(985) = lu(985) - lu(683) * lu(980)
         lu(987) = - lu(684) * lu(980)
         lu(988) = lu(988) - lu(685) * lu(980)
         lu(989) = lu(989) - lu(686) * lu(980)
         lu(991) = lu(991) - lu(687) * lu(980)
         lu(1069) = lu(1069) - lu(679) * lu(1067)
         lu(1071) = lu(1071) - lu(680) * lu(1067)
         lu(1072) = lu(1072) - lu(681) * lu(1067)
         lu(1073) = lu(1073) - lu(682) * lu(1067)
         lu(1074) = lu(1074) - lu(683) * lu(1067)
         lu(1076) = lu(1076) - lu(684) * lu(1067)
         lu(1077) = lu(1077) - lu(685) * lu(1067)
         lu(1078) = lu(1078) - lu(686) * lu(1067)
         lu(1080) = lu(1080) - lu(687) * lu(1067)
         lu(1088) = lu(1088) - lu(679) * lu(1086)
         lu(1090) = lu(1090) - lu(680) * lu(1086)
         lu(1091) = - lu(681) * lu(1086)
         lu(1092) = lu(1092) - lu(682) * lu(1086)
         lu(1093) = lu(1093) - lu(683) * lu(1086)
         lu(1095) = lu(1095) - lu(684) * lu(1086)
         lu(1096) = lu(1096) - lu(685) * lu(1086)
         lu(1097) = lu(1097) - lu(686) * lu(1086)
         lu(1099) = lu(1099) - lu(687) * lu(1086)
      end subroutine lu_fac13
      subroutine lu_fac14( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(696) = 1._r8 / lu(696)
         lu(697) = lu(697) * lu(696)
         lu(698) = lu(698) * lu(696)
         lu(699) = lu(699) * lu(696)
         lu(700) = lu(700) * lu(696)
         lu(701) = lu(701) * lu(696)
         lu(702) = lu(702) * lu(696)
         lu(703) = lu(703) * lu(696)
         lu(704) = lu(704) * lu(696)
         lu(705) = lu(705) * lu(696)
         lu(706) = lu(706) * lu(696)
         lu(707) = lu(707) * lu(696)
         lu(750) = lu(750) - lu(697) * lu(749)
         lu(751) = lu(751) - lu(698) * lu(749)
         lu(752) = lu(752) - lu(699) * lu(749)
         lu(754) = lu(754) - lu(700) * lu(749)
         lu(755) = lu(755) - lu(701) * lu(749)
         lu(756) = lu(756) - lu(702) * lu(749)
         lu(757) = lu(757) - lu(703) * lu(749)
         lu(758) = lu(758) - lu(704) * lu(749)
         lu(759) = lu(759) - lu(705) * lu(749)
         lu(760) = lu(760) - lu(706) * lu(749)
         lu(761) = lu(761) - lu(707) * lu(749)
         lu(893) = lu(893) - lu(697) * lu(892)
         lu(894) = lu(894) - lu(698) * lu(892)
         lu(895) = lu(895) - lu(699) * lu(892)
         lu(897) = lu(897) - lu(700) * lu(892)
         lu(898) = lu(898) - lu(701) * lu(892)
         lu(899) = lu(899) - lu(702) * lu(892)
         lu(900) = lu(900) - lu(703) * lu(892)
         lu(901) = lu(901) - lu(704) * lu(892)
         lu(902) = lu(902) - lu(705) * lu(892)
         lu(903) = lu(903) - lu(706) * lu(892)
         lu(904) = lu(904) - lu(707) * lu(892)
         lu(915) = lu(915) - lu(697) * lu(914)
         lu(916) = lu(916) - lu(698) * lu(914)
         lu(917) = lu(917) - lu(699) * lu(914)
         lu(919) = lu(919) - lu(700) * lu(914)
         lu(920) = lu(920) - lu(701) * lu(914)
         lu(921) = lu(921) - lu(702) * lu(914)
         lu(922) = lu(922) - lu(703) * lu(914)
         lu(923) = lu(923) - lu(704) * lu(914)
         lu(924) = lu(924) - lu(705) * lu(914)
         lu(925) = lu(925) - lu(706) * lu(914)
         lu(926) = lu(926) - lu(707) * lu(914)
         lu(934) = lu(934) - lu(697) * lu(933)
         lu(935) = lu(935) - lu(698) * lu(933)
         lu(936) = lu(936) - lu(699) * lu(933)
         lu(938) = lu(938) - lu(700) * lu(933)
         lu(939) = lu(939) - lu(701) * lu(933)
         lu(940) = lu(940) - lu(702) * lu(933)
         lu(941) = lu(941) - lu(703) * lu(933)
         lu(942) = lu(942) - lu(704) * lu(933)
         lu(943) = lu(943) - lu(705) * lu(933)
         lu(944) = lu(944) - lu(706) * lu(933)
         lu(945) = lu(945) - lu(707) * lu(933)
         lu(965) = lu(965) - lu(697) * lu(964)
         lu(966) = lu(966) - lu(698) * lu(964)
         lu(967) = lu(967) - lu(699) * lu(964)
         lu(969) = lu(969) - lu(700) * lu(964)
         lu(970) = lu(970) - lu(701) * lu(964)
         lu(971) = lu(971) - lu(702) * lu(964)
         lu(972) = lu(972) - lu(703) * lu(964)
         lu(973) = lu(973) - lu(704) * lu(964)
         lu(974) = lu(974) - lu(705) * lu(964)
         lu(975) = lu(975) - lu(706) * lu(964)
         lu(976) = lu(976) - lu(707) * lu(964)
         lu(1069) = lu(1069) - lu(697) * lu(1068)
         lu(1070) = lu(1070) - lu(698) * lu(1068)
         lu(1071) = lu(1071) - lu(699) * lu(1068)
         lu(1073) = lu(1073) - lu(700) * lu(1068)
         lu(1074) = lu(1074) - lu(701) * lu(1068)
         lu(1075) = lu(1075) - lu(702) * lu(1068)
         lu(1076) = lu(1076) - lu(703) * lu(1068)
         lu(1077) = lu(1077) - lu(704) * lu(1068)
         lu(1078) = lu(1078) - lu(705) * lu(1068)
         lu(1079) = lu(1079) - lu(706) * lu(1068)
         lu(1080) = lu(1080) - lu(707) * lu(1068)
         lu(1088) = lu(1088) - lu(697) * lu(1087)
         lu(1089) = lu(1089) - lu(698) * lu(1087)
         lu(1090) = lu(1090) - lu(699) * lu(1087)
         lu(1092) = lu(1092) - lu(700) * lu(1087)
         lu(1093) = lu(1093) - lu(701) * lu(1087)
         lu(1094) = - lu(702) * lu(1087)
         lu(1095) = lu(1095) - lu(703) * lu(1087)
         lu(1096) = lu(1096) - lu(704) * lu(1087)
         lu(1097) = lu(1097) - lu(705) * lu(1087)
         lu(1098) = lu(1098) - lu(706) * lu(1087)
         lu(1099) = lu(1099) - lu(707) * lu(1087)
         lu(1106) = lu(1106) - lu(697) * lu(1105)
         lu(1107) = - lu(698) * lu(1105)
         lu(1108) = lu(1108) - lu(699) * lu(1105)
         lu(1110) = lu(1110) - lu(700) * lu(1105)
         lu(1111) = lu(1111) - lu(701) * lu(1105)
         lu(1112) = lu(1112) - lu(702) * lu(1105)
         lu(1113) = lu(1113) - lu(703) * lu(1105)
         lu(1114) = lu(1114) - lu(704) * lu(1105)
         lu(1115) = lu(1115) - lu(705) * lu(1105)
         lu(1116) = lu(1116) - lu(706) * lu(1105)
         lu(1117) = lu(1117) - lu(707) * lu(1105)
         lu(733) = 1._r8 / lu(733)
         lu(734) = lu(734) * lu(733)
         lu(735) = lu(735) * lu(733)
         lu(736) = lu(736) * lu(733)
         lu(737) = lu(737) * lu(733)
         lu(738) = lu(738) * lu(733)
         lu(739) = lu(739) * lu(733)
         lu(740) = lu(740) * lu(733)
         lu(741) = lu(741) * lu(733)
         lu(752) = lu(752) - lu(734) * lu(750)
         lu(753) = lu(753) - lu(735) * lu(750)
         lu(754) = lu(754) - lu(736) * lu(750)
         lu(755) = lu(755) - lu(737) * lu(750)
         lu(757) = lu(757) - lu(738) * lu(750)
         lu(758) = lu(758) - lu(739) * lu(750)
         lu(759) = lu(759) - lu(740) * lu(750)
         lu(761) = lu(761) - lu(741) * lu(750)
         lu(788) = lu(788) - lu(734) * lu(786)
         lu(789) = lu(789) - lu(735) * lu(786)
         lu(790) = lu(790) - lu(736) * lu(786)
         lu(791) = lu(791) - lu(737) * lu(786)
         lu(793) = lu(793) - lu(738) * lu(786)
         lu(794) = lu(794) - lu(739) * lu(786)
         lu(795) = lu(795) - lu(740) * lu(786)
         lu(797) = lu(797) - lu(741) * lu(786)
         lu(835) = lu(835) - lu(734) * lu(833)
         lu(836) = lu(836) - lu(735) * lu(833)
         lu(837) = lu(837) - lu(736) * lu(833)
         lu(838) = lu(838) - lu(737) * lu(833)
         lu(840) = lu(840) - lu(738) * lu(833)
         lu(841) = lu(841) - lu(739) * lu(833)
         lu(842) = lu(842) - lu(740) * lu(833)
         lu(844) = lu(844) - lu(741) * lu(833)
         lu(895) = lu(895) - lu(734) * lu(893)
         lu(896) = lu(896) - lu(735) * lu(893)
         lu(897) = lu(897) - lu(736) * lu(893)
         lu(898) = lu(898) - lu(737) * lu(893)
         lu(900) = lu(900) - lu(738) * lu(893)
         lu(901) = lu(901) - lu(739) * lu(893)
         lu(902) = lu(902) - lu(740) * lu(893)
         lu(904) = lu(904) - lu(741) * lu(893)
         lu(917) = lu(917) - lu(734) * lu(915)
         lu(918) = lu(918) - lu(735) * lu(915)
         lu(919) = lu(919) - lu(736) * lu(915)
         lu(920) = lu(920) - lu(737) * lu(915)
         lu(922) = lu(922) - lu(738) * lu(915)
         lu(923) = lu(923) - lu(739) * lu(915)
         lu(924) = lu(924) - lu(740) * lu(915)
         lu(926) = lu(926) - lu(741) * lu(915)
         lu(936) = lu(936) - lu(734) * lu(934)
         lu(937) = lu(937) - lu(735) * lu(934)
         lu(938) = lu(938) - lu(736) * lu(934)
         lu(939) = lu(939) - lu(737) * lu(934)
         lu(941) = lu(941) - lu(738) * lu(934)
         lu(942) = lu(942) - lu(739) * lu(934)
         lu(943) = lu(943) - lu(740) * lu(934)
         lu(945) = lu(945) - lu(741) * lu(934)
         lu(967) = lu(967) - lu(734) * lu(965)
         lu(968) = lu(968) - lu(735) * lu(965)
         lu(969) = lu(969) - lu(736) * lu(965)
         lu(970) = lu(970) - lu(737) * lu(965)
         lu(972) = lu(972) - lu(738) * lu(965)
         lu(973) = lu(973) - lu(739) * lu(965)
         lu(974) = lu(974) - lu(740) * lu(965)
         lu(976) = lu(976) - lu(741) * lu(965)
         lu(982) = lu(982) - lu(734) * lu(981)
         lu(983) = lu(983) - lu(735) * lu(981)
         lu(984) = lu(984) - lu(736) * lu(981)
         lu(985) = lu(985) - lu(737) * lu(981)
         lu(987) = lu(987) - lu(738) * lu(981)
         lu(988) = lu(988) - lu(739) * lu(981)
         lu(989) = lu(989) - lu(740) * lu(981)
         lu(991) = lu(991) - lu(741) * lu(981)
         lu(1071) = lu(1071) - lu(734) * lu(1069)
         lu(1072) = lu(1072) - lu(735) * lu(1069)
         lu(1073) = lu(1073) - lu(736) * lu(1069)
         lu(1074) = lu(1074) - lu(737) * lu(1069)
         lu(1076) = lu(1076) - lu(738) * lu(1069)
         lu(1077) = lu(1077) - lu(739) * lu(1069)
         lu(1078) = lu(1078) - lu(740) * lu(1069)
         lu(1080) = lu(1080) - lu(741) * lu(1069)
         lu(1090) = lu(1090) - lu(734) * lu(1088)
         lu(1091) = lu(1091) - lu(735) * lu(1088)
         lu(1092) = lu(1092) - lu(736) * lu(1088)
         lu(1093) = lu(1093) - lu(737) * lu(1088)
         lu(1095) = lu(1095) - lu(738) * lu(1088)
         lu(1096) = lu(1096) - lu(739) * lu(1088)
         lu(1097) = lu(1097) - lu(740) * lu(1088)
         lu(1099) = lu(1099) - lu(741) * lu(1088)
         lu(1108) = lu(1108) - lu(734) * lu(1106)
         lu(1109) = lu(1109) - lu(735) * lu(1106)
         lu(1110) = lu(1110) - lu(736) * lu(1106)
         lu(1111) = lu(1111) - lu(737) * lu(1106)
         lu(1113) = lu(1113) - lu(738) * lu(1106)
         lu(1114) = lu(1114) - lu(739) * lu(1106)
         lu(1115) = lu(1115) - lu(740) * lu(1106)
         lu(1117) = lu(1117) - lu(741) * lu(1106)
         lu(751) = 1._r8 / lu(751)
         lu(752) = lu(752) * lu(751)
         lu(753) = lu(753) * lu(751)
         lu(754) = lu(754) * lu(751)
         lu(755) = lu(755) * lu(751)
         lu(756) = lu(756) * lu(751)
         lu(757) = lu(757) * lu(751)
         lu(758) = lu(758) * lu(751)
         lu(759) = lu(759) * lu(751)
         lu(760) = lu(760) * lu(751)
         lu(761) = lu(761) * lu(751)
         lu(788) = lu(788) - lu(752) * lu(787)
         lu(789) = lu(789) - lu(753) * lu(787)
         lu(790) = lu(790) - lu(754) * lu(787)
         lu(791) = lu(791) - lu(755) * lu(787)
         lu(792) = lu(792) - lu(756) * lu(787)
         lu(793) = lu(793) - lu(757) * lu(787)
         lu(794) = lu(794) - lu(758) * lu(787)
         lu(795) = lu(795) - lu(759) * lu(787)
         lu(796) = lu(796) - lu(760) * lu(787)
         lu(797) = lu(797) - lu(761) * lu(787)
         lu(835) = lu(835) - lu(752) * lu(834)
         lu(836) = lu(836) - lu(753) * lu(834)
         lu(837) = lu(837) - lu(754) * lu(834)
         lu(838) = lu(838) - lu(755) * lu(834)
         lu(839) = lu(839) - lu(756) * lu(834)
         lu(840) = lu(840) - lu(757) * lu(834)
         lu(841) = lu(841) - lu(758) * lu(834)
         lu(842) = lu(842) - lu(759) * lu(834)
         lu(843) = lu(843) - lu(760) * lu(834)
         lu(844) = lu(844) - lu(761) * lu(834)
         lu(895) = lu(895) - lu(752) * lu(894)
         lu(896) = lu(896) - lu(753) * lu(894)
         lu(897) = lu(897) - lu(754) * lu(894)
         lu(898) = lu(898) - lu(755) * lu(894)
         lu(899) = lu(899) - lu(756) * lu(894)
         lu(900) = lu(900) - lu(757) * lu(894)
         lu(901) = lu(901) - lu(758) * lu(894)
         lu(902) = lu(902) - lu(759) * lu(894)
         lu(903) = lu(903) - lu(760) * lu(894)
         lu(904) = lu(904) - lu(761) * lu(894)
         lu(917) = lu(917) - lu(752) * lu(916)
         lu(918) = lu(918) - lu(753) * lu(916)
         lu(919) = lu(919) - lu(754) * lu(916)
         lu(920) = lu(920) - lu(755) * lu(916)
         lu(921) = lu(921) - lu(756) * lu(916)
         lu(922) = lu(922) - lu(757) * lu(916)
         lu(923) = lu(923) - lu(758) * lu(916)
         lu(924) = lu(924) - lu(759) * lu(916)
         lu(925) = lu(925) - lu(760) * lu(916)
         lu(926) = lu(926) - lu(761) * lu(916)
         lu(936) = lu(936) - lu(752) * lu(935)
         lu(937) = lu(937) - lu(753) * lu(935)
         lu(938) = lu(938) - lu(754) * lu(935)
         lu(939) = lu(939) - lu(755) * lu(935)
         lu(940) = lu(940) - lu(756) * lu(935)
         lu(941) = lu(941) - lu(757) * lu(935)
         lu(942) = lu(942) - lu(758) * lu(935)
         lu(943) = lu(943) - lu(759) * lu(935)
         lu(944) = lu(944) - lu(760) * lu(935)
         lu(945) = lu(945) - lu(761) * lu(935)
         lu(967) = lu(967) - lu(752) * lu(966)
         lu(968) = lu(968) - lu(753) * lu(966)
         lu(969) = lu(969) - lu(754) * lu(966)
         lu(970) = lu(970) - lu(755) * lu(966)
         lu(971) = lu(971) - lu(756) * lu(966)
         lu(972) = lu(972) - lu(757) * lu(966)
         lu(973) = lu(973) - lu(758) * lu(966)
         lu(974) = lu(974) - lu(759) * lu(966)
         lu(975) = lu(975) - lu(760) * lu(966)
         lu(976) = lu(976) - lu(761) * lu(966)
         lu(1071) = lu(1071) - lu(752) * lu(1070)
         lu(1072) = lu(1072) - lu(753) * lu(1070)
         lu(1073) = lu(1073) - lu(754) * lu(1070)
         lu(1074) = lu(1074) - lu(755) * lu(1070)
         lu(1075) = lu(1075) - lu(756) * lu(1070)
         lu(1076) = lu(1076) - lu(757) * lu(1070)
         lu(1077) = lu(1077) - lu(758) * lu(1070)
         lu(1078) = lu(1078) - lu(759) * lu(1070)
         lu(1079) = lu(1079) - lu(760) * lu(1070)
         lu(1080) = lu(1080) - lu(761) * lu(1070)
         lu(1090) = lu(1090) - lu(752) * lu(1089)
         lu(1091) = lu(1091) - lu(753) * lu(1089)
         lu(1092) = lu(1092) - lu(754) * lu(1089)
         lu(1093) = lu(1093) - lu(755) * lu(1089)
         lu(1094) = lu(1094) - lu(756) * lu(1089)
         lu(1095) = lu(1095) - lu(757) * lu(1089)
         lu(1096) = lu(1096) - lu(758) * lu(1089)
         lu(1097) = lu(1097) - lu(759) * lu(1089)
         lu(1098) = lu(1098) - lu(760) * lu(1089)
         lu(1099) = lu(1099) - lu(761) * lu(1089)
         lu(1108) = lu(1108) - lu(752) * lu(1107)
         lu(1109) = lu(1109) - lu(753) * lu(1107)
         lu(1110) = lu(1110) - lu(754) * lu(1107)
         lu(1111) = lu(1111) - lu(755) * lu(1107)
         lu(1112) = lu(1112) - lu(756) * lu(1107)
         lu(1113) = lu(1113) - lu(757) * lu(1107)
         lu(1114) = lu(1114) - lu(758) * lu(1107)
         lu(1115) = lu(1115) - lu(759) * lu(1107)
         lu(1116) = lu(1116) - lu(760) * lu(1107)
         lu(1117) = lu(1117) - lu(761) * lu(1107)
         lu(788) = 1._r8 / lu(788)
         lu(789) = lu(789) * lu(788)
         lu(790) = lu(790) * lu(788)
         lu(791) = lu(791) * lu(788)
         lu(792) = lu(792) * lu(788)
         lu(793) = lu(793) * lu(788)
         lu(794) = lu(794) * lu(788)
         lu(795) = lu(795) * lu(788)
         lu(796) = lu(796) * lu(788)
         lu(797) = lu(797) * lu(788)
         lu(836) = lu(836) - lu(789) * lu(835)
         lu(837) = lu(837) - lu(790) * lu(835)
         lu(838) = lu(838) - lu(791) * lu(835)
         lu(839) = lu(839) - lu(792) * lu(835)
         lu(840) = lu(840) - lu(793) * lu(835)
         lu(841) = lu(841) - lu(794) * lu(835)
         lu(842) = lu(842) - lu(795) * lu(835)
         lu(843) = lu(843) - lu(796) * lu(835)
         lu(844) = lu(844) - lu(797) * lu(835)
         lu(896) = lu(896) - lu(789) * lu(895)
         lu(897) = lu(897) - lu(790) * lu(895)
         lu(898) = lu(898) - lu(791) * lu(895)
         lu(899) = lu(899) - lu(792) * lu(895)
         lu(900) = lu(900) - lu(793) * lu(895)
         lu(901) = lu(901) - lu(794) * lu(895)
         lu(902) = lu(902) - lu(795) * lu(895)
         lu(903) = lu(903) - lu(796) * lu(895)
         lu(904) = lu(904) - lu(797) * lu(895)
         lu(918) = lu(918) - lu(789) * lu(917)
         lu(919) = lu(919) - lu(790) * lu(917)
         lu(920) = lu(920) - lu(791) * lu(917)
         lu(921) = lu(921) - lu(792) * lu(917)
         lu(922) = lu(922) - lu(793) * lu(917)
         lu(923) = lu(923) - lu(794) * lu(917)
         lu(924) = lu(924) - lu(795) * lu(917)
         lu(925) = lu(925) - lu(796) * lu(917)
         lu(926) = lu(926) - lu(797) * lu(917)
         lu(937) = lu(937) - lu(789) * lu(936)
         lu(938) = lu(938) - lu(790) * lu(936)
         lu(939) = lu(939) - lu(791) * lu(936)
         lu(940) = lu(940) - lu(792) * lu(936)
         lu(941) = lu(941) - lu(793) * lu(936)
         lu(942) = lu(942) - lu(794) * lu(936)
         lu(943) = lu(943) - lu(795) * lu(936)
         lu(944) = lu(944) - lu(796) * lu(936)
         lu(945) = lu(945) - lu(797) * lu(936)
         lu(968) = lu(968) - lu(789) * lu(967)
         lu(969) = lu(969) - lu(790) * lu(967)
         lu(970) = lu(970) - lu(791) * lu(967)
         lu(971) = lu(971) - lu(792) * lu(967)
         lu(972) = lu(972) - lu(793) * lu(967)
         lu(973) = lu(973) - lu(794) * lu(967)
         lu(974) = lu(974) - lu(795) * lu(967)
         lu(975) = lu(975) - lu(796) * lu(967)
         lu(976) = lu(976) - lu(797) * lu(967)
         lu(983) = lu(983) - lu(789) * lu(982)
         lu(984) = lu(984) - lu(790) * lu(982)
         lu(985) = lu(985) - lu(791) * lu(982)
         lu(986) = lu(986) - lu(792) * lu(982)
         lu(987) = lu(987) - lu(793) * lu(982)
         lu(988) = lu(988) - lu(794) * lu(982)
         lu(989) = lu(989) - lu(795) * lu(982)
         lu(990) = lu(990) - lu(796) * lu(982)
         lu(991) = lu(991) - lu(797) * lu(982)
         lu(1072) = lu(1072) - lu(789) * lu(1071)
         lu(1073) = lu(1073) - lu(790) * lu(1071)
         lu(1074) = lu(1074) - lu(791) * lu(1071)
         lu(1075) = lu(1075) - lu(792) * lu(1071)
         lu(1076) = lu(1076) - lu(793) * lu(1071)
         lu(1077) = lu(1077) - lu(794) * lu(1071)
         lu(1078) = lu(1078) - lu(795) * lu(1071)
         lu(1079) = lu(1079) - lu(796) * lu(1071)
         lu(1080) = lu(1080) - lu(797) * lu(1071)
         lu(1091) = lu(1091) - lu(789) * lu(1090)
         lu(1092) = lu(1092) - lu(790) * lu(1090)
         lu(1093) = lu(1093) - lu(791) * lu(1090)
         lu(1094) = lu(1094) - lu(792) * lu(1090)
         lu(1095) = lu(1095) - lu(793) * lu(1090)
         lu(1096) = lu(1096) - lu(794) * lu(1090)
         lu(1097) = lu(1097) - lu(795) * lu(1090)
         lu(1098) = lu(1098) - lu(796) * lu(1090)
         lu(1099) = lu(1099) - lu(797) * lu(1090)
         lu(1109) = lu(1109) - lu(789) * lu(1108)
         lu(1110) = lu(1110) - lu(790) * lu(1108)
         lu(1111) = lu(1111) - lu(791) * lu(1108)
         lu(1112) = lu(1112) - lu(792) * lu(1108)
         lu(1113) = lu(1113) - lu(793) * lu(1108)
         lu(1114) = lu(1114) - lu(794) * lu(1108)
         lu(1115) = lu(1115) - lu(795) * lu(1108)
         lu(1116) = lu(1116) - lu(796) * lu(1108)
         lu(1117) = lu(1117) - lu(797) * lu(1108)
         lu(836) = 1._r8 / lu(836)
         lu(837) = lu(837) * lu(836)
         lu(838) = lu(838) * lu(836)
         lu(839) = lu(839) * lu(836)
         lu(840) = lu(840) * lu(836)
         lu(841) = lu(841) * lu(836)
         lu(842) = lu(842) * lu(836)
         lu(843) = lu(843) * lu(836)
         lu(844) = lu(844) * lu(836)
         lu(897) = lu(897) - lu(837) * lu(896)
         lu(898) = lu(898) - lu(838) * lu(896)
         lu(899) = lu(899) - lu(839) * lu(896)
         lu(900) = lu(900) - lu(840) * lu(896)
         lu(901) = lu(901) - lu(841) * lu(896)
         lu(902) = lu(902) - lu(842) * lu(896)
         lu(903) = lu(903) - lu(843) * lu(896)
         lu(904) = lu(904) - lu(844) * lu(896)
         lu(919) = lu(919) - lu(837) * lu(918)
         lu(920) = lu(920) - lu(838) * lu(918)
         lu(921) = lu(921) - lu(839) * lu(918)
         lu(922) = lu(922) - lu(840) * lu(918)
         lu(923) = lu(923) - lu(841) * lu(918)
         lu(924) = lu(924) - lu(842) * lu(918)
         lu(925) = lu(925) - lu(843) * lu(918)
         lu(926) = lu(926) - lu(844) * lu(918)
         lu(938) = lu(938) - lu(837) * lu(937)
         lu(939) = lu(939) - lu(838) * lu(937)
         lu(940) = lu(940) - lu(839) * lu(937)
         lu(941) = lu(941) - lu(840) * lu(937)
         lu(942) = lu(942) - lu(841) * lu(937)
         lu(943) = lu(943) - lu(842) * lu(937)
         lu(944) = lu(944) - lu(843) * lu(937)
         lu(945) = lu(945) - lu(844) * lu(937)
         lu(969) = lu(969) - lu(837) * lu(968)
         lu(970) = lu(970) - lu(838) * lu(968)
         lu(971) = lu(971) - lu(839) * lu(968)
         lu(972) = lu(972) - lu(840) * lu(968)
         lu(973) = lu(973) - lu(841) * lu(968)
         lu(974) = lu(974) - lu(842) * lu(968)
         lu(975) = lu(975) - lu(843) * lu(968)
         lu(976) = lu(976) - lu(844) * lu(968)
         lu(984) = lu(984) - lu(837) * lu(983)
         lu(985) = lu(985) - lu(838) * lu(983)
         lu(986) = lu(986) - lu(839) * lu(983)
         lu(987) = lu(987) - lu(840) * lu(983)
         lu(988) = lu(988) - lu(841) * lu(983)
         lu(989) = lu(989) - lu(842) * lu(983)
         lu(990) = lu(990) - lu(843) * lu(983)
         lu(991) = lu(991) - lu(844) * lu(983)
         lu(1073) = lu(1073) - lu(837) * lu(1072)
         lu(1074) = lu(1074) - lu(838) * lu(1072)
         lu(1075) = lu(1075) - lu(839) * lu(1072)
         lu(1076) = lu(1076) - lu(840) * lu(1072)
         lu(1077) = lu(1077) - lu(841) * lu(1072)
         lu(1078) = lu(1078) - lu(842) * lu(1072)
         lu(1079) = lu(1079) - lu(843) * lu(1072)
         lu(1080) = lu(1080) - lu(844) * lu(1072)
         lu(1092) = lu(1092) - lu(837) * lu(1091)
         lu(1093) = lu(1093) - lu(838) * lu(1091)
         lu(1094) = lu(1094) - lu(839) * lu(1091)
         lu(1095) = lu(1095) - lu(840) * lu(1091)
         lu(1096) = lu(1096) - lu(841) * lu(1091)
         lu(1097) = lu(1097) - lu(842) * lu(1091)
         lu(1098) = lu(1098) - lu(843) * lu(1091)
         lu(1099) = lu(1099) - lu(844) * lu(1091)
         lu(1110) = lu(1110) - lu(837) * lu(1109)
         lu(1111) = lu(1111) - lu(838) * lu(1109)
         lu(1112) = lu(1112) - lu(839) * lu(1109)
         lu(1113) = lu(1113) - lu(840) * lu(1109)
         lu(1114) = lu(1114) - lu(841) * lu(1109)
         lu(1115) = lu(1115) - lu(842) * lu(1109)
         lu(1116) = lu(1116) - lu(843) * lu(1109)
         lu(1117) = lu(1117) - lu(844) * lu(1109)
      end subroutine lu_fac14
      subroutine lu_fac15( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(897) = 1._r8 / lu(897)
         lu(898) = lu(898) * lu(897)
         lu(899) = lu(899) * lu(897)
         lu(900) = lu(900) * lu(897)
         lu(901) = lu(901) * lu(897)
         lu(902) = lu(902) * lu(897)
         lu(903) = lu(903) * lu(897)
         lu(904) = lu(904) * lu(897)
         lu(920) = lu(920) - lu(898) * lu(919)
         lu(921) = lu(921) - lu(899) * lu(919)
         lu(922) = lu(922) - lu(900) * lu(919)
         lu(923) = lu(923) - lu(901) * lu(919)
         lu(924) = lu(924) - lu(902) * lu(919)
         lu(925) = lu(925) - lu(903) * lu(919)
         lu(926) = lu(926) - lu(904) * lu(919)
         lu(939) = lu(939) - lu(898) * lu(938)
         lu(940) = lu(940) - lu(899) * lu(938)
         lu(941) = lu(941) - lu(900) * lu(938)
         lu(942) = lu(942) - lu(901) * lu(938)
         lu(943) = lu(943) - lu(902) * lu(938)
         lu(944) = lu(944) - lu(903) * lu(938)
         lu(945) = lu(945) - lu(904) * lu(938)
         lu(970) = lu(970) - lu(898) * lu(969)
         lu(971) = lu(971) - lu(899) * lu(969)
         lu(972) = lu(972) - lu(900) * lu(969)
         lu(973) = lu(973) - lu(901) * lu(969)
         lu(974) = lu(974) - lu(902) * lu(969)
         lu(975) = lu(975) - lu(903) * lu(969)
         lu(976) = lu(976) - lu(904) * lu(969)
         lu(985) = lu(985) - lu(898) * lu(984)
         lu(986) = lu(986) - lu(899) * lu(984)
         lu(987) = lu(987) - lu(900) * lu(984)
         lu(988) = lu(988) - lu(901) * lu(984)
         lu(989) = lu(989) - lu(902) * lu(984)
         lu(990) = lu(990) - lu(903) * lu(984)
         lu(991) = lu(991) - lu(904) * lu(984)
         lu(1074) = lu(1074) - lu(898) * lu(1073)
         lu(1075) = lu(1075) - lu(899) * lu(1073)
         lu(1076) = lu(1076) - lu(900) * lu(1073)
         lu(1077) = lu(1077) - lu(901) * lu(1073)
         lu(1078) = lu(1078) - lu(902) * lu(1073)
         lu(1079) = lu(1079) - lu(903) * lu(1073)
         lu(1080) = lu(1080) - lu(904) * lu(1073)
         lu(1093) = lu(1093) - lu(898) * lu(1092)
         lu(1094) = lu(1094) - lu(899) * lu(1092)
         lu(1095) = lu(1095) - lu(900) * lu(1092)
         lu(1096) = lu(1096) - lu(901) * lu(1092)
         lu(1097) = lu(1097) - lu(902) * lu(1092)
         lu(1098) = lu(1098) - lu(903) * lu(1092)
         lu(1099) = lu(1099) - lu(904) * lu(1092)
         lu(1111) = lu(1111) - lu(898) * lu(1110)
         lu(1112) = lu(1112) - lu(899) * lu(1110)
         lu(1113) = lu(1113) - lu(900) * lu(1110)
         lu(1114) = lu(1114) - lu(901) * lu(1110)
         lu(1115) = lu(1115) - lu(902) * lu(1110)
         lu(1116) = lu(1116) - lu(903) * lu(1110)
         lu(1117) = lu(1117) - lu(904) * lu(1110)
         lu(920) = 1._r8 / lu(920)
         lu(921) = lu(921) * lu(920)
         lu(922) = lu(922) * lu(920)
         lu(923) = lu(923) * lu(920)
         lu(924) = lu(924) * lu(920)
         lu(925) = lu(925) * lu(920)
         lu(926) = lu(926) * lu(920)
         lu(940) = lu(940) - lu(921) * lu(939)
         lu(941) = lu(941) - lu(922) * lu(939)
         lu(942) = lu(942) - lu(923) * lu(939)
         lu(943) = lu(943) - lu(924) * lu(939)
         lu(944) = lu(944) - lu(925) * lu(939)
         lu(945) = lu(945) - lu(926) * lu(939)
         lu(971) = lu(971) - lu(921) * lu(970)
         lu(972) = lu(972) - lu(922) * lu(970)
         lu(973) = lu(973) - lu(923) * lu(970)
         lu(974) = lu(974) - lu(924) * lu(970)
         lu(975) = lu(975) - lu(925) * lu(970)
         lu(976) = lu(976) - lu(926) * lu(970)
         lu(986) = lu(986) - lu(921) * lu(985)
         lu(987) = lu(987) - lu(922) * lu(985)
         lu(988) = lu(988) - lu(923) * lu(985)
         lu(989) = lu(989) - lu(924) * lu(985)
         lu(990) = lu(990) - lu(925) * lu(985)
         lu(991) = lu(991) - lu(926) * lu(985)
         lu(1075) = lu(1075) - lu(921) * lu(1074)
         lu(1076) = lu(1076) - lu(922) * lu(1074)
         lu(1077) = lu(1077) - lu(923) * lu(1074)
         lu(1078) = lu(1078) - lu(924) * lu(1074)
         lu(1079) = lu(1079) - lu(925) * lu(1074)
         lu(1080) = lu(1080) - lu(926) * lu(1074)
         lu(1094) = lu(1094) - lu(921) * lu(1093)
         lu(1095) = lu(1095) - lu(922) * lu(1093)
         lu(1096) = lu(1096) - lu(923) * lu(1093)
         lu(1097) = lu(1097) - lu(924) * lu(1093)
         lu(1098) = lu(1098) - lu(925) * lu(1093)
         lu(1099) = lu(1099) - lu(926) * lu(1093)
         lu(1112) = lu(1112) - lu(921) * lu(1111)
         lu(1113) = lu(1113) - lu(922) * lu(1111)
         lu(1114) = lu(1114) - lu(923) * lu(1111)
         lu(1115) = lu(1115) - lu(924) * lu(1111)
         lu(1116) = lu(1116) - lu(925) * lu(1111)
         lu(1117) = lu(1117) - lu(926) * lu(1111)
         lu(940) = 1._r8 / lu(940)
         lu(941) = lu(941) * lu(940)
         lu(942) = lu(942) * lu(940)
         lu(943) = lu(943) * lu(940)
         lu(944) = lu(944) * lu(940)
         lu(945) = lu(945) * lu(940)
         lu(972) = lu(972) - lu(941) * lu(971)
         lu(973) = lu(973) - lu(942) * lu(971)
         lu(974) = lu(974) - lu(943) * lu(971)
         lu(975) = lu(975) - lu(944) * lu(971)
         lu(976) = lu(976) - lu(945) * lu(971)
         lu(987) = lu(987) - lu(941) * lu(986)
         lu(988) = lu(988) - lu(942) * lu(986)
         lu(989) = lu(989) - lu(943) * lu(986)
         lu(990) = lu(990) - lu(944) * lu(986)
         lu(991) = lu(991) - lu(945) * lu(986)
         lu(1076) = lu(1076) - lu(941) * lu(1075)
         lu(1077) = lu(1077) - lu(942) * lu(1075)
         lu(1078) = lu(1078) - lu(943) * lu(1075)
         lu(1079) = lu(1079) - lu(944) * lu(1075)
         lu(1080) = lu(1080) - lu(945) * lu(1075)
         lu(1095) = lu(1095) - lu(941) * lu(1094)
         lu(1096) = lu(1096) - lu(942) * lu(1094)
         lu(1097) = lu(1097) - lu(943) * lu(1094)
         lu(1098) = lu(1098) - lu(944) * lu(1094)
         lu(1099) = lu(1099) - lu(945) * lu(1094)
         lu(1113) = lu(1113) - lu(941) * lu(1112)
         lu(1114) = lu(1114) - lu(942) * lu(1112)
         lu(1115) = lu(1115) - lu(943) * lu(1112)
         lu(1116) = lu(1116) - lu(944) * lu(1112)
         lu(1117) = lu(1117) - lu(945) * lu(1112)
         lu(972) = 1._r8 / lu(972)
         lu(973) = lu(973) * lu(972)
         lu(974) = lu(974) * lu(972)
         lu(975) = lu(975) * lu(972)
         lu(976) = lu(976) * lu(972)
         lu(988) = lu(988) - lu(973) * lu(987)
         lu(989) = lu(989) - lu(974) * lu(987)
         lu(990) = lu(990) - lu(975) * lu(987)
         lu(991) = lu(991) - lu(976) * lu(987)
         lu(1077) = lu(1077) - lu(973) * lu(1076)
         lu(1078) = lu(1078) - lu(974) * lu(1076)
         lu(1079) = lu(1079) - lu(975) * lu(1076)
         lu(1080) = lu(1080) - lu(976) * lu(1076)
         lu(1096) = lu(1096) - lu(973) * lu(1095)
         lu(1097) = lu(1097) - lu(974) * lu(1095)
         lu(1098) = lu(1098) - lu(975) * lu(1095)
         lu(1099) = lu(1099) - lu(976) * lu(1095)
         lu(1114) = lu(1114) - lu(973) * lu(1113)
         lu(1115) = lu(1115) - lu(974) * lu(1113)
         lu(1116) = lu(1116) - lu(975) * lu(1113)
         lu(1117) = lu(1117) - lu(976) * lu(1113)
         lu(988) = 1._r8 / lu(988)
         lu(989) = lu(989) * lu(988)
         lu(990) = lu(990) * lu(988)
         lu(991) = lu(991) * lu(988)
         lu(1078) = lu(1078) - lu(989) * lu(1077)
         lu(1079) = lu(1079) - lu(990) * lu(1077)
         lu(1080) = lu(1080) - lu(991) * lu(1077)
         lu(1097) = lu(1097) - lu(989) * lu(1096)
         lu(1098) = lu(1098) - lu(990) * lu(1096)
         lu(1099) = lu(1099) - lu(991) * lu(1096)
         lu(1115) = lu(1115) - lu(989) * lu(1114)
         lu(1116) = lu(1116) - lu(990) * lu(1114)
         lu(1117) = lu(1117) - lu(991) * lu(1114)
         lu(1078) = 1._r8 / lu(1078)
         lu(1079) = lu(1079) * lu(1078)
         lu(1080) = lu(1080) * lu(1078)
         lu(1098) = lu(1098) - lu(1079) * lu(1097)
         lu(1099) = lu(1099) - lu(1080) * lu(1097)
         lu(1116) = lu(1116) - lu(1079) * lu(1115)
         lu(1117) = lu(1117) - lu(1080) * lu(1115)
         lu(1098) = 1._r8 / lu(1098)
         lu(1099) = lu(1099) * lu(1098)
         lu(1117) = lu(1117) - lu(1099) * lu(1116)
         lu(1117) = 1._r8 / lu(1117)
      end subroutine lu_fac15
      subroutine lu_fac( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
      call lu_fac01( lu )
      call lu_fac02( lu )
      call lu_fac03( lu )
      call lu_fac04( lu )
      call lu_fac05( lu )
      call lu_fac06( lu )
      call lu_fac07( lu )
      call lu_fac08( lu )
      call lu_fac09( lu )
      call lu_fac10( lu )
      call lu_fac11( lu )
      call lu_fac12( lu )
      call lu_fac13( lu )
      call lu_fac14( lu )
      call lu_fac15( lu )
      end subroutine lu_fac
      end module mo_lu_factor
